Does backreaction enforce the averaged null energy condition in semiclassical gravity? 
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The expectation value, (Tat) , of tlie renormalized stress-energy tensor of quantum fields generically 
violates the classical, local positive energy conditions of general relativity. Nevertheless, it is possible 
that (Tab) may still satisfy some nonlocal positive energy conditions. Most prominent among these 
nonlocal conditions is the averaged null energy condition (ANEC), which states that J {Tab)k°'k'' dX > 
along any complete null geodesic, where fc" denotes the geodesic tangent, with afhne parameter 
A. If ANEC holds, then traversable wormholes cannot occur. However, although ANEC holds in 
Minkowski spacetime, it is known that ANEC can be violated in curved spacetimes if one is allowed to 
choose the spacetime and quantum state arbitrarily, without imposition of the semiclassical Einstein 
equation. Gat = S,TT{Tab)- In this paper, we investigate whether ANEC holds for self-consistent 
solutions of the semiclassical Einstein equation. We study a free, linear, massless scalar field with 
arbitrary curvature coupling in the context of perturbation theory about the flat spacetime/vacuum 
solution, and we modify the perturbed semiclassical equations by the "reduction of order" procedure 
to eliminate spurious solutions. We also restrict attention to the limit in which the lengthscales 
determined by the state and metric are much larger than the Planck length. At first order in the 
metric and state perturbations, and for pure states of the scalar field, we find that the ANEC integral 
vanishes, as it must for any positivity result to hold. For mixed states, the ANEC integral can be 
negative. However, we prove that if we average the ANEC integral transverse to the geodesic, using 
a suitable Planck scale smearing function, a strictly positive result is obtained in all cases except 
for the trivial flat spacetime/vacuum solution. Similar results hold for pure states at second order 
in perturbation theory, when we additionally specialize to the situation where incoming classical 
gravitational radiation does not dominate the first order metric perturbation. These results suggest 
— in agreement with conclusions drawn by Ford and Roman from entirely independent arguments — 
that if traversable wormholes do exist as self-consistent solutions of the semiclassical equations, they 
cannot be macroscopic but must be "Planck scale". In the course of our analysis, we investigate a 
number of more general issues relevant to doing perturbative expansions of the semiclassical equations 
off of flat spacetime, including an analysis of the nature of the semiclassical Einstein equation and 
of prescriptions for extracting physically relevant solutions. A large portion of our paper is devoted 
to the treatment of these more general issues. 
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I. INTRODUCTION AND SUMMARY 

A. Brief overview of the issues addressed here 

A characteristic feature of general relativity is that it 
provides a framework for understanding many objects 
and phenomena in which spacctimc behaves in ways that 
are qualitatively completely different from our everyday 
experience and intuition. For example, solutions of Ein- 
stein's equation could in principle exist which describe 
the creation of closed timelike curves or whose topolo- 
gies are nontrivial. Whether or not such solutions exist 
depends on the nature of the matter that inhabits space- 
time. 

For the typos of matter normally considered to be phys- 
ically realistic, all observers measure locally non-negative 
energy densities in the approximation where matter is 
treated classically. This condition that TabU°'u^ > 
for all timelike (known as the "weak energy condi- 
tion") as well as other, similar, local positive energy con- 
ditions arc sufficient to strongly constrain the space of 
solutions to Einstein's equation. In particular, macro- 
scopic traversable wormholes are forbidden when this 
condition is satisfied [1,2]. Moreover, the positivity of 
locally mcasiircd energy density plays a key role in the 
positive energy theorems [3] and the singularity theorems 
[4,5], which predict that general relativity must break 
down at the endpoint of gravitational collapse. 

However, it is well known that quantum fields violate 
all the local, pointwise energy conditions [6,7]. For ex- 
ample, the Casimir vacuum for the electromagnetic field 
between two perfectly conducting plates has a negative 
local energy density; indirect effects of this have been ob- 
served experimentally [8]. Squeezed states of light also 
violate the energy conditions [9] and also have been pro- 
duced experimentally [10]. Energy condition violations 
are also fundamental to the evaporation of black holes, 
and also to particle production in a gravitational field 
(such as that sometimes hypothesized to seed galaxy for- 
mation in the early Universe) [11]. 

These ubiquitous violations of energy conditions have 
led people to consider the possibility that the semiclassi- 
cal equations could admit solutions that are qualitatively 
very different from classical solutions, such as solutions 
with negative ADM mass or solutions in which gravita- 
tional collapse occurs without the formation of singular- 
ities. In particular, in recent years there has been con- 
siderable speculation that semiclassical solutions could 
exist which contain macroscopic traversable wormholes, 
and, perhaps, even describe the creation of closed time- 
like curves in an initially causally well behaved spacetime 
[1,12,13]. 

Are such objects allowed in semiclassical gravity? 
There are three different types of possibilities [14,15]: 

• The semiclassi- 
cal equations might forbid traversable wormholes. 



the creation of closed timelike curves, and negative 
mass objects. The space of solutions of the semi- 
classical equations would then not be very different 
qualitatively from that of the classical equations. 

• The semiclassical equations might allow such ob- 
jects, but only in such a way that they always lie 
outside the domain of validity of the semiclassi- 
cal theory, either because the curvature scales are 
Planckian somewhere in the corresponding space- 
times, or because the quantum fluctuations in the 
stress tensor are comparable to its expected value. 

• The semiclassical equations might allow such ob- 
jects in situations where the semiclassical theory is 
a good approximation and the objects are "macro- 
scopic" in size (as opposed to Planck-scale). 

In the last several years, a variety of evidence has ac- 
cumulated that indicates against the third of these possi- 
bilities, and in the direction of either the first or second. 
In particular, the following evidence has been adduced 
against the possibility of creating "time machines" via 
macroscopic, traversable wormholes: First, it has been 
argued that appropriate nonlocal energy conditions may 
hold, which prevent traversable wormholes (no less time 
machines) from being produced [14-26]; see Sec. IB be- 
low. Second, it has been argued [27] that, for a wide va- 
riety of states in fiat spacetime, whenever the expected 
value of the energy density is negative, then the fiuctua- 
tions in the stress tensor are comparable to the expected 
value [28] . This suggests that the semiclassical equations 
should not be trusted in the case of solutions which de- 
pend in a crucial way on energy condition violations, such 
as in the case of traversable wormholes. Finally, it has 
been argued that even if traversable wormholes could be 
produced, quantum field effects near a chronology hori- 
zon will result in a singular {Tab), which could prevent 
the occurrence of closed timelike curves [29-31]. 

One of the principal purposes of this paper is to 
present additional evidence that nonlocal energy condi- 
tions which are sufficiently strong to rule out the occur- 
rence of macroscopic, traversable wormholes may hold in 
semiclassical gravity. We shall investigate the validity of 
averaged null energy condition (ANEC) in perturbation 
theory off of Minkowski spacetime. The key new ingredi- 
ent in our analysis is that we will impose the semiclassical 
Einstein equation 

Gab[9cd]=^T^{fab[9cd]), (1-1) 

on the spacetime and quantum state. Although we shall 
find that ANEC can be violated even for solutions of 
Eq. (1.1), we shall show that in perturbation theory, a 
transverse smearing over several Planck lengths of the 
ANEC integral is sufficient to ensure positivity. Our re- 
sults thus suggest that violations of ANEC in semiclassi- 
cal gravity may be confined to the Planck scale, where the 
semiclassical approximation itself is suspect. In particu- 
lar, since violations of ANEC are necessary for producing 



2 



traversable wormholcs [2], at the very least, it should be 
necessary for traversable wormholes to have a "Planck 
scale structure" . 

Our analysis applies only to non-self-interacting quan- 
tum fields. Thus, although the full semiclassical theory 
we are considering is nonlinear due to the coupling to the 
classical metric, the quantum portion of the theory is lin- 
ear. It is possible that semiclassical solutions for inter- 
acting fields could be qualitatively different from semi- 
classical solutions for free fields [32]. However, we are 
not aware of any evidence which suggests that this is the 
case, provided, of course, that the energy conditions for 
the interacting fields are satisfied classically. 

An additional principal purpose of this paper is to in- 
vestigate the nature of the semiclassical Einstein equa- 
tion (1.1) and its solutions. In particular, this equation 
has a character that is very similar to the radiation reac- 
tion equation for a classical charged point particle. Equa- 
tion (1.1) contains time derivatives of order higher than 
two, and , correspondingly, there exist, in effect, "extra 
degrees of freedom" in its solutions, including so-called 
runaway solutions which grow exponentially in time. We 
build on recent work of Simon [33-36] , and discuss in de- 
tail the pathologies that arise and possible resolutions. 
Our conclusion is that in the special case of perturbation 
theory about fiat space, it is possible to resolve the diffi- 
culties by a "reduction of order" prescription, but that in 
general there are still open questions with respect to the 
extraction of physical predictions from the semiclassical 
equations. 

B. Nonlocal constraints on the stress-energy tensor 

We now briefiy discuss, as background, the status of 
nonlocal energy conditions in relativity; see Yurtsever 
[19] for a recent review. Let (M, gab) be a globally hyper- 
bolic spacetime, let be a quantum field on this space- 
time, and consider the expected stress tensor, (Tab) on 
all states of this field. Although at any given point in the 
spacetime, we may choose the state so as to make the en- 
ergy density be arbitrarily negative [19], there can exist 
nonlocal constraints on the stress tensor — i.e., quan- 
tum field theory does seem to restrict the amount and 
nature of energy condition violations. A complete un- 
derstanding of these nonlocal constraints is not yet in 
hand, and the search for such an understanding is one of 
the key, active areas of research in semiclassical gravity 
[19]. Nevertheless, the results that have been obtained 
to date [15-26,14] suggest that nonlocal constraints on 
stress tensors may play a key role in restricting the space 
of solutions of the semiclassical equations. The present 
paper will present additional evidence in this direction. 

The nonlocal constraints have the following general 
form [37]. Let ^^{x) be a tensor distribution on the 
fixed spacetime, {M,gab), such that the quantity 



£ = j d^x^r\x)T,b{x) (1.2) 

is classically positive. Denote by fmin[/''^,5cd] the min- 
imum over all quantum states of the expected value of 
the quantity (1.2). There arc now three different possi- 
bilities. First, it is possible that Smini f'^'' , 9 cd] = — oo, so 
that quantum field theory does not restrict the value of 
£. This will be the case, for example, when /"'' is pro- 
portional to a four dimensional delta function, so that 
£ depends only on the value of the stress tensor at one 
point. The second possibility is that fmin [/'*'', 5cd] is fi- 
nite and negative, so that 

y d^xV^r^xXTabix)) > £min[r\9cd] (1-3) 

for all quantum states. An interesting possibility - which 

appears worthy of further investigation is that this may 
be the case whenever /"'' is smooth and of compact sup- 
port. A specific conjecture of the form (1.3) has also been 
suggested by Yurtsever [15] (see below). 

Ford and Roman [20-26,14] have derived a number of 
results of the form (1.3) in both flat and curved spacetime 
quantum field theory, which they call "quantum inequali- 
ties" . For example, Ford showed that the flux AE of neg- 
ative energy through some surface in flat spacetime, when 
averaged over a time At, must satisfy AE > —h/ At, a 
result reminiscent of the time-energy uncertainty relation 
except for the minus sign [21]. Similar results can also 
be derived for the spatial average of energy density over 
a length AL in two dimensions [38] . More recently. Ford 
and Roman have derived constraints on the average over 
time of the energy density measured at a particular point 
by inertial observers [24] in flat spacetime, and they have 
argued that their results can be extrapolated to curved 
spacetime so as to constrain certain types of traversable 
wormhole spacetimcs to be "Planck scale" [14]. 

The third possibility with respect to the quantity 
fled] is that it vanish (or be positive), so that 

j d^x^gr\x){Tab{x)) > (1.4) 

for all quantum states. Inequalities of the form (1.4) 
are usually called "averaged energy conditions" [39] . An 
example of a constraint of this type is the well known fact 
that in Minkowski spacetime, the integral of the energy 
density over a constant time slice (i.e. the Hamiltonian) 
is a positive operator. 

A particular averaged energy condition — upon which 
much attention has been focused — is the averaged null 
energy condition (ANEC), which states that 

[ {Tab)k''k''dX>0, (1.5) 

where the integral is along any complete, achronal null 
geodesic 7, fc" denotes the geodesic tangent, and A is an 
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afRne parameter [40]. The reason that this and other 
similar conditions (with null replaced by timelikc) arc 
useful is that they dovetail nicely with the methods used 
to prove global results about spacetimes in general rel- 
ativity. Many of the standard global results that were 
originally proved to hold when pointwise energy condi- 
tions are satisfied, can be shown to also hold under the 
weaker assumption that the stress tensor satisfies ANEC. 
For example, in spacetimes in which ANEC is satisfied, 
the topological censorship theorem of Friedman, Schle- 
ich and Witt [2] rules out traversable wormholes. Under 
the same hypotheses, the Pcnrose-Sorkin-Wolgar positive 
mass theorem shows that the asymptotic mass must be 
positive [41]. Finally, the positivity of the ANEC integral 
along future complete null geodesies is sufficient to prove 
some singularity theorems [40] . 

The averaged null energy condition is therefore of con- 
siderable interest. Is it enforced by quantum field theory? 
Early investigations showed that it holds in Minkowski 
spacetimc for free scalar fields and electromagnetic fields 
[16-18], and also in generic, curved, 2D spacetimes [18]. 
However, it has been shown that it can be violated in 
generic, curved four dimensional spacetimes [18,42], even 
if the spacetime is nearly flat. 

The failure of the ANEC condition in general space- 
times does not, however, sound a death knell for the pro- 
gram of deriving global results in semiclassical gravity, 
since there are some modifications of the original ANEC 
conjecture that may give rise to nontrivial constraints on 
solutions. One idea, suggested by Yurtsever [15], is sim- 
ply to weaken the conjecture from being an inequality of 
the type (1.4) to one of the type (1.3), in analogy with 
the quantum inequalities of Roman and Ford. In other 
words, a modified ANEC conjecture would be that the 
quantity fmin [/"'', 5ab] is always finite and not — oo, when 
the distribution is chosen such that the quantity £ 
is the ANEC integral along a null geodesic. Yurtsever 
shows that if this is true, then reasonable assumptions 
about the dependence of fmin on the spacetime geometry 
lead to the conclusion that macroscopic, static wormholes 
are excluded; only Planck-scale wormholes are (possibly) 
allowed. 

In this paper, however, we shall follow a different path 
by investigating the validity of ANEC when the space- 
time and quantum state are constrained by the semi- 
classical Einstein equation (1.1), since any violations of 
ANEC occurring when this equation fails to hold would 
not be physically relevant. In order to analyze generic 
solutions to Eq. (1.1), we will be forced to resort to 
perturbation theory about the trivial solution, namely, 
Minkowski spacetime with the quantum field in the vac- 
uum state. We use the "reduction of order" procedure 
to eliminate the unphysical solutions of the perturbative 
semiclassical Einstein equation. We make the additional 
approximation that "wavelengths are large compared to 
the Planck scale", and for the portions of our analysis 
involving second order perturbations, we also will need 
to assume that incoming classical gravitational radiation 



does not dominate the metric perturbation at first order. 
In the "note in proof" section of Ref. [18], violations of 
ANEC for pure states were obtained at first order in de- 
viation from flatness. A key result of our analysis is that 
this type of counterexample is eliminated by imposing the 
semiclassical equation: When Eq. (1.1) holds, the ANEC 
integral always vanishes for pure states at first order in 
deviation from flatness. This result has the side con- 
sequence that we must go to second order perturbation 
theory in order to give a complete analysis of the posi- 
tivity properties of the ANEC integral for pure states in 
nearly flat spacetimes. 

As will be described in more detail in the next subsec- 
tion, we shall show that ANEC can be violated. However, 
a suitable transversely smeared ANEC integral is always 
non-negative in the context of our perturbation expan- 
sions. The condition that a smeared ANEC integral al- 
ways be non-negative in general spacetimes is clearly a 
much weaker condition than the usual ANEC condition. 
Nevertheless, when the width of the smearing function is 
of the order of the Planck length as it is in our analyses, 
the positivity of a smeared ANEC integral would be suf- 
ficient to derive interesting constraints on the spacetime 
geometry. For example, suppose that a spacetime con- 
tains a macroscpoic traversable wormhole. Then it must 
contain one geodesic 7 for which the ANEC integral is 
negative. However, the transversely smeared ANEC in- 
tegral centered on 7 must be positive. There are now two 
qualitatively different possibilities — either the compen- 
sating positive contribution to the smeared ANEC inte- 
gral comes from within a few Planck lengths, or it comes 
from a macroscopic distance away, corresponding to the 
tail of the smearing function. 

In the first case the stress tensor and Einstein tensor 
must vary significantly over lengthscales of the order of 
the Planck length, and therefore the spacetime presum- 
ably lies outside the domain of validity of semiclassical 
gravity. In the second case, there can be macroscopic 
regions of spacetime which ANEC is violated. The exis- 
tence of violations of ANEC of this type was suggested 
by some results of Visser [43], and in Appendix E we 
present an explicit example of an approximate self con- 
sistent solution which violates ANEC in this way. In this 
second case, however, the positivity of smeared ANEC 
would restrict the amount of violation to be incredibly 
small compared to the distant, positive mass. We argue 
in Appendix E below that such violations of ANEC would 
be far to small to allow macroscopic traversable worm- 
holes. Analogous arguments apply to spacetimes with 
negative asymptotic mass and with compactly generated 
chronology horizons. This provides evidence in favor of 
the second (or first) of the three possibilities discussed 
in Sec. I A. Consequently, if the semiclassical equations 
were to enforce the positivity of a transversely smeared 
ANEC integral in general spacetimes, with a smearing 
width of order of the Planck length, this would provide 
almost as strong a constraint on physical possibilities as 
the positivity of the ANEC integral itself. 
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Our positivity result for the transversely smeared 
ANEC integral in perturbation theory is the first nonlo- 
cal constraint on stress tensors that has been proved in a 
generic class of four-dimensional, curved spacetimes. Our 
perturbative result suggests that something similar may 
be true in general spacetimes, and consequently that the 
behavior of solutions in semiclassical gravity within 
the domain of validity of that theory — may be qualita- 
tively similar to classical solutions. 

C. Brief summary of our assumptions and results 

We consider a massless scalar field with arbitrary cou- 
pling, ^, to the scalar curvature. We wish to consider 
a one-parameter family (with parameter denoted by e) 
of spacetimes (M, and quantum states satisfy- 

ing Eq. (1.1), with the spacetime reducing to Minkowski 
spacetime and the quantum state reducing to the vacuum 
state when e = 0. It is somewhat awkward and overly re- 
strictive to attempt to describe the one-parameter family 
of quantum states as though they were states in a single, 
fixed, Hilbert space, since in a general, curved spacetime 
(M, 5ob), there is no unique preferred Hilbert space of 
possible states, and, in general, there is no unique, pre- 
ferred way of identifying the states occurring in Hilbert 
space constructions in different spacetimes. It is much 
more useful to adopt the algebraic approach, wherein 
one characterizes a state by its n-point distributions on 
spacetime. We shall adopt this philosophy here and 
write the one-parameter family of states as io{s). Wc; 
shall denote the expected stress-energy in the state u by 
{Tab)uj- In fact, since the expected stress-energy tensor 
is directly determined from only the 2-point distribution 
G(x,y) = {<^{x)^{y))u; of the quantum field, the higher 
order correlation functions will play no role in our anal- 
ysis. Thus, for the purposes of our analysis a "state" 
may be viewed as synonymous with a 2-point distribu- 
tion on spacetime satisfying the wave equation in each 
variable, as well as the positivity and Hadamard condi- 
tions (see subsection II C below). However, little harm 
would be done in most of our analysis below by pretend- 
ing that u){e) corresponds to a one-parameter family of 
density matrices, /5(e) in some fixed Hilbert space, with 

{Tab)u, = tr[p fafc]. 

It is useful to characterize the state u){e) by the behav- 
ior of its correlation functions in the asymptotic past. 
Assuming that suitable asymptotic conditions hold on 
the spacetime {M,gab{£)), and state, w(e), we may as- 
sociate with co{e) a state a;in(e) on Minkowski space- 
time which agrees with lli{s) in the asymptotic past 
under an appropriate identification of (M, gatis)) with 
Minkowski spacetime. In particular, the 2-point distri- 
bution, G{x,y;e), of uj{s) can then be characterized by 
the function Fin{x,y;e) — Giri{x,y;e) — Gin.oix,y) on 
Minkowski spacetime, where Ginix,y'-e) is the 2-point 
function of u)in{s) and Gin,o{x,y) is the 2-point function 



of the Minkowski vacuum state, Lo-m.o- For Hadamard 
states, -Fin (a;, y, &) will be a smooth bisolution of the wave 
equation in Minkowski spacetime, whose initial data at 
past null infinity, J'~ , may be viewed as the freely speci- 
fiable initial data for the state (which, however, is subject 
to the positivity constraints discussed in Section II C be- 
low). 

It is well known that the semiclassical Einstein equa- 
tion is of a "higher derivative" character than the corre- 
sponding classical equation, and that consequently there 
exist new — presumably spurious — "degrees of free- 
dom" in semiclassical gravity, closely analogous to the 
"run-away" solutions which occur for the dynamics of 
a point charge in classical electrodynamics when radi- 
ation reaction effects are included. Thus, in order to 
extract any physical predictions from the semiclassical 
equations, we need a prescription either for extracting 
the "physical solutions" to these equations or for modi- 
fying the equations themselves so that the spurious solu- 
tions no longer arise. We investigate this issue in depth 
in Sec. IV. We conclude that — at least in the context of 
perturbation theory about flat spacetime — the "reduc- 
tion of order" algorithm for modifying the equations [36] 
provides a satisfactory means of eliminating the spurious 
solutions without (significantly) sacrificing accuracy at 
"long wavelengths", i.e., in the regime where the dom- 
inant lengthscales in the solution are much larger than 
the Planck scale. The validity of ANEC is investigated in 
the context of solutions to the reduced order perturbative 
semiclassical Einstein equation. 

We expand the one-parameter family of "in" -states as 

Win(e) = Win,0 + £W^n^ + ^V'^ + 0{e^) (1.6) 

and we expand the metric as 

Qabie) = Vab + eh''^^ + e^h^2 + 0{e% (1.7) 

where rjab is the flat, Minkowski metric. The metric per- 
turbation h^^^ can be written as 

/^iV=^i^'" + A/.il\ (1.8) 

where /i^^^ satisfles the linearized Einstein equation in 
vacuum and represents incoming classical gravitational 
waves at J~ . The remaining portion A/i^^^ is deter- 
mined from the first order state perturbation uj^^ via the 
reduced order, perturbed semiclassical equations. There 
is a decomposition similar to Eq. (1.8) for the second or- 

der metric perturbation h}^f^ ■ Note that to second order, 
solutions to the reduced order, perturbative semiclassical 
equations can be characterized by their "initial data" at 
past null infinity Sf~ , consisting of the incoming gravita- 
tional radiation hf^^^^^ and h''^^ as well as the incom- 
ing, freely propagating (with respect to rjab) pieces F^^^ 
and -F/j^ ' of the two-point functions associated with uj^]^ 
and ujU (see Figure 1). 
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FIG. 1. The "scattering picture" for solutions of the semi- 
classical equations. The spacetime metric is determined by 
(i) the incoming first and second order metric perturbations 



'■ab 



and h 



(2) .in 
ab 



at J' , which describe freely propagating 
gravitational radiation, and (ii) the two-point functions F-j^^ 
and F>^' at J of the first and second order perturbations 



,(1) 



and (aj^'^^ to the incoming state a;in(e) of the scalar field. 
These determine, via the semiclassical equations, correspond- 
ing outgoing quantities at . However, we choose to pa- 
rameterize the solutions in terms of the incoming quantities 
at J-. 



The "transverse smearing" of the ANEC integral 
referred to in the previous subsection is defined in 
Minkowski spacetime as follows. Let 7 be a null geodesic, 
let A, a;^, 
given by C = 
is 



X be coordinates such that the geodesic is 
x"^ = 0, A = 1,2, and such that the metric 



ds^ = -2d\dC + {dx^f + {dx 



2\2 



(1.9) 



Now let 



dX / d'^x'^Txx{X,C = 0,x\x'^)S{x\x'^), (1.10) 



where S is some positive function that is peaked at 
a;'^ = and that falls off at large x^. This is essentially 
an ANEC type integral, with an additional averaging in 
two spatial directions transverse to the geodesic. The 
definition (1.10) can be generalized in a natural way to 
general, curved spacetimes by using Fermi- Walker type 
coordinates [see Sec. VA below for details], provided, of 



course, that the smearing function, S, vanishes outside 
the region where such coordinates are well defined. 

Let 7(e) be a one-parameter family of null geodesies 
in (M, gabis)), and let /5(e) denote the smeared ANEC 
integral along 7(e). We expand as Is = e/i^^ +e^/i^^ + 
0{e^). Using the reduced-order semiclassical equations 
we calculate the dependence of the perturbations of the 
ANEC integral on the initial data: 



j(i) 

1(2) 



(1.11) 



We restrict attention to the case where the incoming two- 
point functions F^^^ and F^^^ satisfy the regularity condi- 
tions discussed in Appendix C. We choose the transverse 
smearing function S'(x) that enters into the definition 
(1.10) oils to be 



S'(x) oc 



1 



l-KxVAf, ' 



(1.12) 



where x = (x^,a;^), and At is greater than a certain 
critical length, AT,crit, of the order of the Planck length 
[see Sec. VI below]. We also specialize to the limit in 
which the lengthscales determined by the incoming state 
arc much larger than the Planck length. With these as- 
sumptions, our results may be summarized as follows (see 
table I below): 

At first order in perturbation theory, we find that (i) 

Is^^ vanishes for pure incoming states, (ii) /i^^ is posi- 
tive for mixed incoming states if the transverse smearing 
length At is chosen to be greater than Ax.crit! and (iii) 
/i^^ can be made negative if At is chosen to be suffi- 
ciently small, and in particular, if there were no trans- 
verse smearing at all. Thus, the transverse smearing is a 
crucial ingredient in our analysis. Note also that result 
(i) is a necessary condition for Ig to be always nonneg- 
ativc; if li^"' were nonvanishing for some incoming pure 
state;, it could be made to have either sign by choosing the 
sign of the first order perturbations appropriately. [See 
Sec. IIIB below for an explanation of why pure states 
and mixed states behave differently in this regard.] As 
was already mentioned in the previous subsection, the 
analog of result (i) fails to hold when the semiclassical 
Einstein equation is not imposed. 

Because the smeared ANEC integral vanishes at first 
order for pure states, it is necessary to go to second 
order in perturbation theory to see if the positivity of 
this integral can be violated for pure states. A com- 
plete calculation of all second order effects would have 
required us to derive a formula for the complete correc- 
tions to the stress-energy tensor of the quantmn field 
valid to second order in deviation from fiatness. This 
would have been a major undertaking in its own right, 
and we did not attempt to do this. Instead, we spe- 
cialized to a limit in which not only are the lengthscales 
determined by the incoming state much larger than the 
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Planck length, but, in addition, the incoming classical 
gravitational radiation does not dominate the first or- 
der metric perturbation. Under these conditions, the un- 
known term in the second order stress energy tensor is 
negligible. In Sec. VII, we calculate the three remaining 
terms in [F^^^ , F^^^ , /i^^^ ''"] , and we show that in this 

(2) 

limit Is is always positive. As before, the positivity 
only holds with transverse smearing. 

One unsatisfactory feature of our analysis is the follow- 
ing. The positivity of the transversely smeared ANEC in- 
tegral holds in our perturbation expansion only when the 
transverse smearing function S falls off like a;""* or more 
slowly. In particular, the smeared ANEC integral with 
transverse smearing of compact support is not always 
positive. However, in a curved spacetime, the Fermi- 
Walker type coordinates needed for the generalization of 
the definition (1.10) will be well defined only in a neigh- 
borhood of the null geodesic in question. Thus, although 
we prove the positivity of a smeared ANEC integral in 
the context of perturbation theory, we do not even have 
an obvious candidate for a "smeared ANEC conjecture" 
outside of this context. Nevertheless, we interpret our 
results as having the physical implications described in 
the previous subsection. 

Finally, it is worth mentioning why we need to con- 
sider general, mixed incoming states instead of just pure 
states. In most situations, whatever is true for pure states 
will generalize trivially to mixed states. Here, however, 
it turns out that pure states and mixed states behave 
qualitatively very differently. One reason for this (which 
may be viewed as an artifact of semiclassical theory) is 
that the linearity of the evolution law for states breaks 
down as a result of the coupling to the classical metric. 
This arises only at second order in perturbation theory. 
However, even at first order in perturbation theory where 
/i^^ [i^in''] depends hnearly on OJ^^^ it is not true that pos- 
itivity of Is^^ for all pure states would imply positivity 
for mixed states. To sec this, let V denote the space of 
states, and let Vp denote the space of pure states. Then 
the space of allowed state perturbations ool^^ is not V but 
is the tangent space T^g^dT^) to V at ojq = |0) (0|. Simi- 
larly, the linear space of allowed pure state perturbations 
is the tangent space Tvac('E'p) to T>p at wq. Now the 
key point is that, although V is the convex hull of X'p, 
Tva.c{T^) is not the convex hull of T^a.ci'Dp) but is larger 
than it. [The convex hull of T^adT^p) is just itself since 
it is a linear space.] Therefore, results for pure states 
do not generalize to mixed states. Thus, the differences 
between pure and mixed states arise in our analysis be- 
cause of our working in perturbation theory. Roughly 
speaking, we find that mixed states at first order in per- 
turbation theory behave very similarly to pure states at 
second order. 



D. Conventions 

We use the metric signature (— , -|-, -|-, -|-), and the sign 
conventions of Rcfs. [44] and [5], as well as the abstract 
index notation explained in Ref. [5]. The following is 
our convention for Fourier transforms, li F = F{x) is 
a function on Minkowski spacetime, then we define the 
Fourier transform to be 

F{k) = j (fxe-'^-'' F{x). (1.13) 

Similarly if / = /(x) is a function on three dimensional 
Euclidean space, we define its Fourier transform to be 

/>) = j d'xe-^'^-^fi^). (1.14) 

We use gravitational units in which the speed of light c 
and Newton's gravitational constant G are unity, so that 
U = Lp, where Lp is the Planck length. An index of 
notation is given at the end of the paper in table II. 



II. SEMICLASSICAL GRAVITY WITH 
BACKREACTION 

Difficulties arise in the calculation of backreaction ef- 
fects in semiclassical gravity for several reasons. First, 
as we discuss further in section IV, there are problems 
associated with spurious, runaway solutions to the equa- 
tions. Second, the dependence of the renormalized stress 
tensor on the spacetime geometry is nonlocal and compli- 
cated, making calculations outside of linear perturbation 
theory prohibitively difficult. These difficulties do not 
appear in two dimensional models, for the reasons ex- 
plained in Ref. [45], and in recent years there have been 
a variety of calculations of backreaction effects on two di- 
mensional black hole backgrounds [46]. However, in four 
dimensions, most backreaction calculations have been re- 
stricted to linear order perturbation theory off some fixed 
background (except for analyses of various conformally 
flat spacetimes [36,47]). Our backreaction analysis is ap- 
parently the first to go beyond linear order perturbation 
theory for a generic class of four dimensional spacetimes. 
We shall build upon and extend the work of Horowitz 
[48], who considered semiclassical gravity in perturba- 
tion theory about flat spacetime to linear order, without 
allowing perturbations to the incoming quantum state. 

A. The classical equations 

We now describe in more detail the model of gravity 
coupled to a scalar field that we will treat in our analy- 
sis. We consider a scalar field $ of arbitrary curvature 
coupling ^ and arbitrary mass m, so that the Lagrangian 
is 
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- m2$2 - ^m^}, (2.1) 

where k = 1/{8ttG). We will specialize later to the mass- 
less case. The classical equations of motion are 



KGabigcd] = {Tab[gcd])u 



(2.6) 



1^ Gab — Tabt 



(2.2) 



and 



(2.3) 



[□ - rn^ - ^R] $ = 0, 
where the stress tensor is 

Tab = Va$V6$ - ^9ab{V^? ' IgabmH^ 
+ C [Gab^^ - 2Va($V6$) + 2g„bV^(Wc$)] . (2.4) 

In the classical theory, this stress tensor violates the 
pointwise mill energy condition for ^ / 0, but satisfies 
the averaged null energy condition for all negative val- 
ues of ^, provided that suitable asymptotic fall-off con- 
ditions hold for More precisely, for any null geodesic 
2.0 _ (A) with affine parameter A and tangent vector 
A« = {d/d\Y, it follows from Eqs. (2.2) and (2.4) that 



jG.bX''\'dX = jd\^^^ 



+4^2 / 



(K-e$2)2' 



(2.5) 



provided that ^ < and ($$')/(« - C*^) ^ as |A| 
00, where primes denote derivatives with respect to A. 
However, for at least some positive values of ^, it seems 
likely that ANEC can be violated [49]. 

The failure of the classical stress-energy tensor (2.4) 
to satisfy pointwise energy conditions when ^ 7^ 0, and 
its possible failure to satisfy the classical ANEC condi- 
tion for positive values of ^ indicates that this class of 
field theories may not be suitable for a general study 
of the validity of ANEC in quantum field theory in 
curved spacetime. However, in this paper, we shall be 
concerned only with a perturbative analysis of solutions 
about Minkowski spacetime. In this case, it follows im- 
mediately from Eq. (2.5) that, for all values of ^, the 
classical ANEC condition holds for all solutions near 
the trivial solution of Minkowski spacetime with $ = 0. 
Thus, this class of models should provide a good testing 
ground as to whether, in the context of perturbation the- 
ory, quantum fields can attain violations of ANEC which 
are classically forbidden. 



B. The status of the semiclassical equations and the 
local curvature ambiguity 

It is natural to postulate the semiclassical equation 



to describe the backreaction of the quantized matter 
degrees of freedom on spacetime. However, as is well 
known, the exact status and domain of applicability of 
this equation is far from clear, due mainly to the fact 
that we do not have a complete, more fundamental the- 
ory of quantum gravity from which Eq. (2.6) could be 
derived. In addition there exist pathological runaway- 
type solutions to these equations, as analyzed in detail 
by Horowitz [48] for the linearized equations. Proposals 
for dealing with these unphysical solutions have been put 
forward by Simon [33-36], and we will discuss Simon's 
suggestions in detail in Sec. IV below. In this subsection, 
we critically examine the origin, status and uniqueness 
of the semiclassical equations. Specifically, we consider 
the following two basic issues: 

• How is the "semiclassical approximation" derived 
and what is its domain of applicability, i.e., what 
is the class of states uj and classical metrics Qab for 
which the semiclassical approximation is good? 

• How unique is the semiclassical Einstein equation 
(2.6) itself, i.e., how unique is the prescription to 
obtain the expected stress energy tensor (Tab)? 

With regard to the first issue, it should first be noted 
that the classical metric, gab^ appearing in the semiclas- 
sical equations presumably must correspond from the 
vantage point of a complete quantum theory of gravity 
— to the expected value of a quantum metric operator 
in some state. However, we wish to emphasize here that 
the "average value" of a metric which has nonzero am- 
plitudes to correspond to different spacetime geometries 
is an intrinsically gauge non-covariant concept; it cannot 
even be defined unless an algorithm is given which com- 
pletely fixes all gauge freedom, i.e., which rigidly fixes 
a coordinate system for each spacetime geometry. This 
phenomenon can be seen even in the simple case of a 
classical probability distribution which assigns probabil- 
ity 1/2 to the spacetime {M^g^^j^) and probability 1/2 
(2'] 

to the spacetime {M.g^^,'). We could say that the ex- 
pected metric on M is (gab) = {g^lb + 9^ab )l'^- However, 
we could equally well have represented the second space- 
time as (M, 'fp*9ab )' where tp : M ^ M is any diffeomor- 
phism. The expected metric would then be computed to 

be (g'ab) = (ffil^ + rgT)/"^- However, (g,,) and {g'J 
will not, in general, be gauge equivalent (nor is there any 
guarantee that either of them will even define Lorentz 
metrics!). Thus, unless one eliminates all gauge freedom, 
the notion of an "expected metric" will make sense only 
in a limit where the fluctuations of the spacetime geom- 
etry about its mean value are negligibly small. Such a 
limit undoubtedly would have to be contemplated in any 
case in order to justify a semiclassical approximation, but 
it is worth bearing in mind that the quantity one is trying 
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to calculate in the semiclassical theory is ill defined (with- 
out complete gauge fixing) except in this limit. Note that 
this same phenomenon would occur in the calculation of 
the expected value of any gauge dependent quantity in 
any non-abelian gauge theory. 

There are at least two different approximation schemes 
by which the semiclassical equation (2.6) can be formally 
derived, and these give rise to distinct viewpoints on the 
range of validity of this equation [50] . In the first scheme 
one formally expands a quantum metric and quantum 
scalar field about Minkowski spacetime (or, more gen- 
erally, any classical vacuum solution [53]), and derives 
the equation of motion (2.6) for the expected metric by 
keeping only the "tree diagrams" for gravitons and the 
"tree diagrams" and "one-loop" terms in the scalar field 
[54,55]. [As is well known, predictions can be obtained 
at one loop order despite the fact that the theory is non- 
renormalizable, at the expense of having to introduce two 
new, undetermined, coupling constants into the theory; 
sec Eq. (2.9) below.] This loop expansion formally cor- 
responds to an expansion in powers of ?i, so keeping only 
the "one-loop" terms corresponds to keeping only the 
lowest order correction in h to the equations of motion. 
There is no obvious mathematical or physical justifica- 
tion for dropping the one-loop "graviton" terms, since 
the quantum effects of the metric field are, a priori, just 
as important as those of a scalar field [56]. It is possi- 
ble that the neglect of the one-loop graviton terms could 
be jTistified for certain choices of the state of the scalar 
and gravitational fields — wherein the expected scalar 
stress-energy tensor dominates the corresponding effec- 
tive graviton contribution, and the fiuctuations in the 
scalar stress-energy tensor are suitably small — but we 
are not aware of any analysis demonstrating this. How- 
ever, in any case, the effects of the graviton loops would 
be expected to be qualitatively similar to the effects of 
the scalar field loops, so, at the very least, Eq. (2.6) can 
be justified as a simplified model of the exact equations 
resulting from keeping all one-loop terms. Note that 
the one-loop graviton terms could be handled within this 
approximation by treating the metric perturbation as a 
linear field propagating in a background classical space- 
time. However, the phenomenon described in the pre- 
vious paragraph would then manifest itself by the fact 
that, at second order in the perturbed metric, the differ- 
ence between the Einstein tensor of the expected metric 
and the expected value of the Einstein tensor (which ef- 
fectively acts as a "graviton stress energy tensor" ) would 
be gauge non-covariant. 

Thus, in this first scheme, the semiclassical Einstein 
equation (2.6) would be, at best, an approximate equa- 
tion valid for certain states of the scalar and gravita- 
tional fields, and, at worst, be a "model equation" whose 
properties should be qualitatively similar to the equa- 
tion resulting from a complete one-loop approximation. 
In either case, higher loop contributions would modify 
Eq. (2.6) by terms proportional to quadratic and higher 
powers of ?l. In the context of perturbation theory off of 



Minkowski spacetime, a formal expansion in powers of h 
is equivalent to an expansion in powers of l/C^ (keeping 
constants of nature such as Ti fixed), where C denotes 
a typical lengthscale associated with a solution. Thus, 
from the viewpoint of this first scheme we would expect 
there to be modifications to Eq. (2.6) which are of order 
^ Lp/C"^, and consequently the domain of validity of the 
equation is restricted to C ^ Lp, where Lp is the Planck 
length. 

In the second scheme the semiclassical Einstein equa- 
tion (2.6) is formally derived from a fully quantum treat- 
ment by imagining that there are N decoupled scalar 
fields present — all of which are in the same quantum 
state — and then taking the limit iV ^ oo, with GN = 
(constant) (see, e.g., Ref. [58]). In this "1/A^ expansion", 
the graviton loops are suppressed relative to the matter 
loops simply because there are N scalar fields but only 
one graviton field. Note that since the scalar fields are 
free (i.e., non-self- interacting), only one-loop terms in the 
scalar fields arise. Furthermore, in the iV — s- oo limit, the 
fiuctuations in the expected total stress-energy tensor of 
the scalar fields becomes negligible. Thus, one formally 
obtains Eq. (2.6) exactly in the N ^ oo limit. Correc- 
tions to this equation should be of order or higher. 
Thus, one might expect that if N is sufficiently large, 
Eq. (2.6) would be a good approximation up to the cur- 
vature scale corresponding to the effective Planck length 
in the rescaled theory, C ~ Lp,eff = VNLp. Indeed, es- 
timates of the order of magnitude of successive terms in 
the loop expansion support this viewpoint. However, it 
is far from clear that the loop expansion would provide a 
good approximation in this regime, so we do not believe 
that there are solid grounds for believing that "graviton 
effects" can be neglected for any finite value of N when 
C ^ ip.off- In other words, we would expect effects which 
are non-perturbative in H could be important at these 
scales. One piece of evidence that this is the case is that 
this seems to be the only way to escape the conclusion 
that flat spacetime is unstable [58] . As we shall see in de- 
tail in Appendix A below, pathological solutions of the 
linearized version of Eq. (2.6) exist on scales C ~ ip,eff, 
and we will not regard these solutions as being physical. 

We comment that the existence of these two different 
points of view on the status of the semiclassical equa- 
tion has given rise to some controversies in the literature. 
Some implications of the first point of view have been 
discussed at length by Simon [33-36] (see Sec. IV below 
for further discussion). On the other hand, Suen [59] 
disagreed with Simon's analysis and argued that there 
should not be any corrections to Eq. (2.6). In our view, 
the existence of higher order corrections to Eq. (2.6) de- 
pends on whether or not one justifies that equation in 
terms of a one-loop approximation or by the invocation 
of a limit. Certainly, corrections to Eq. (2.6) will 
appear in the physically realistic case of finite N. 

The viewpoint will shall adopt in this paper is the fol- 
lowing: We will appeal to the 1/N limit to give a mathe- 
matically clean justification for ignoring graviton contri- 
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butions to Eq. (2.6). Thus, in an h (or, cquivalently, a 
"long wavelength" ) expansion, we shall regard Eq. (2.6) 
as valid to all orders in Ti (or, equivalently, to all orders 
in inverse wavelength), although in most of our analy- 
ses we will not make use of its validity beyond order h"^. 
However, as already indicated above, we will not view 
Eq. (2.6) (or any other known equation) as adequate for 
describing phenomena where wavelengths comparable to 
Lp occur. This will justify our modification of Eq. (2.6) 
in Sec. IV below using the reduction of order procedure, 
so that it still describes "long wavelength" behavior ac- 
curately, but no longer predicts pathological behavior at 
the Planck scale. 

We now turn to a discussion of the second issue, i.e., 
the uniqueness of the right hand side of Eq. (2.6). If we 
assume that the prescription for obtaining the expected 
stress t(;nsor satisfies the axioms discussed in Ref. [45], 
then the expected stress tensor is unique up the addition 
of local, conserved curvature tensors. These are tensors 
which are functionals of the metric, and whose value at a 
point depend only on the geometry in an arbitrarily small 
neighborhood of that point. This ambiguity is assumed 
to be only a two parameter ambiguity [60,61], because in 
general spacetimes there are only two independent con- 
served local curvature tensors of dimensions (length)"^, 
which are explicitly given by [62] 

+ ^RRab + {RcdR'^ab - ^RacbdR"^ (2.7) 

and 

= 2VaVbR - 2UR g^h +Ir^ 9ab - 2i? Rub- (2.8) 

Thus, if we denote by (T'a6)point split the stress tensor 
given by the point splitting algorithm (briefly reviewed 

in the next subsection), then the appropriate right hand 
side for Eq. (2.6) must be of the form 

(Tab) = {Tab) point split + aAab+pBab. (2.9) 

Here a and /3 are dimensionless coefficients which can 
be regarded as free parameters. Different renormaliza- 
tion schemes predict different values of a and /3. Indeed, 
a and /? may be viewed as new "coupling constants" 
which must be introduced into the theory as a result of 
the non-renormalizability of quantum gravity (coupled 
to matter) at one-loop order. In quantum gravity, at 
higher (graviton) loop orders, additional new "coupling 
constants" would have to be introduced as coefficients of 
the conserved local curvature terms of the appropriate di- 
mension for that order; these would constitute a portion 



of the 0{h^) and higher corrections to the semiclassical 
equation referred to above. 

An important feature of (Tab) is that it has an anoma- 
lous behavior under a scaling of the spacetime metric and 
the corresponding scaling transformation of the state. 
To see this more explicitly, consider, for simplicity, the 
case of a massless field. Under the scaling transforma- 
tion gab f^^gab, we must scale the 2-point function, G, 
of the quantum field as G{x, y) ^~'^G{x, y) in order 
to preserve Eq. (2.15) below. Physically, this scaling of 
the quantum state can be interpreted as preserving the 
"particle content" of that state, so that the new state 
in the new metric fJ^gab corresponds simply to increas- 
ing the wavelength of all of the particles by the factor 
jjL. Note that this required scaling of G contrasts sharply 
with the situation in the classical theory, where the am- 
plitude, of the scalar field may be scaled in an arbi- 
trary manner independently of the scaling of the metric 
tensor [63] . In the classical theory, the stress-energy ten- 
sor also scales in a straightforward manner under any 
combined scaling transformation of gab and However, 
the situation is quite different in the semiclassical the- 
ory because one is forced to introduce a lengthscale, Aq, 
in the prescription for defining {Tab)- In particular, in 
the point-splitting algorithm (reviewed briefly in the next 
subsection), a lengthscale implicitly enters into the log- 
arithmic term in the local Hadamard subtraction term 
G^(x, y) [45]; a lengthscale similarly enters all other reg- 
ularization prescriptions as a "renormalization point" or 
"cutoff". Under a change of this lengthscale Aq — » /xAq, 
we have [64,42] 

{Tab) [gcd, /"Ao] = {Tab) [9cd, Ao] 

-F47r logM [aAab + bBab] , (2.10) 

where a, b are speciflc numerical coefficients that depend 
on the curvature coupling ^ [see Eq. (3.34) below]. Equiv- 
alently, if we keep Aq fixed but scale the metric and state 
via 

9ab l^^9ab 

G{x,y) ^ fi-^Gix,y). (2.11) 

we obtain 

{Tab) [M^PcdiAo] = IJ.~'^{Tab) [fl'cdjAo] 

-47r/i-^logM [aAab + hBab] • (2.12) 

Note that the ambiguity (2.10) in {Tab) resulting from 
the need to introduce a lengthscale is subsumed by the 
more general ambiguity given by Eq. (2.9). Indeed, one 
way of describing the above anomalous scaling behavior 
(2.10) is to say that a particular linear combination of 
the two new dimensionless "coupling constants" a and /3 
is a "running coupling constant", i.e., in effect, its value 
depends upon the scale one is considering. 

The scaling behavior given by Eq. (2.10) has two im- 
portant consequences. First, it shows that (at least part 
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of) the ambiguity occurring in Eq. (2.9) for a mass- 
less field cannot be eliminated by any criterion arising 
from the study of the quantum field theory of that field 
propagating in a fixed classical background spacetime, 
since that theory does not have a preferred lengthscale, 
whereas any prescription fixing a and /? would have the 
effect of determining a lengthscale. (Although a massive 
field does have an associated lengthscale - namely 1 /m 
- using this lengthscale to fix a and /3 would give rise to 
singular behavior in the m ^ limit.) The second conse- 
quence is that the scaling behavior (2.10) will affect the 
nature of solutions in the "long wavelength" limit. As 
we will discuss in detail in Sec. Ill E below, when we per- 
form an expansion in a "wavelength parameter" Lp/C, 
we will need to introduce terms which vary as powers of 
ln[Lp/£](Lp/£^) as well as powers of Lp/C"^. 

Although the numerical values of the dimensionless co- 
efficients a and (3 are not known - and, indeed, are not 
determinable without a more complete theory we shall 
assume below that their values are of order unity "at the 
Planck scale" . More precisely, if we define (ra6)point split 
by choosing the lengthscale Aq arising in that prescription 
to be Lp, then we shall assume that the correct formula 
for (Tab) is given by Eq. (2.9) with a and /3 of order unity. 
However, our results concerning ANEC will be valid for 
all values of a and /3. 

C. The semiclassical equations 

Consider now the semiclassical theory where the met- 
ric gab is treated classically, but where the scalar field <& 
is treated as a quantum field. The metric and quantum 
state, ui are required to satisfy the semiclassical Einstein 
equation (2.6). The coupled, evolving degrees of freedom 
in the semiclassical theory consist of (i) the metric, and 
(ii) all the obscrvablcs associated with the scalar field. 
These field observables may be taken to be the n-point 
correlation functions ($(a;i)...$(x„))aj in the given quan- 
tum state, Lu. However, as we discuss further below, the 
expected value of the stress tensor in any state is deter- 
mined via the point splitting procedure from a knowledge 
of only the two-point function, 

G{x,y) = i^xMy))^. (2.13) 

Moreover, in free field theory (which we are considering 
here) the evolution of the two-point function is decou- 
pled from that of the higher n-point functions. Thus, 
if we are only interested in the metric and not in other 
observables depending on the state of the scalar field, we 
can regard the semiclassical equations as a set of coupled 
equations for the metric gab and the distributional bi- 
solution G{x,y) to the scalar wave equation (2.3). This 
is a key feature which simplifies our analysis. Prom this 
point of view, states which differ only in their n point 
functions for n 7^ 2 are effectively identical as far as semi- 
classical gravity is concerned. 



The appropriate set of bi-distributional solutions for 
a given, fixed, globally hyperbolic spacetime (M, Qab) 
can be characterized as follows [65,45]: Let S{M) de- 
note the space of smooth solutions of the Klein-Gordon 
equation (2.3) with initial data of compact support on 
Cauchy surfaces, and let (M) be the space of smooth 
test functions of compact support on spacetime. Define 
the usual Klein-Gordon-like symplectic product on pairs 
F,G& S{M) 

n{F,G) = - J^FVaGdi:'', (2.14) 

where E is any Cauchy surface. Let G{f, g) denote 
the two-point bidistribution evaluated on ("integrated 
against") test functions f,g & C^{M). Then G{f,g) 
must be of the form 

G{f,g) = ^G«(/,5) + \mf,Eg), (2.15) 

where Ef denotes the advanced minus retarded solution 
with source /. Furthermore, the symmetric part G(^)/2 
of G must satisfy the positivity conditions G^^\f, /) > 
and 

G(i)(/,/)G(i)(3,5) > n{Ef,Egf. (2.16) 

Moreover, in order that the stress tensor be well de- 
fined, the two-point function G{x, y) must be of so-called 
Hadamard form [65]. 

To specify more explicitly the coupled evolution equa- 
tions for the metric gab{x) and the two-point function 
G{x,y), we need to discuss the point splitting prescrip- 
tion for calculating the stress tensor. Let G^{x, y) denote 
the locally constructed Hadamard bidistribution given 
by the algorithm described in Refs. [65,45,66], as speci- 
fied/modified in the following way. Use the differential 
operator 

V = U-rr? -S,R. (2.17) 

Choose the lengthscale, Aq, implicitly appearing in the 
logarithmic term in G^{x,y) to be the Planck length 
Lp, choose w7o = in the notation of Rcf. [45], and trun- 
cate the series expansion after three terms. Then the 
regulated two-point function 

f{x,y) = G{x,y)-G"{x,y) (2.18) 

will be well defined in a neighborhood of the "diagonal" 
X = y of M X M and will be at least for Hadamard 
states [65]. The expected value of the stress tensor will 

be given by 

(T'a6(a;))point split = lim Vabf{x, y) + Q{x)gab{x), (2.19) 

3/— >x 

where Vab is a particular second order differential op- 
erator and (5 is a particular local curvature invariant 
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[45]. Note that if we modify the prescription for cal- 
culating by truncating the series after say four terms 
instead of three terms, then the regularized two-point 
fimction f{x,y) will be altered, but the value (2.19) of 
(rah(.T))point split will not be changed. As explained in 
the previous subsection, we shall assume that the correct 
value of {Tab{x)) is given by Eq. (2.9), with a and /3 of 
order unity. 

In summary, the independent variables in the semi- 
classical evolution equations consist of a smooth metric 
Qabix) and a bidistribution G{x, y). G{x, y) is required to 
satisfy the wave equation (2.3) in each variable, as well as 
Eqs. (2.15), (2.16) and the Hadamard condition. Finally, 
Qab and G are required to satisfy the scmiclassical Ein- 
stein equation (2.6), with {Tab) given by Eqs. (2.9) and 
(2.19). The gauge freedom in this formulation of semi- 
classical gravity simply consists of the diffeomorpliisms 
if : M ^ M, under which gab and G{x,y) get trans- 
formed by the natural action of 

Although the above formulation of scmiclassical grav- 
ity is fully satisfactory mathematically, the unknown 
variable G{x,y) has a distributional character, and it 
would be more convenient to specify the independent de- 
grees of freedom in terms of a smooth function. This can 
be done as follows in the physically relevant case of a 
spacetime (M, gab) which becomes flat in the asymptotic 
past. First, we assume that the state is sufficiently reg- 
ular and that the approach to flatness of the spacetime 
occurs at a sufficiently rapid rate that G{x,y) asymp- 
totically approaches the two-point function Gin{x,y) of 
a state, LUin, in Minkowski spacetime {M^rjab)- For suffi- 
ciently regular states of a massless field, asymptotic flat- 
ness of (M, gab) at null infinity should suffice for this to 
hold, with G{a{x,y) being the two-point function of the 
state in Minkowski spacetime with the same initial data 
on past null infinity J~ as G(x, y) [under a suitable iden- 
tification of J~ for {M,gab) with J~ for {M,rjab)]- For 
a massive field, it is less clear precisely what asymptotic 
conditions on {M,gab) would sufiice, but the necessary 
conditions presumably would be qualitatively similar to 
those for a massless field. Note that the positivity con- 
dition will hold for G{x, y) if and only if it holds for 
Gin(x,i/). This is because, in general, the positivity con- 
dition (2.16) can be expressed as a condition on initial 
data on E X S, where S is any Cauchy surface [65], and 
will be preserved under evolution. The two-point func- 
tions G and Gin have the same initial data at , and 
are both evolved forward using the homogeneous wave 
equation with respect to the appropriate metric [gab for 
G, -qab for Gin). 

Now, let Ginfi{x,y) denote the two-point function of 
the ordinary vacuum state, Win.o, in Minkowski space- 
time {M,riab), and let Go{x,y) be two-point function 
of the corresponding state, wq, in {M,gab) which ap- 
proaches Gin,o{x,y) in the asymptotic past; in other 
words, let Go{x,y) be the two-point function of the "in 
vacuum state" in (M, (/<,;,). For any Hadamard state, co, 
in {M,gab) with two-point function G{x,y), we define 



F{x,y) = G{x,y)-Go{x,y). (2.20) 

Then, clearly, F contains the same information as the 
bidistribution G. However, the Hadamard condition on 
G now corresponds simply to the statement that F{x, y) 
is a smooth fimction on M x M . Furthermore, Eq. (2.15) 
is equivalent to the statement that F is symmetric and 
real-valued. In addition, F{x, y) satisfies the wave equa- 
tion (2.3) in each variable 

V^F{x, y) = VyF{x, y) = 0, (2.21) 

by virtue of the fact that the bi-distributions G{x, y) and 
Ga{x,y) each satisfy this equation. The only other re- 
striction on F is the positivity condition arising from the 
corresponding condition (2.16) on G. 

Since the expected stress tensor {Tab)ij has a linear 
dependence on the two-point distribution, G, associated 
with oj, we have, for any Hadamard state, ui, 

{Tab)u: = {Tab)u,o + {Tab) F ■ (2.22) 

Here as above ojq denotes the "in" vacuum state, and 

{Tab)F=\^vtiVabF{x,y), (2.23) 

with Vab being the same difii'erential operator as appeared 
in Eq. (2.19). The vacuum polarization term (^ab)uo is 
functional of the spacetime metric gab alone. For massless 
fields, it has been evaluated by Horowitz [48] to first order 
in the perturbed metric about fiat spacetime, and we will 
make use of Horowitz' formula in our analysis below. 

Thus, in our reformulation of scmiclassical theory, 
the independent variables consist of a smooth metric 
gab{x) and a smooth, symmetric, real- valued function 
F{x, y) which satisfies the positivity condition arising 
from Eq. (2.16). The dynamical evolution equations are 
simply that F{x,y) satisfy the wave equation (2.3) in 
each variable with respect to the metric gab, and that gab 
and F satisfy the semiclassical Einstein equation (2.6) 
with stress-energy tensor given by Eqs. (2.22) and (2.23). 
We remark that a purely classical theory, where gravity 
is coupled to a statistical ensemble of scalar field con- 
figurations, would differ from this version of semiclas- 
sical gravity only in the following two respects. First, 
in the classical theory, the term {Tab)ujo would be ab- 
sent from Eq. (2.22). One may view this term, which 
is proportional to h, as describing the vacuum polariza- 
tion effects occurring in the quantum theory. Second, 
the quantum mechanical positivity condition on F aris- 
ing from Eq. (2.16) (which can be viewed as a restriction 
on the space of allow(xi initial data for _F on E x E for 
any Cauchy surface E) is less restrictive than the corre- 
sponding classical positivity condition. 

Under our above assumptions about the asymptotic 
behavior of states, it follows from Eq. (2.20) that in the 
asymptotic past F{x,y) approaches the smooth function 
Fin{x,y) on Minkowski spacetime defined by 
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Fin{x, y) = Gin{x, y) - Gi„,o(a;, y). (2.24) 

Note that Fin is just the usual "regularized" two-point 
function of the state Win in Minkowski spacetime. It 
satisfies a positivity condition which ensures that Gin is 
the two-point distribution of a state in Minkowski space- 
time. We may view Fi^{x^y) and a corresponding quan- 
tity describing the incoming classical gravitational radi- 
ation [sec Eq. (1.8) above] as the freely specifiable "ini- 
tial data" for semiclassical gravity. In the next section, 
we will develop a systematic perturbation expansion for 
semiclassical solutions in terms of this "initial data" . 

Wc conclude this section by reminding the reader of 
our notation for the four different states under considera- 
tion here: (i) the state of interest, u), on the curved space- 
time (M, Qab) (ii) the "incoming state" Win on Minkowski 
spacetime {M,r]ab), such that the n-point hmctions of to 
approach those of Win in the asymptotic past (iii) the 
vacuum state, Win,o in Minkowski spacetime, and (iv) 
the corresponding "in-vacuum" state luq in the curved 
spacetime {M,gab), such that the n-point functions of wo 
approach those of a;in,o in the asymptotic past. The two- 
point functions of these states are, respectively, G, Gin, 
Gin,o, and Gq. 

III. PERTURBATION THEORY ABOUT FLAT 
SPACETIME 

A. Derivation of the equations 

In this section we derive the explicit form of the pertur- 
bation expansion of the semiclassical equations off of the 
Minkowski spacetime/ vacuum solution, valid to second 
order in deviation from flatness. Let (M, ^^^(e), w(£)) be 
a smooth one parameter family of solutions of the semi- 
classical equations discussed in the last section, such that 
{M,gab{0) = r]ab) is Minkowski spacetime and w(0) is the 
Minkowski vacuum state, Wi„,o. We expand all relevant 
quantities about £ = as follows: 

gab{e)=Vab + eh';^^+6'h'^^+0{6^) 

Win = Win,0 + + 0{£^) 

F = £FW+£2f(2)+0(£3) 

Fin=£FiW+£^Fi(^)+0(£^) 

D = r'(o) +C)(e2) 
V,b = v'^2+^T^ab+0{e% (3.1) 

where Win, F, Fjn, 2?, and Dab were defined in the previ- 
ous subsection. For the remainder of this paper the op- 
erators Vabi "D^ab ^^"i ^ifc^ ^i^l implicitly understood 
to include the operation of taking the coincidence limit 
y ^ x; these operators thus act on functions on M x M 
and take values in the space of tensors on M. Note that 
it would not make sense to introduce an expansion in e 



of the family of states u}{e), as for different values of e, 
these states act on the different algebras of observables 
corresponding to the different spacetimes. For this rea- 
son, we view the state u){e) as a functional of the "in" 
state LOin{s) on Minkowski spacetime and of the spacetime 
metric gabi^:), i.e., to = u)[u)in, gcd], and we then expand 
Win(e). Note that for a fixed metric gab, the state w is a 
linear function of Win- 
We write 

(T'abbcd], Win) = {Tab[9cd]) ojlu>in,gcd]- i^-"^) 

This defines the tensor Tab on the left hand side as a lin- 
ear map on the space of states on Minkowski spacetime 
which takes values in the space of conserved tensors on 
the spacetime {M,gab)- Using the expansion (3.1) for 
gabie), wc express this linear map on Minkowski space- 
time states in the form as 

F„,[.9c.(e)]=ri°)-feT«[/.«] 

{T!,l\h('^]+T^l^[h('\h('^]}+0{e% (3.3) 

where Tj°' is the usual stress tensor operator in 
Minkowski spacetime. Equation (3.3) defines the ten- 
sors T^l^ and T^^^ which act on metric perturbations and 
pairs of metric perturbations respectively. 

The perturbativc semiclassical Einstein equations are 
obtained by inserting the expansions (3.1) and (3.3) into 
Eq. (2.6). We obtain at first order, 

«G«[/.(i)] = (Fi°\ w«) + (T«[/.(i)],Win,o) (3.4) 

and at second order, we get 

.G«[/.(^)] + .Gi^,)[/.W,/.«] = (Tr,wL^)) 

+(T(,^)[/i(i),/i(i)],Win,o). (3.5) 
Here the tensors G^^^ and G^^^ are defined by the identity 

Gablv + ah] = aGi\^ [h] + a^G^^^ [h, h] + 0{a''), (3.6) 

which holds for any tensor hcd, i-e., they are the linear 
and quadratic parts of the Einstein tensor. Explicit ex- 
pressions for G^-^^ are given in for example Ref. [5], and 
for G(2) in Ref. [67]. 

Before discussing the explicit form of the terms appear- 
ing on the right sides of Eqs. (3.4) and (3.5), it may be 
useful to write the equations in the notation appropriate 
to the case where the one parameter family of incom- 
ing states Win(e) corresponds to a one parameter family 
of density matrices /5(e) on the usual Fock space Ti. of 
states on Minkowski spacetime. If this family of density 
matrices is expanded as 

p{e) = |0)(01 + + + 0{s% (3.7) 
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then Eqs. (3.4) and (3.5) can be rewritten as 
/.G«[/zW] = (0|T«[/.W]|0)+tr [p^T^ 



(3.8) 



and 



5(i)Ti,^)[,(i) 



= tr 



^ -'■ab 



tr 



+mii\h^'^m + {o\T^i\h^'\h^'^m. (3.9) 

In Eqs. (3.4) and (3.5), the terms involving Win,o are 
the "vacuum polarization" terms, corresponding to the 
first term on the right hand side of Eq. (2.22). The term 

(-^afc^[^'"^'']''*^in.o) has been computed by Horowitz [48] in 
the massless case, and we will review Horowitz' results in 

subsection HID below. The term {T^l\h''^\ h''^'^],LJin,o) 
has not been computed, and, to avoid having to do so, 
wc will eventually pass to an approximation in which this 
term can be neglected (sec subsection HI E below) . 

We now derive explicit expressions for the non- 
vacuum-polarization terms that appear in Eqs. (3.4) and 
(3.5). Prom Eqs. (2.23) and (3.1) we obtain 

(3.10) 

Furthermore, from Eqs. (2.21) and (3.1), we obtain 



I?WFW(x,y)=0 



(3.11) 



together with similar equations involving ^/-derivatives. 
Since F Fin in the asymptotic past, the solutions to 
these equations are 



in 



(3.12) 



Here the quantity £[sx,Sy] is defined for any sources 
Sx{x,y) and Sy{x,y) satisfying V^^h^ = V^^hy, by 

£[sx,Sy]{x',y') = / d'^xGi°l{x',x)sx{x,y') 
Jm 

+ [ d'yGi%',y)sy{x',y) 
Jm 

- [ d'x j d^yGill{x',x) 
JM Jm 

xGilliy',y)s{x,y), (3.13) 

where § = Vy^^s^ = V^^Sy, and G^^] is the retarded 
Greens function for the differential operator 'D^'^\ Com- 
bining all of the above results together with Eqs. (2.22), 
(3.4) and (3.5) yields the fairly obvious relations 



X-'-ab ' 
(rp(0) 



ab '^in 



(2)\ _ 



T,iO) p(l) 
^ab ^in 

(0) p(2) 
ab ^in ' 



(3.14) 



together with 



-Vi'^FB\ 



(1) 



(3.15) 



Finally, we note that the perturbative semiclassical 
equations have the following structure: we can spec- 
ify arbitrarily the incoming state perturbations W;^^ and 

(2,) 

wU (or, more precisely, just their two-point functions), 
as well as the incoming metric perturbations h^^^ '™ and 



,(2) an 
'■ab 



[see Eq. (1.8) above]. We then may solve Eqs. (3.4) 

and (3.5) to obtain the metric perturbations hl^^j and /i^^^ 
[68]. 



B. Pure states versus mixed states 

As discussed in the introduction, our results concerning 

ANEC depend crucially on whether the incoming state 
Win is pure or mixed to first order in e. Consider first 
the case where all the states, Win(£), correspond to den- 
sity matrices, p{e), in the usual Fock space. Then the 
perturbed state will be pure to first order if and only if 

p{e) = \^{e)){^{e)\ + 0{e'), (3.16) 

where 

[*(£)) = \0)+e\i,f^+O{s'). (3.17) 
Thus the state is pure to first order if and only if 

p(i) = |0)(V'|« + W|V)(0|, (3.18) 

for some jV*)^^^ S H with {0\tp)^'^^ = 0. By contrast, the 
most general first order density matrix perturbation is of 
the form 



^(1) 



|0)(V,|« + «|V)(0|+P-(trP)|0)(0|, 



(3.19) 



where P is a positive, Hermitian trace class operator on 
7i such that P|0) = 0. Thus, the perturbed state is pure 
to first order if and only if P = 0. 

By inspection, when Eq. (3.18) holds, it can be seen 

that the two-point function F-j^^ will have the property 

that its mixed- frequency part [i.e., the part that is posi- 
tive frequency with respect to one variable and negative 
frequency with respect to the other; see Eq. (6.1) below] 
vanishes. Moreover, the converse is also true, since if the 
mixed-frequency part vanishes we have 



tr $_(u)$+(u)l = 



(3.20) 
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for any test function u, where we have used the decom- 
position 

$ = + (3.21) 

of the field operator into its positive and negative fre- 
quency parts. [In a conventional mode expansion, <E>_|_ 
would consist of the annihilation operators and the 
creation operators.] Now since the operator P is trace 
class, it is compact. Thus there will exist an orthonor- 
mal basis of the space of states orthogonal to |0) 
such that 

oo 

P = J2P^\^j){i^j\ (3.22) 

3=0 

for some pj > 0, j = 0,1,2..., where the convergence 
is in the operator norm topology. Equations (3.19) and 
(3.20) now imply that 

Y,Pj\\K{u)\^j)\f = 0. (3.23) 
i 

However, if $+(w) |V'j) = for all u, then \tpj) = |0), 
which contradicts the fact that is orthogonal to |0). 
Therefore pj = for all j, and the incoming state is pure 
to first order. 

In addition, it follows immediately from Eq. (3.19) that 

if p(^) is a possible first order state perturbation, then 
—p^^^ will be an allowable first order state perturbation 
if and only if the state is pure to first order. This has 
the implication, which we discussed in the introduction, 
that if the ANEC integral is to be non-negative gener- 
ally, it must vanish at first order for pure states, but not 
necessarily for mixed states. 

For general, algebraic states there is a more abstract 
notion of purity, which defines a state lo to be pure if it 
cannot be written in the form cu>i + {l—c)u)2 where ui and 
C02 are distinct states and < c < 1 (see, e.g., Ref. [45]). 
Thus, the pure states arc extreme boundary points of the 
convex linear space of all states. It should follow that for 
general algebraic states, if lu^^^ is a possible first order 

state perturbation off of a pure state, then —ui^^^ will be 
an allowable first order state perturbation if and only if 
the perturbed state is pure to first order. In addition, the 
positivity condition, Eq. (2.16), applied to both 0^111,0 + 

ewjjl^ and Win,o — ^'^in'' should then imply the vanishing 
of the mixed-frequency part of the two-point function 
for general first-order pure states. Thus, the conclusions 
of the previous paragraph should continue to hold for 
general, algebraic states, although we have not attempted 
to give a rigorous proof of these results. 

Our analysis of ANEC given in Sec. VI C below will 
divide into two cases, depending upon whether or not the 
mixed frequency part of the perturbed two point function 
vanishes. In the remainder of this paper, we shall use 
the terminology "pure perturbed two point function" to 



mean a perturbed two point function for which the mixed 
frequency part vanishes, and "mixed perturbed two point 
function" to mean one for which the mixed frequency 
part does not vanish. For states in the usual Fock space 
- and, presumably, also for general algebraic states - a 
perturbed state will be pure to first order if and only if 
its perturbed two-point function is pure. 

C. Gauge freedom and transformations 

In this subsection we analyze the gauge freedom in the 
perturbation equations. The gauge freedom in the one 
parameter family of exact solutions (M, ^^^(e), a;(e)) con- 
sists simply of one parameter families of diffeomorphisms 
(fe'.M^ M, where (po is the identity map [69]. Here, 
these diffeomorphisms act simultaneously on Qat and the 
two-point function G of w. 

It can be shown [70] that for an arbitrary one- 
parameter family of diffeomorphisms with ipo being 
the identity map, there exist unique vector fields ^"(^^ 
and on M, such that to order 

-I?^(2)(e)oI?^(i,(£V2). (3.24) 

Here T>r{X) : M ^ M denotes the one parameter group 
of diffeomorphisms generated by t° The vector fields 
^(1) and C"^^^ are given by the following formulae in 
terms of their actions on test functions / G C^{M): 

r(')V„/ = ^(/o^,)|,=o (3.25) 

and 

= ^(/ ° <^-)k=o - -C^d, £^(1, /, (3.26) 

where C denotes the Lie derivative. Thus, the gauge 
freedom in the second order perturbation equations can 
be parameterized by pairs of vector fields on M . 

Now let T(e) be any one parameter family of tensor 
fields on M (we suppress tensor indices), which has the 
expansion 

T{e) = + eT(i) + £^^(2) ^ q(^3)_ (3 27) 

Then from Eq. (3.24) we can calculate the transformation 
properties of the expansion coefficients T^°\ T^^) etc. We 
find 

<^*r(£) = r(°) + eT^i) + e^T^^) + 0(e3), (3.28) 
where 

=t(i)+£j(i)T(°), (3.29) 

and 

f (2) = T(2) + l£^(,)TW + 

+£j(i)T(i). (3.30) 
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Let us denote gauge transformed quantities by overbars. 
Then we find 



+Cl,,F^'\x,y). 
We use these formulae in Sec. VII below. 



(3.31) 



D. The linearized stress tensor and the expUcit form 
of the first order perturbation equation 

In this subsection we analyze the vacuum polarization 
term appearing in the first order perturbation equation 

(3.4) in the massless case, m = 0. This has been ealcu- 
lated by Horowitz [48] using an axiomatic approach, who 
obtained 



+aA«(x)+/3BW(x). 



(3.32) 



Here A^^^ and b[1^ are the Unearized versions of the local 
curvature tensors (2.7) and (2.8): 

B^ab = -'^Vab^R^^^ + 2VaVbR^^^ , (3.33) 

where i?^^ is the linearized Ricci scalar and the deriva- 
tive operators V° and □ are the zeroth order derivative 
operators associated with the flat metric 770b. The coeffi- 
cients a and b in Eq. (3.32) are constants which depend 
on the curvature coupling ^. These coefficients were given 
by Horowitz for the cases ^ = and ^ = 1/6; the general 
formulae can be derived from point splitting and are 

1 



47r(9607r2) 

(1 - 6g)^ 
47r(5767r2) ' 



(3.34) 



Note that the coefficient a is positive for all values of 
^, and moreover the corresponding coefficient for other 
fields such as Maxwell and neutrino fields is also positive 
[48]. This fact will be be relevant in our analysis below. 
The quantities a, (3 in Eq. (3.32) are free parameters, 



c.f., the discussion after Eq. (2.9). The parameter A is a 
free parameter with dimensions of length corresponding 
to the lengthscale Aq discussed in subsection II B. The 
quantity Hx is a distribution with support on the past 
light cone, rather like the retarded Greens function solu- 
tion to the massless wave equation. An explicit formula 
for Hx is 

Hx{x) = lim [-5'{G-a)@+{x) + 2-K\n[-a/X^]5'^{x)], 

(3.35) 

where a = XaX°'/2, ©+(x) takes the value 1 inside the 
past light cone and vanishes elsewhere, and it is under- 
stood that the limit a ^ 0~ is taken after integrating 
against a test function. See Horowitz [48] or Jordan [54] 
for more details. 

The expression (3.32) is essentially a special case of 
the general formula (2.9), but linearized and specialized 
to the incoming vacuum state. As such it contains a lin- 
earized version of the two parameter local curvature am- 
biguity described by the parameters a and (3. Also the 
logarithmic scaling described by Eq. (2.10) has a counter- 
part in Eq. (3.32): the distribution (3.35) has the prop- 
erty that 



Hy{x) - Hx{x) = 47rln(A/A')^^(a;). 



(3.36) 



Therefore the free parameters a, /3 and A are not inde- 
pendent. In Sec. II B we chose to make A = Lp, thus 
fixing the values of a and /3. In the linearized analysis 
here and below, we follow Horowitz [48] and choose that 
value of A which makes a = 0. Thus, the two indepen- 
dent free parameters in the linearized stress tensor are /? 
and A. Note that our assumption discussed in Sec. II B 
that a and (3 are "of order unity at the Planck scale" 
translates into the assumption that 



A ~ Lp 
/3~1. 



(3.37) 



Formula (3.32) completes the explicit specification of 
all the terms appearing in the linearized semiclassical 

Einstein equation (3.4). Thus, using Eq. (3.14) and set- 
ting a = 0, the complete, explicit form of this equation 
is 



nG^m = rW^F^^+hpB^ix) 



(0) 



+n I d^y{ Hx{x - y) 



aA^^^{y) + hB)^>{y) 



(3.38) 



For later convenience, we have explicitly inserted the fac- 
tors of h appearing in this equation [71]. It can be seen 
that the right hand side of Eq. (3.38) contains terms in- 
volving fourth derivatives of the perturbed metric, so the 
linearized semiclassical Einstein equation has the nature 
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of a fourth order intcgro-difFcrcntial equation rather than 
a second order differential equation. As previously re- 
marked, it can be solved by specifying the source term 



Sab = in 



(0) 
b ' 



(0) 



(3.39) 



and solving for the metric perturbation h^'^jj . Its exact 
solutions have been discussed in detail by Horowitz [48], 
in the special case of the homogeneous version of the 
equation, without the source term (3.39). In Appendix 
A wc obtain all solutions to Eq. (3.38) whose spatial 
Fourier transforms exist, thereby generalizing the analy- 
sis of Horowitz to allow for a non- vacuum incoming quan- 
tum state. 

As found in Appendix A and in Ref . [48] , the linearized 
semiclassical Einstein equation has, in effect, more de- 
grees of freedom than the corresponding classical equa- 
tion, and new "run-away" solutions exist. Thus, presum- 
ably, not all of the exact solutions should be regarded as 
physical. Section IV below will be devoted to the issue of 
how to extract physical information from Eq. (3.38) and 
equations of a similar character. 



E. Explicit form of the second order perturbation 
equation in the "long-wavelength" limit 

As we have previously indicated, we will find that the 
ANEC integral vanishes identically for pure states to first 
order in deviation from flatness. This is a key necessary 
condition for the validity of ANEC, but it does not, by 
itself, provide strong evidence that any version of ANEC 
actually holds for pure states. In order to investigate this 
issue further, it is necessary to go (at least) to second or- 
der perturbation theory. However, we are unable to do 
this for general perturbations because we do not have 
an expHcit expression for the term (T^^-* [/i(^\ ft,^^^], Win.o) 
appearing in Eq. (3.5). Nevertheless, as we now shall 
describe, there are limiting circumstances under which 
this unknown term will be negligible compared with the 
other terms appearing in Eq. (3.5). These limiting cir- 
cumstances correspond to the case of "long wavelengths" 
(compared with the Planck scale) of the field and per- 
turbed metric, together with the condition that the first 
order perturbed metric not be dominated by incoming 
gravitational radiation. These conditions should encom- 
pass a wide range of physically interesting and potentially 
achievable situations. In Section VII below, we will per- 
form an analysis of the validity of ANEC for pure states 
to second order in perturbation theory in this limit. 

We now give a precise description of the long wave- 
length limit in terms of one-parameter families of solu- 
tions to the semiclassical equations, in which the charac- 
teristic lengthscale £ of a solution satisfies £ ^ oo, while 
the length scales that determine the semiclassical theory, 
Lp and A [the lengthscale appearing in Eq. (3.38)], are 
kept fixed. However, we remark that this limit is equiv- 
alent to a limit in which ?i = — > with jC and X/Lp 



fixed, i.e., the long- wavelength limit is equivalent to a 
h ^ limit, as long as the lengthscale A is taken to scale 
proportionally to Lp. 

As previously discussed in subsection II B, under the 
scaling gab oi^gab of the spacetime metric, the natu- 
ral scaling of the two-point function (corresponding to 
keeping the particle content fixed but going to longer 
wavelengths by the factor a) is 



G(x,y) 



(3.40) 



Thus, in the context of perturbation theory off of flat 
spacetime, the transformation corresponding to keeping 
the incoming particle content the same but increasing 
the wavelength of the particles by the factor a is given 
by ?7a6 a'^Vab and -Fm^"-' a~^Fin''"' for all n. Equiv- 
alently, by applying the diffeomorphism ax^ to 

this transformation, we see that the "long wavelength 
limit" should correspond to the large a limit of solu- 
tions to the perturbative semiclassical equations on a 
fixed background Minkowski spacetime, {M,r]ab), with 
a one-parameter family of states whose nth order (in s) 
perturbed incoming initial data varies as 



Fj^\x,y;a) = a-^F-J^\x/a,y/a), 



(3.41) 



where Fin^"''(a;, y) is the perturbed incoming initial data 
of some fixed (i.e., a independent) state 

Win(e) = Win,o + £ wL'^ + + 0{e^). (3.42) 

Without loss of generality, we may choose to have a 
characteristic lengthscale, Cq, of order the Planck length, 
£o ^ Lp. We may then interpret the parameter a as 
measuring the characteristic lengthscale L of the initial 
state in units of the Planck length. 
The quantity 



Sab = K 



b ' 



(3.43) 



acts as a source term for the first order metric pertur- 
bation, h^^jj {x; a), in the linearized semiclassical Einstein 
equation (3.4). For the family of incoming states given 

by Eq. (3.41), .Sab scales as 



Sab{x; a) = a 



,(1) 



){x/a). 



(3.44) 



In order to examine the behavior of Eq. (3.4) under scal- 
ing, it is convenient to make the change of variables 
U^\x;a) = h^'^^ax;a). If we substitute Eq. (3.44) into 
Eq. (3.4) and use Eqs. (3.32) and (3.36), we find that the 
explicit a dependence of the resulting equation is given 

by 



(i)\ 



(T«[/.(i)],c.i„,o) 



In 



a 



(3.45) 



Here 
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(3.46) 



is the Unearized anomalous scaling contribution to the 
vacuum polarization discussed in subsections II B and 
HID above; see Visser [42] for extensive further discus- 
sion. Note that the a dependence in Eq. (3.45) exactly 
mirrors the h dependence in Eq. (3.38), provided that we 
assume that X oc Lp [see Eq. (3.36) above]. This is as ex- 
pected from our remark above concerning the equivalence 
of the ?i. ^ and long-wavelength limits, since 



c?' 



(3.47) 



The solution to the first order perturbation equation 
with source (3.44) can be written as [c.f., Eq. (1.8) above] 



(1)/ 



(3.48) 



where /i^^^ '™ is the homogeneous solution to Eq. (3.45) where 



with the same initial data as h^-^^ (x, a) at J~ and Ah^^^ 
is the retarded solution of that equation with source 
(3.44). [Note, however, that on account of the higher 
derivative nature of Eq. (3.45), there will be more initial 

data to specify for /i^^'(a;,a) than occurs for the clas- 
sical linearized Einstein equation. This situation will 
be rectified when we replace Eq. (3.45) with the re- 
duced order equation (4.20) in subsection IV D below.] 
To specify a one-parameter family of solutions, we are 
free to chose any scaling with a we wish for the ini- 
tial data for /i^^^ '™ [72] . However, a natural choice of 
scaling is suggested by the following considerations. If 
we were solving the classical linearized Einstein equation 
with source (3.44), the retarded solution would scale as 
A/i'^'n6(x; a) = a~'^ Ah'^i^{x / a) . The retarded solution 

^^ab ^^"^ linearized scmiclassical Einstein equation 
will not scale this simply on account of the anomalous 
scaling appearing in Eq. (3.45), but the dominant scal- 
ing behavior will still be of this form [sec Eqs. (3.50) and 
(3.52) below]. Thus, if we scale the initial data for the 
incoming gravitational radiation so that near J'~ we have 



(x; a) 



"■ab 



\x/a), 



(3.49) 



then from Eq. (3.48) the first order semiclassical metric 
perturbation will not be dominated by incoming gravi- 
tational radiation in the long wavelength limit a ^ oo. 
This is the situation we wish to consider with regard to 
the second order perturbation equations. 

Before considering the second order equations, it is use- 
ful to make the change of variables 



Prom Eq. (3.45), we obtain 

In 



(3.50) 



+i{Til\x^%^i^,o). (3.51) 



With the above ansatz for the form of the incoming grav- 
itational radiation, the a dependence of x'^-* will be given 

by 

= ^(1.0) + i!^^(M) + 1_ ii,2) ^ o[{lnar/a% 
a'' 

(3.52) 

where the x^^'-'-' are independent of a. 

We now turn to the second order perturbation equation 
(3.5). Note that once h^^ij has been obtained, this equa- 

(2) 

tion as an equation for h^f^ has exactly the same structure 
as the first order equation, except that the source term 
now includes additional pieces constructed from the first 
order quantities: 



^G^ab\h('^] = + {T^l\h(%u,,^,o), (3.53) 



(2) _ /^(O) , ,(2) 
^ab 



+ (Ti,^)[/.(i),ft(i)],a;i„,o) - ^G^^^[h('\h(% (3.54) 

Consider now the behavior of the source term (3.54) un- 
der the scaUng given by Eqs. (3.41), (3.42) and (3.49). If 
we make the change of variables 



4 f2). 



(3.55) 



(2) 

(which is chosen to make sj^j, independent of a in the 
large a limit), then we find that 



+^{T^\x^'\x^%^in,o). (3.56) 

(2) 

Here the quantity Z^^^' is just the second order part of 
the anomalous scaling tensor 



Zab = 47r [aAab + bBab] ; 



(3.57) 



see Eqs. (2.10), (3.33) and (3.46). In our calculation 
of the second order perturbation to the ANEC inte- 
gral in Sec. VII below, we shall work to 0(1/ a^) be- 
yond leading order, and consequently we can drop the 
last two, vacuum polarization terms from the right hand 
side of Eq. (3.56). In particular, the uncalculated term 

(^a?[x^^')X^^^])'*^in,o) may be dropped. Moreover, in this 

equation we can replace x*^^'' by its leading order approx- 
imation x^^''^'' in the large a limit, given by 



(0) 



(3.58) 



(see Eqs. (3.51) and (3.52) above). Making the change of 
variables h^'^\x;a) = h^'^\ax; a), the resulting equation 
is 
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->S^[x^^'^X^^'°^]+0(^). (3.59) 

IV. DEALING WITH THE PATHOLOGICAL 
SOLUTIONS: SIMON'S PRESCRIPTION(S). 

As seen in the previous section, a key feature of the hn- 
earized semiclassical Einstein equation is that it fails to 
remain a second order partial differential equation; rather 
it contains fourth order time derivatives of the metric. 
Consequently, the semiclassical equations have more "de- 
grees of freedom" than the classical equations, i.e., ad- 
ditional free functions — namely, the second and third 
time derivatives of the metric can be specified as initial 
data. Directly related to the presence of these additional 
degrees of freedom is the existence — for all values of 
the free parameters ^, j3 and A — of new, "pathological" 
solutions of the semiclassical equations, which grow ex- 
ponentially in time on a timescale of order the Planck 
time (see Appendix A and Ref. [48]). In addition, for all 
^, /3, and A except for ^ = 1/6 and A > Acrit) there exist 
solutions with spatial wavelength Lp which oscillate 
in time at a frequency of order the Planck frequency. 

What attitude should one take to these solutions? If 
one justifies the semiclassical equations via the 1/N ap- 
proximation (see Sec. II B above) and formally takes the 
limit N ^ oo, one could view the semiclassical equa- 
tions as holding exactly. In that case, one should take 
all of its solutions seriously and conclude that flat space- 
time is unstable in semiclassical gravity. However, as al- 
ready indicated in Sec. II B above, we do not adopt this 
view here. Rather, although we view the semiclassical 
equations as accurate - for sufficiently large N — to ar- 
bitrarily high order in a long wavelength expansion, we 
nevertheless view them as approximate equations, with 
domain of validity is C ^ Lp. Hence, in our view, the 
solutions which grow exponentially on a timescale of or- 
der the Planck time or oscillate in time (with significant 
amplitude) at frequencies of order the Planck frequency 
lie outside the domain of validity of the approximation. 
From this point of view, the additional degrees of free- 
dom admitted by the semiclassical equations are merely 
artifacts of the semiclassical approximation. In partic- 
ular, from this point of view the exponentially growing 
solutions to the semiclassical equations are spurious, and 
are not indicative of any physical instability of flat space- 
time. A similar attitude toward the additional degrees of 
freedom admitted by the semiclassical equations would, 
of course, result from viewing the semiclassical equations 
as approximate equations, with unknown higher order 
correction terms of relative magnitude ~ Lp/C^ as dis- 
cussed in Sec. II B above and in Ref. [34]. 



However, two nontrivial issues result from this view- 
point towards the "pathological" solutions, with regard 
to the extraction of physical predictions from the semi- 
classical equations. First, if the domain of validity of the 
semiclassical equations is £ ^ Lp, then since the correc- 
tions to the classical equations which appear in the semi- 
classical equations are of relative magnitude Lp/£'^, it 
might seem that solutions to the semiclassical equations 
cannot accurately describe any situation where the devi- 
ation from a classical solution becomes large. In other 
words, in any circumstance where semiclassical theory 
makes a prediction signiflcantly different from that of the 
classical theory, it might be expected to be highly inaccu- 
rate. (Note that this difficulty would appear to be even 
more severe if one takes the view that the semiclassical 
equations have unknown corrections of order ~ Lp/C^.) 
If that were the case, there would be little point in study- 
ing the detailed properties of solutions to the semiclas- 
sical equations. However, we shall argue in subsection 
IV B below that this is not the case: Semiclassical theory 
should be able to accurately describe phenomena where 
the deviations from classical behavior are locally small, 
but where, nevertheless, long term cumulative effects re- 
sult in very large global deviations from classical solu- 
tions. Furthermore, this ability of semiclassical theory to 
accurately describe such phenomena should the hold even 
in the case of finite N when unknown correction terms 
of order ~ Lp/ C'^ appear in the semiclassical Einstein 
equation. While this viewpoint on the domain of validity 
of the semiclassical Einstein equation is fairly widespread 

e.g., it is commonly assumed that semiclassical gravity 
will give an accurate description of the black hole evapo- 
ration process until the stage where the black hole mass 
approaches the Planck mass — we will devote subsection 
IV B to giving a clear justification for it and to distin- 
guishing this use of the semiclassical equations from using 
them to obtain "approximate perturbative solutions" . 

The second issue that arises is that if the semiclassical 
equations admit spurious solutions, by precisely what cri- 
teria do we determine whether a given solution is "phys- 
ical" or not? This issue has recently been addressed by 
Simon in a scries of papers [34 36] . Simon actually makes 
a number of independent suggestions with regard to the 
semiclassical equations. In this section we will discuss 
all these suggestions in detail. We will argue that, in 
a general context, all of the suggestions have shortcom- 
ings. However, we also will argue that in the special 
case of perturbation theory about Minkowski spacetime, 
Simon's "reduction of order" proposal [36] yields a sat- 
isfactory prescription for extracting physical predictions 
from the semiclassical equations. In the remaining sec- 
tions of this paper, we then will investigate the validity of 
ANEC for solutions to the "reduced order" semiclassical 
equations. 

The above mathematical and physical issues which 

arise from the "higher derivative" nature of the semi- 
classical equations are not unique to this context. In- 
deed, in this regard, the semiclassical Einstein equation 
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is closely analogous to the Abraham-Lorentz equation of 
motion for a classical charged point particle including 
radiation reaction, which also has exponentially growing, 
"runaway" solutions [36]. For completeness and also to 
aid our analysis of the semiclassical Einstein equation, 
we discuss this (much simpler and well studied) exam- 
ple in subsection IV A. In subsection IV B, we explain 
how solutions to the semiclassical equations can be ac- 
curate even when they predict large (global) deviations 
from classical theory. The final two subsections critically 
analyze the "perturbative solutions" and the "reduction 
of order" proposals discussed by Simon [33-36]. 



A. The analogy to radiation reaction of point 
particles 

We begin this subsection with a discussion of the na- 
ture and range of validity of the Abraham-Lorentz equa- 
tion, analogous to our discussion in Sec. II B of the status 
of the semiclassical Einstein equation. It is well known 
that the idealization of classical, charged point particles 
is inconsistent, and that this inconsistency is the source 
of the well known difficulties with the radiation reaction 
equation. More specifically, consider a finite distribution 
of charge with some physical size ^ C. In the nonrela- 
tivistic limit the radiation reaction force on the charge 
distribution will be given by [73] 



■ react 

m 



Ta-FO[Ta(£/c)], 



(4.1) 



where a is the acceleration, m is the mass, q is the charge, 
c is the speed of light and 



T= -(gVmc^). 



(4.2) 



One might think that the "point particle limit" can sim- 
ply be obtained by letting £ — > 0, yielding the usual radi- 
ation reaction force Freact = ^nra, without any unknown 
correction terms. The difficulty with this is of course that 
all physical, finite distributions of charge satisfy 



(4.3) 



assuming only that the electromagnetic self energy ~ 
I a cannot exceed the externally measured mass-energy 
m<? (i.e., that the "bare mass" for a finite distribution 
of charge is always positive). Therefore the limit £ — > 
for fixed mass m and charge q is unphysical. The limit 
£ — > 0, g ^ with (f' j C fixed presumably will exist and 
presumably will yield the unique ratio Freact/?^ given by 
Eq. (4.1) (without correction terms) in the limit. How- 
ever, for fixed mass m and (nonzero) charge g, there will 
always be an ambiguity in the radiation reaction force 
that depends on the structure of the particle, which is of 
fractional magnitude > r/r*, where r* is the timescale 
over which the acceleration is changing. 



It should be noted that the restriction (4.3) and, cor- 
respondingly, the presence of unknown correction terms 
of order O(t^) in the radiation reaction force (4.1) — 
should occur in any theory where classical electromag- 
netic fields are coupled to some other degrees of freedom 
with a continuous distribution of charge. In particular, 
in semiclassical QED, a one-particle, nonrelativistic elec- 
tron state is effectively a finite distribution of charge with 
width £ ^ cr/a, where a is the fine structure constant, 
since the expected value of the current operator will not 
be concentrated at a point but will typically have a width 
of order £. Hence, any "derivation" of Eq. (4.1) for elec- 
trons from semiclassical QED should also give rise to cor- 
rection terms of the order indicated. 

Thus, the status of the equation of motion for a non- 
relativistic, charged particle in a given, fixed external 
electric field E, 



X = — E(x, i) -Frx, 
m 



(4.4) 



is closely analogous to that of the semiclassical Einstein 
equation. First, the domain of validity of this equation 
is limited to the regime r* r, similar to the domain of 
validity £ Lp of the semiclassical Einstein equation. 
Second, as in the case of the semiclassical Einstein equa- 
tion, the small parameter appearing in Eq. (4.4) multi- 
plies a term containing higher-order time derivatives than 
originally appeared in the equation of motion, thereby ef- 
fectively increasing the "number of degrees of freedom" 
of the system. This higher-order time derivative term 
TX is responsible for the existence of so-called runaway 
solutions, which grow exponentially in time. 

Consider now the space of solutions of the radiation 
reaction equation (4.4). If the electric field E is of com- 
pact support in time, the runaway solutions will have 
the property that x oc exp(t/T) for large t. Solutions 
which do not manifest this run- away behavior form a 6- 
dimensional submanifold of the 9-dimensional manifold 
of solutions to Eq. (4.4), since the non-runaway solutions 
have vanishing acceleration at late times (for E of com- 
pact support in time). Indeed, when the electric field is 
bounded above, it can be shown that given any initial 
position and velocity, there exists a unique initial accel- 
eration generating a non-runaway solution. In particular, 
in the special case where the electric field is independent 
of position x, the general solution of Eq. (4.4) for an 
arbitrary initial acceleration ao at time t = is given by 



a(t) = e' 



f-OO 

clo / dsE(rs)e~ 

/•CXD 

-h— / (isE(t + Ts)e" 



(4.5) 



where a = x. It is clear that the choice of initial ac- 
celeration, ao, which makes the first term vanish is the 
unique choice which generates a non-runaway solution. 
Moreover, as is well known, the non-runaway solutions 
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can be characterized as those solutions which satisfy the 
integro-differential equation given in Jackson [73] 

ir=—f ds'E,(x(t + Ts),t + Ts)e-\ (4.6) 

The runaway solutions quickly evolve into a regime 
where the unknown corrections to the radiation reaction 
force will become as large as the radiation reaction force 
itself, i.e., where t* ^ r in the notation used above. 
Therefore these solutions lie outside the domain of valid- 
ity of the Abraham-Lorentz equation and are normally 
deemed to be "unphysical" . It is conventional to take the 
space of "physical solutions" to be the six-dimensional 
space of non-runaway solutions satisfying Eq. (4.6). How- 
ever, these solutions have the unphysical property that 
"prc-acceleration" is required at early times (before the 
electric field is turned on) in order to avoid run-away 
behavior at late times [74]. 

In summary, there arc close parallels between the ra- 
diation reaction equation and the scmiclassical Einstein 
equation. In our discussions in the remainder of this 
section, we will use the radiation reaction equation as a 
simple example and model for the issues that arise. We 
will return to the issue of obtaining "physical solutions" 
to the radiation reaction equation and of the semiclassi- 
cal Einstein equation in the last two subsections of this 
section. However, we first address an important issue 
concerning the accuracy of solutions to such equations. 

B. The physical relevance of solutions to the 
semiclassical equations 

As we have discussed, the equations we are considering 
have only a limited domain of validity. More specifically, 
Eq. (4.4) holds only when r* > t, whereas Eq. (2.6) [as 
well as its linearized version Eq. (3.38)] has the domain 
of validity C ^ Lp. Indeed, the situation with regard to 
Eq. (4.4) is even worse in that there are unknown correc- 
tion terms of order in that equation. [As we have pre- 
viously discussed, similar unknown correction terms also 
would appear in Eq. (2.6) if one justifies that equation via 
the one-loop approximation or takes N to be finite in the 
approximation: sec Sec. II B above.] Since Eq. (4.4) 
differs from the corresponding equation without radia- 
tion reaction only by a term of order t/t* , and Eq. (2.6) 
differs from the corresponding classical equation only by 
terms of order (Lp/£)^, it might seem that Eqs. (4.4) or 
(2.6) can never be valid in any circumstance where they 
predict large deviations from the corresponding classical 
behavior, i.e., in any circumstance where they have the 
potential to make dramatically new predictions. 

Indeed, it might appear that, at best, the only useful 
and reliable information that ever could be justifiably ex- 
tracted from sohitions to c^qiiations like the ones we arc 
considering would be the information contained in "ap- 
proximate, perturbative solutions" to some finite order. 



To explain what is meant by this, let us focus attention 
on the radiation reaction equation (4.4). Suppose that 
we attempt to perturbatively solve this equation order 
by order in an expansion in the small parameter r ap- 
pearing in that equation. Thus, we seek approximate 
solutions of the form 

J 

x^(^) = ^x(^)(^)T^ (4.7) 

where each x^^^ (t) will satisfy a differential equation that 
is second order in time with sources constructed from the 
externally applied fields and from the x^*^) with k < j. 
Note that there is no problem with runaway effects in 
constructing the approximate solutions yij{t), and, fur- 
thermore, for any given order J, an initial position and 
velocity will determine a unique approximate solution x j 
(unlike the situation with the exact solutions) [75] . 

The key question here is: How can we be justified in 
keeping the higher order corrections, x(-'^(t) for j > 2, 
to the solutions, when the equation of motion itself is 
ambiguous at 0(t^)? More specifically, the equations for 
x^^\t) for all j > 1 are completely ambiguous because 
of the unknown correction terms appearing in Eq. (4.4). 
Thus, there is apparently no justification for going be- 
yond the lowest order approximate solution 

xi(i)=xW(i)+rxW(i), (4.8) 

since unknown terms of order unity appear in the equa- 
tion for x^^^ . If we stop with the approximate solution 
(4.8), the difficulties with the existence of additional de- 
grees of freedom and the presence of "pathological" solu- 
tions to the radiation reaction equations would not arise, 
since the perturbative equations arc not problematical. 
However, solutions of the type (4.8) would not be of much 
interest, since - whatever precise form the exact radiation 
reaction equations take - one could not expect Eq. (4.8) 
to be a good approximation whenever rx^^^t) becomes 
of order x'-°\t), since then the unknown higher order cor- 
rections, x^^\t)T^ with J > 1, to the solution should also 
be comparably large. In other words, Eq. (4.8) should be 
a poor approximation whenever radiation reaction has 
any significant effect upon the motion of the particle. 

Completely parallel remarks would apply to the case 
of the semiclassical Einstein equation (2.6) if one justifies 
that equation via the one-loop approximation, so that 
there are unknown corrections of order h?' [76]. How- 
ever, even when we take our viewpoint of justifying this 
equation via the approximation and thereby treat 
Eq. (2.6) as being valid to all orders in h, the situation 
is not significantly different. Although we would be for- 
mally justified in going beyond the first order (in Ti) ap- 
proximate perturbative solution, there apparently would 
be little point in doing so, since the equation is valid 
only when C ^ Lp, and, in that regime, the higher or- 
der perturbative solutions should merely make tiny addi- 
tional corrections to the first order perturbative solution. 
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Again, however, first order approximate perturbative so- 
lutions would be of little interest. 

Simon [34-36] has argued that one is not justified in go- 
ing beyond approximate solutions of the type (4.8), and 
that this therefore solves the problem of the pathological 
solutions [77]. As we now discuss, we disagree with this 
conclusion. With regard to Eq. (4.4) it is indeed true that 
in generic situations the higher order corrections, x^-'^ (t) 
for j > 2, will be sensitive to the unknown higher order 
corrections to the equation of motion, and that there- 
fore only the approximate solution (4.8) will be physically 
meaningful. However, as we now explain, there are situa- 
tions where the higher order corrections x^-'-' for all j > 2 
are not sensitive to the higher order corrections to the 
equation of motion, and where, correspondingly, physi- 
cally meaningful solutions to Eq. (4.4) can be obtained 
which go well beyond the approximation (4.8) and pre- 
dict large radiation reaction effects. Similarly, there are 
situations where Eq. (2.6) may predict large deviations 
from classical behavior even though C ^ Lp everywhere. 
We now explain these comments in detail. 

These situations where the higher order corrections to 
the solutions to Eq. (4.4) are not sensitive to the higher 
order corrections to the equation of motion (4.4) itself 
arise when radiation reaction effects are "locally small" 
but accumulate secularly, so that they become large at 
late times. This occurs, in particular, when the evolution 
timescale is set by radiation reaction. A good example 
is the case of an electromagnctically bound particle in 
a Coulomb field undergoing a radiation-reaction-driven 
inspiral. Clearly the approximate solution (4.8) will pro- 
vide a poor description of the motion once the radius of 
the orbit has shrunk by a factor of two. However, an ac- 
curate description of the motion will be provided by an 
appropriate solution of the exact equation (4.4), provided 
only that the timescale of the inspiral is much larger than 
the orbital period, so that the radiation reaction effects 
are locally small. 

To see this more concretely, suppose that the exact 
equation of motion (including all higher order correc- 
tions) were of the form 



a = a^xt + na + aiT2 a. 



(4.9) 



Here aext is the acceleration due to the externally applied 
electric field, ai is an unknown numerical coefficient of 
order unity, and ti = T2 = t. (We have temporarily 
distinguished the r's that appear in the radiation reac- 
tion acceleration, and in the next order correction to this 
acceleration, to aid the following discussion.) Then it is 
clear that effects that are quadratic and higher order in 
Ti will be important in describing the inspiral. However, 
contributions of order r| to the solution will give rise 
to small corrections to the inspiral of relative magnitude 
< r/xp, where Tp ^ t is the initial orbital period. This 
can be seen by solving the exact equation with ai = 0, 
and checking a posteriori that the Q:iT|a term is always 
a small correction to the equation of motion when eval- 
uated on this solution [78]. 



Indeed, the above type of situation - where non- 
perturbative effects in a small parameter are large but 
can be reliably calculated even though the equation is 
known only to first order in that parameter - actually oc- 
curs quite commonly in physics. A good example is pro- 
vided by Newtonian hydrodynamics. Dissipative terms 
in the hydrodynamic equations normally have an effect 
on the fluid motion that is smaller than the effects of the 
non-dissipative terms by a factor ^ e, where e is the ra- 
tio of a microscopic lengthscale to a macroscopic length- 
scale. Derivations of the hydrodynamic equations from 
statistical mechanics throw away small correction terms 
of order e^. However, effects that are non-perturbative in 
e in solutions to the equations will be meaningful when 
the macroscopic evolution timescale is determined by dis- 
sipative effects. In this case the dissipative terms and the 
higher order corrections are, in effect, boosted from be- 
ing 0{e) and O(e^) respectively, to being 0(1) and 0(e) 
respectively (relative to the non-dissipative terms). 

To illustrate this claim with a simple, concrete exam- 
ple, consider the one dimensional heat equation 



dT 
'dt 



(4.10) 



where a = /^/r, and / and t are some microscopic length 
and time scales. First, we note that it is clear that effects 
which are non-perturbative in the "small parameter" a 
are very significant for the solutions. Indeed, if we choose 
initial data of compac;t support, then at later times the 
temperature T will be nonzero outside the support of 
the initial data. However, the approximate perturbative 
solutions analogous to (4.7) (generated by expanding in 
a) will be nonzero only in the support of the initial data. 
Therefore none of the approximate perturbative solutions 
are even qualitatively accurate; non perturbative effects 
are vitally important. 

Now suppose that there were a correction term to 
Eq. (4.10) of the form -al(Z^/r)5*^/5a;^ where ai is 
an unknown numerical coefficient of order unity. Suppose 
that the initial data is of the form T(x, t = Q) = f{x/C), 
where £ is the macroscopic lengthscale over which the 
initial data varies. We make the following rescaling of 
variables: let x — Cp^ and let t = Ts, where T = {C/1)'^t 
is the macroscopic evolution timescale associated with 
the heat conduction. Then the modified heat equation 
takes the form 



dT 
5^ 



(4.11) 



where e = [l/C)"^ 1, and where the initial condi- 
tions are T{p,0) = f{p). From the form of this equa- 
tion it is clear that the solutions will have important 
non-perturbative contributions from the first term on the 
right hand side, but that the second term (for a\ > 0) 
will be a small correction of order e. 

The situation with regard to the semiclassical Einstein 
equation is closely analogous. Approximate perturbative 
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solutions similar to Eq. (4.7) to any finite order in h will 
not be adequate to describe such phenomena as the evap- 
oration of a black hole over timescales long enough for the 
black hole to lose a significant fraction of its initial mass. 
To describe this process, it will be necessary to consider 
effects that are non pcrtiirbative in h. It should be pos- 
sible to calculate these effects reliably from the semiclas- 
sical equations provided only that TZ ^ Lp everywhere 
in the region of interest, where TZ is the local radius of 
curvature. Even if there were unknown corrections to the 
semiclassical equations that are higher order in h, these 
corrections should be qualitatively unimportant in the 
black hole evaporation process except near the singu- 
larity and near the final moments of evaporation, where 
these unknown corrections become locally large. 

To summarize, given an equation of motion whose 
range of validity restricts a "correction term" in that 
equation to be locally very small, there, nevertheless, can 
be a wide range of circumstances where this equation can 
reliably predict phenomena where this correction term is 
responsible for producing large deviations from the un- 
corrected motion. The approximate, perturbative solu- 
tions (4.7) to any finite order are completely inadequate 
for describing such phenomena. Even if the equation be- 
ing considered is itself valid only to first order in some 
small parameter, it can occur that solutions which are 
non-perturbative in this parameter are physically mean- 
ingful. Therefore, it is of critical importance to have a 
means of determining which solutions to these equations 
should be viewed as "physically relevant" and which so- 
lutions should be deemed to be "spurious" . The next 
two subsections examine two proposals for extracting the 
physically relevant solutions. 

C. Extraction of a preferred subclass of "physical 
solutions" 

A possible method for dealing with the additional de- 
grees of freedom of the modified equations is to iden- 
tify a preferred subclass of the space of exact solutions. 
This can be done for the radiation reaction equation by 
simply discarding the run-away solutions, although the 
remaining solutions have the unphysical feature of "pre- 
accclcration" . However, it is less clear what should be 
done when some of the additional degrees of freedom are 
associated with oscillatory solutions instead of exponen- 
tial solutions, as occurs in the linearized semiclassical 
equations. It is even less clear what should be done in 
the case of nonlinear equations (such as the full semiclas- 
sical equations), where the solutions might not cleanly 
separate into subclasses of the "correct size" on the basis 
of their late and/or early time behavior. 

One proposal for identifying a preferred subclass of 
"physical solutions" is to admit only those solutions 
which are "perturbatively expandable" [33] in the ap- 
propriate small parameter e (where e would be r in the 



radiation reaction case, and h or, equivalently, l/a^, in 
the semiclassical case). By "perturbatively expandable", 
it is meant that the solution can be expressed as a con- 
vergent power series of the form (4.7) (with J = oo), with 
each x^-^^ satisfying the appropriate jth order perturba- 
tion equation. Equivalently, the requirement is that the 
solution should be connected to a solution with e = by 
a one parameter family of exact solutions with parame- 
ter e which is analytic in e. Note that it is essential that 
analyticity in e be imposed, since it should be possible to 
connect every solution at finite e to a solution with e = 
with a one-parameter family which is merely smooth in e 
[80]. 

This proposal would appear to be of the correct char- 
acter, since the perturbative equations have the correct 
number of degrees of freedom, and the "pathological solu- 
tions" do not have analytic behavior in the small param- 
eter e at e = 0. This subsection is devoted to a critical 
examination of this proposal. 

To begin, consider the radiation reaction equation, and 
suppose that the electric field E(x, t) is analytic and in- 
dependent of position. In this case it can be seen from 
Eq. (4.4) that the series generated by solving order by 
order in r is 

= E y^m. (4.12) 

n=0 

This is precisely the expansion in r of the non-runaway 
solution [the second term in Eq. (4.5)]. Thus, when the 
series converges, the two coincide. It seems plausible that 
the "perturbative expandability" criterion also will select 
the non-runaway solutions in the more general case of an 
analytic E which is position dependent. 

However, the criterion of "perturbative expandability" 
fails in the case of smooth but non- analytic E. In partic- 
ular, under the circumstances where Eq. (4.5) is applica- 
ble, it can be seen that no are no solutions which are ana- 
lytic in T. In this case the series (4.12) will not converge. 
Indeed, when the electric field is smooth and of compact 
support in time (and hence non-analytic), each term in 
the perturbative expansion will consist of straight line 
motion both before and after E is "turned on". There- 
fore the summed series, if it converges, must also have 
this property. However, the summed series must also 
satisfy Eq. (4.4), and all exact solutions of this equation 
which have vanishing acceleration at late times will ex- 
hibit pre-acceleration at early times. Consequently, the 
series which attempts to define the "perturbatively ex- 
pandable" solutions cannot converge. 

The criterion of "perturbative expandability" appears 
to fail much more dramatically when (at least some of) 
the additional degrees of freedom are oscillatory in char- 
acter, as is the case for the semiclassical equations. As a 
simple model of this phenomenon, consider the differen- 
tial equation 

(e2^ + l)5(i)=p(i), (4.13) 
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in the limit e ^ 0. Note that the linearized semiclassical 
Einstein equation for individual spatial Fourier modes 
of the metric perturbation is closely analogous to this 
equation with g of the form g{t) = f {t) + WQf{t). The 
general solution of Eq. (4.13) is given by 



9it) = \ j ds sin(|s|)p(i + se) 

+A{e) sin(f/e) + B{e) cos{t/e). 



(4.14) 



The inhomogeneous, first term in this solution can also 
be written as 



fl'inhom(i) = P-V. 
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rP(w), 



(4.15) 



where P.V. means "the principal value of" . This term is 
just the average of the advanced and retarded solutions. 

If p{u!) is of compact support in ui — which is a much 
stronger requirement than p{t) being an analytic func- 
tion of t then it is easy to show from Eq. (4.15) that 
fl'inhom(i;e) IS analytic in e. It follows that when p is of 
compact support, Eq. (4.15) yields the unique "pertur- 
batively expandable" solution to Eq. (4.13). However, 
when p fails to be of compact support in u), it appears 
that there do not exist any perturbatively expandable so- 
lutions to Eq. (4.13). In essence, the fact that one has 
a pole on the real axis of the w-plane in the integrand 
of Eq. (4.15) - which is associated with the existence of 
new oscillatory modes - makes the analyticity behavior 
in e at e = much worse than in the case where the new 
degrees of freedom correspond to exponentially growing 
or decaying modes (i.e., when the poles occur away from 
the real u> axis). Although we have not obtained a com- 
plete proof that no perturbatively expandable solutions 
exist when p(w) fails to be of compact support, we have 
verified that the series 
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71=0 



(4.16) 



which defines the perturbatively expandable solutions 

fails to converge for some simple, analytic, very well be- 
haved choices of p, including Gaussian behavior in t. 

Thus, it appears that the criterion of perturbative ex- 
pandability is of very limited applicability. Even in cases 
where the additional degrees of freedom have an exponen- 
tially growing and/or decaying character, the criterion 
may fail. However, when the additional degrees of free- 
dom have an oscillatory character, it appears that per- 
turbatively expandable solutions will exist only in very 
exceptional cases. 

One might seek some other criterion which would sin- 
gle out a preferred subclass of "physical solutions" . In 
the case where the additional degrees of freedom have 
an exponentially growing and/or decaying character, the 
non-runaway solutions are, of course a natural candidate 
for this preferred subclass [81], although even in this case, 
these "physical solutions" have unphysical features like 



"pre-acceleration" [74]. However, when the additional 
degrees of freedom have an oscillatory character, there 
seems little hope of singling out a preferred subclass of 
solutions on any physical grounds. The difficulties en- 
countered in doing this can be seen in our above model 
(4.13): The issue of picking out a preferred "physical so- 
lution" is essentially equivalent to picking out a preferred 
Greens function for the diff'erential operator appearing 
in Eq. (4.13). Although it is possible to mathematically 
identify preferred Green's functions (e.g., the retarded 
Green's function, the advanced Green's function, or their 
average) there does not appear to be any grounds for ar- 
guing that any one of these is "better behaved" or "more 
physical" than the others. 

Note that the different solutions obtained by choosing 
different Greens functions will all be tangent to the same 
approximate perturbative solution (in the sense of having 
the same derivatives with respect to e at e = 0). They 
will differ by a function which is smooth in e, but which 
is also non-perturbative in e in the sense that all of its 
derivatives with respect to e vanish at e = 0. For exam- 
ple, in our simple model (4.13), the difference between 
the advanced and retarded solutions is just 
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(4.17) 



which is smooth in e as e — > if p{t) is smooth. 

We now turn our attention to a quite different idea: 

the modification of the equations themselves so that all 
of their solutions will be "physical" . 



D. Reduction of order — modifying the original 
equation 

In this section, we analyze the method of "self- 
consistent reduction of order" [36] as a means for ob- 
taining physical predictions from the radiation reaction 
or semiclassical equations. Instead of seeking to identify 
a subset of "physical solutions" to the given equation, 
this approach generates a modified, second order equa- 
tion, which is "as accurate" or "nearly as accurate" as 
the original equation, but whose solutions are all well be- 
haved and can be interpreted as being "physical". The 
idea of reduction of order is quite old — it has been ad- 
vocated in the context of the radiation reaction equation 
by Landau and Lifshitz [82], Teitelboim [83], and Ford 
et. al. [84]. It is also a standard procedure that is used 
in the derivation of post-Newtonian equations of motion 
in classical relativity, see, e.g., Ref. [85]. It has been used 
in the context of classical, higher derivative theories of 
gravity by Bel et. al. [86] , and more recently it has been 
discussed in detail in a wide variety of contexts, and in 
particular advocated in the context of semiclassical grav- 
ity by Simon [36]. 

The justification for this method can be understood 
as follows. We are given an equation of motion which 
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is believed to accurately describe phenomena with suffi- 
ciently large length and/or time scales, e.g., r* ^ r for 
Eq. (4.4) or £ > Lp for Eqs. (2.6) and (3.38). However, 
the given equation (presumably) does not predict even 
qualitatively correct behavior outside of its range of va- 
lidity. Now, generically, any solution to the given equa- 
tion will have some non-vanishing Fourier components 
which lie outside the equation's domain of validity. For 
some solutions, these Fourier components behave in such 
a pathological manner that the entire solution is domi- 
nated by the qualitatively incorrect behavior, as occurs 
for the "run-away" solutions. However, a good remedy 
for this difficulty would be to modify the given equation 
so that it is equivalent — to the desired accuracy — to 
the given equation at large length and/or time scales but 
does not predict any pathological behavior at small scales 
(and, thus, presumably, is at least qualitatively correct 
in this regime). In the situations where it is applicable, 
the reduction of order method achieves this goal. 

The reduction of order algorithm for an ordinary or 
partial differential equation may be stated as follows. We 
start with an equation (or system of equations) for the 
unknown variable x of the general form 



(4.18) 



where r is a "small parameter" . Here we assume that P 
contains terms involving no more that (n—1) time deriva- 
tives of X, but that Q contains terms involving m > n 
time derivatives of x, so that the "small correction", tQ, 
introduces time derivatives of the same or higher differ- 
ential order as appeared in the original equation. To 
apply the reduction of order procedure, we differentiate 
Eq. (4.18) (m — n) times with respect to t, and substi- 
tute the resulting formula for d"^x/dt"^ into the expres- 
sion for Q. We then discard the resulting terms in Q 
which are quadratic and higher order in r. We thereby 
obtain a new equation which is formally equivalent to 
Eq. (4.18) to order r and which has the same general 
form as Eq. (4.18), but for which the term which plays 
the role of Q now contains at most (m — 1) time deriva- 
tives of X. Wo then continue to iterate this procedure 
until the maximum number of time derivatives of x ap- 
pearing in Q is reduced to (n — 1), at which point no 
further reduction of differential order of the time deriva- 
tives can be achieved. The resiilting equation is then of 
the same differential order in time as the original equa- 
tion Eq. (4.18) with r = 0. Thus, we end up with an 
equation which, formally, is "as accurate as" Eq. (4.18) 
to order r, but which does not contain any new "degrees 
of freedom" . Note that this final, reduced order equation 
is uniquely determined by the requirements: (i) that it 
should contain only terms that are zeroth order or first 
order in r, (ii) that it should be formally equivalent to 
Eq. (4.18) to 0{t^), and (iii) that it should be of the 
form (4.18) where the right hand side does not contain 
any derivatives with respect to time of order higher than 
n- 1. 



However, although the reduction of order algorithm is 
uniquely defined for any equation of the form (4.18) for 
a given choice of variables, it should be noted that some 
ambiguities in the algorithm can be introduced by mak- 
ing a T-dependent change of variables: If one introduces a 
new variable y = y{x; r), re- writes Eq. (4.18) as an equa- 
tion for y, and then neglects the terms of order and 
higher, the resulting reduced order equation for y need 
not be precisely equivalent to the reduced order equation 
for X [87]. However, this inequivalence of the equations 
can occur only at order and higher, and, thus, should 
not have an important effect on the behavior of solutions 
in regimes where reduction of order can be justified [c.f., 
the discussion in Sec. IV B above]. 

The radiation reaction equation (4.4) provides a good 
illustration of how this procedure works and of its jus- 
tification. By following the steps described above, one 
obtains the equation [82,83] 



m 



d'E 

E(x, t) + T^(x, i) + r (x • V) E(x, t) 



. (4.19) 



This modified equation of motion is formally equivalent 
to Eq. (4.4) up to order 0(t), and differs from it at or- 
der O(t^). Since the unmodified equation (4.4) has un- 
known corrections at order O(r^), the modified radiation 
reaction equation (4.19) gives a description of the mo- 
tion whose expected accuracy in the regime t* ^ r is 
just as high as that of the original equation (4.4). [In- 
deed, Eq. (4.19) differs from Eq. (4.4) at order 0(r2) 
merely by the term r^a, which is of the same order of 
magnitude as the expected corrections to Eq. (4.4) dis- 
cussed in Sec. IV A above due to the finite size effects 
[88].] However, Eq. (4.19) suffers from none of the prob- 
lems of the original equation. The modified equation of 
motion is second order in time, so there are no "new de- 
grees of freedom" present. All of its solutions are well 
behaved, i.e, there are no runaway solutions nor any pre- 
acceleration effects. We conclude that, in this case, the 
problem of pathological solutions can therefore be over- 
come by adopting Eq. (4.19) as the equation of motion. 

A refined version of the reduction of order method 
can be applied when the original equation is known to 
higher than first order in the small parameter r, as occurs 
in semiclassical gravity when the equations are justified 
via the approximation, as discussed in Sec. II B 

above. To illustrate this refined version, consider, again, 
Eq. (4.18), but in the case where this equation is known 
to O(t^), so that there is possibly an additional explicit 
correction term present of the form t^Q'. We wish the 
reduced order equation also to be valid to O(r^). To 
achieve this, we eliminate the higher order derivatives 
from the term tQ exactly as before, except that we now 
discard only the new terms which are cubic or higher 
order in r. The resulting equation will then be of the 
desired forni at order r, but there will remain a term of 
the form t'^Q, where Q contains higher derivatives. How- 
ever, these higher derivative contributions to Q can then 
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be eliminated by applying the same procedure to Q as 
was previously applied to Q. Clearly, this procedure can 
be generalized to any finite order in r. 

Consider now the application of the reduction of order 
procedure to the linearized semiclassical Einstein equa- 
tion (3.38), where the small parameter is fi. Equation 
(3.38) is an integro- differential equation for the metric 
perturbation rather than a local partial differential equa- 
tion. In general, the reduction of order procedure could 
be ambiguous for such equations, since it may be possible 
to alter the apparent differential order of terms in an inte- 
grand via integration by parts. However, since the right 
hand side of Eq. (3.38) involves only derivatives of the 
linearized Einstein tensor, there is an obvious procedure 
in this case for obtaining a reduced order equation valid 
to order h [or, equivalently, to order = (Lpj C)^, 

c.f., Eqs. (3.38) and (3.45) above]: We merely substitute 
G^^l = on the right hand side of Eq. (3.38). However, 
the resulting equation is too trivial in that it does not 
incorporate any of the effects of the curved spacetime. 
In order to see these effects - and, thus, the dominant 
curvature-related contributions to ANEC at small curva- 
tures and long wavelengths we must go to second order 
in h. To do so, we apply to Eq. (3.38) the above "re- 
fined version" of the reduction of order algorithm to or- 
der . This corresponds to substituting hsat/n for G^^^ 
everywhere on the right hand side of Eq. (3.38), where 

Sob = '^fb^in^- Using Eq. (3.36), the resulting reduced 
order equation is 

kG';^)^{x) = hsab + hih[aAabix) + bBab{x)] 

+ h^{(3Bab{x)+ [ d*x'Hx/L^{x-x') 
I Jm 

X [aAabix') + bBab{x')] I 
+ 0[n^{\nn)% (4.20) 

where Aab and Bab are given by Eqs. (3.33) with G^^^ 
replaced by Sab/n. As discussed in Sec. HIE above, the 
quantity X/ Lp in Eq. (4.20) is a ?i-indcpendcnt constant. 
This modified equation is second order in time and simply 
has the form of the classical linearized Einstein equation 
with a given source. Thus, it has no new "degrees of 
freedom" nor does it admit any solutions with patholog- 
ical behavior. Furthermore, the exact solutions to this 
reduced order equation (4.20) will fail to satisfy the un- 
modified linearized semiclassical equation (3.38) only by 
terms of order 0{h^). Thus, we shall adopt Eq. (4.20) as 
the equation of motion for linearized semiclassical gravity 
in our subsequent analysis. 

Several facts should be noted concerning the above re- 
duction of order of the first order semiclassical equation. 
First, solutions to the reduced order equation (4.20) ac- 
tually correspond precisely to the second order "approxi- 
mate perturbative solutions" of subsection IV B [i.e., the 
approximate perturbative solutions obtained by retain- 



ing terms of 0(1), 0{h\nh), and 0{h)]. This very spe- 
cial situation arises because all terms of order 0{h) which 
involve /i^^^ in the unmodified equation (3.38) are propor- 
tional to G|j^\ and, thus, vanish in the classical limit. In 
more general situations, even for linear equations, solu- 
tions to reduced order equations will differ significantly 
from approximate perturbative solutions. As discussed 
above, in situations where the solutions do differ and 
where the reduction of order procedure can be justified, 
solutions of the reduced order equations should give a 
much better description of physical phenomenon than 
approximate perturbative solutions. 

Second, although we have formally treated h as the 
small parameter, we could equivalently have started from 
Eq. (3.45) and used l/a^ instead: c.f., Eq. (3.47) above 
and associated discussion. A closely analogous equiva- 
lence applies to the radiation reaction equation: If we 
consider a one parameter family of electric fields given 
by E(x, t;a) = a^'^Eo(x, i/a) and define X(t;a) = 
x{at;a), then X satisfies the differential equation 

X=^Eo(X,t) + -X. (4.21) 
m a 

It is clear that reducing order treating l/a as the small 
parameter is equivalent to treating t as the small param- 
eter. 

Third, we can only justify going to order (or higher) 
in the reduction of order procedure in the context of the 
limit. This is because, in the case of finite N, there 
will be unknown corrections to Eq. (3.38) at the same 
order [0(7i^)] as terms that we have retained. The 1/A'' 
limit is still necessary even if we specialize to the situ- 
ation, discussed in Sec. HIE above, that the incoming 
gravitational radiation does not dominate the first order 
metric perturbation. To see this, let us write the un- 
known corrections to the right hand side of Eq. (3.38), 
for finite N, as 

^Labih^')] + ^Kabiu^l:^] + 0{h\ 1/N% (4.22) 

where Lab and Kab are linear but otherwise unknown 
functionals of h^^^ and u)^^^ respectively. Now if we as- 
sume that the incoming gravitational radiation satisfies 
the condition (3.49), then for all solutions of the equa- 
tion, h^^^ cx fi [c.f., Eq. (3.52) above], and the first term 
in Eq. (4.22) can be neglected. However, the second term 
will still be present unless we let N oo. 

Our final remark is that, in the large N limit, the orig- 
inal equation (3.38) is formally known to all orders in h, 
whereas the reduced order equation is valid only to or- 
der h'^. Therefore, the reduced order equation (4.20) is 
slightly less accurate than the original equation. How- 
ever, this slight loss of accuracy is unimportant since - 
for the reasons explained in Sec. IV B above - the effect of 
0{Ti^) corrections should be negligible in the long wave- 
length limit. Note that if more accurate equations were 
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needed, it would be straightforward to iterate the redue- 
tion of order proeedure to obtain an equation accurate 
to any desired order in 7i. In Appendix B below we shall 
effectively carry out reduction of order to arbitrarily high 
order in h. 

We now consider the second order semiclassical Ein- 
stein equation (3.5). As we have previously noted, the 
explicit form of this equation is not known, since the term 
{T^^\h^-^\ hS^^.LJiafi) has not been evaluated. However, 
in subsection III E above, we derived the explicit approx- 
imate form (3.59) of this equation, which is valid for long 
wavelengths and when the incoming gravitational radia- 
tion does not dominate the first order metric perturba- 
tion. The approximate equation (3.59) also has the char- 
acter of possessing higher derivative terms multiplied by 
a small parameter (namely, 1/a^), and has unknown cor- 
rection terms of order 0(lna/a^). Thus, we may apply 
the reduction of order algorithm directly to this equa- 
tion to obtain an equation which should be as accurate 
as Eq. (3.59) at long wavelengths but which has none 
of the pathological behavior at short wavelengths. We 
obtain 

->S^[x^-\x^^'°^]+o((^). (4.23) 

Here x^^'*^^ denotes the retarded solution to the equation 

«G«[x(^-°)] = (rr,4f). (4.24) 

We shall use Eq. (4.23) in our analysis of the validity of 
ANEC to second order in e in Sec. VII below. Note that 
solutions to Eq. (4.23) also coincide with second order 
approximate pcrturbative solutions. 

As we have just argued, the reduction of order proce- 
dure is applicable to our perturbation analysis and, in 
a completely satisfactory manner, it solves the problem 
of the existence of extra degrees of freedom and patho- 
logical solutions possessed by the unmodified equations. 
However, in general, the method of reduction of order 
has some important limits to its applicability, and we 
now briefly mention two of these. 

First, the method is directly applicable only to local, 
ordinary or partial, differential equations, although we 
were able to extend it in a natural manner to the intcgro- 
differential equation (3.38). However, the full, nonlin- 
ear semiclassical Einstein equation is a highly nonlocal 
equation, which is not known to be even of an integro- 
differential type. Thus, it is not obvious if and/or how 
the reduction of order procedure could be applied to the 
full, nonlinear semiclassical Einstein equation. 

Second, although the reduction of order algorithm can 
be applied to any system of local differential equations of 



the form (4.18), in the case of partial differential equa- 
tions, the procedure is guaranteed only to reduce the dif- 
ferential order of the time derivatives, not of the spa- 
tial derivatives. Consider for example the equation in 
Minkowski spacetime 

□ $ = p + eH°^'=VaVfcVc$ (4.25) 

for a scalar field where H"-'"^ is a fixed tensor and e is 
a small parameter. The reduction of order procedure can 
be used to eliminate the third order time derivative of $ 
from the equation, but it does not eliminate the third or- 
der spatial derivatives or all of the third order mixed spa- 
tial and time derivatives. In particular, the resulting re- 
duced order equation is not hyperbolic, and presumably 
would not have a well posed initial value formulation. 
Furthermore, in circumstances where this happens, the 
reduction procedure will, in general, necessitate breaking 
Lorentz covariance, i.e., one will obtain inequivalent re- 
duced order equations by carrying out the procedure with 
respect to different choices of time coordinate. Thus, it 
is only in the happy circumstance — such as in the case 
of Eq. (3.38) — where the reduction of order procedure 
simultaneously eliminates all of the higher order time and 
space derivatives that this procedure is likely to yield an 
equation with satisfactory mathematical properties. 

Fortunately, the above difficulties need not concern us 
here. As we have discussed above, Eqs. (4.20) and (4.23) 
appear to provide a completely satisfactory solution to 
the problem of extracting physical predictions from per- 
turbative semiclassical gravity. The remainder of this 
paper will be devoted to investigating whether ANEC 
holds for solutions to these reduced order equations. 

V. THE GENERALIZED ANEC INTEGRAL 

We now turn to the second of the two principal pur- 
poses of this paper, which is to analyze the positiv- 
ity of the ANEC integral in solutions of the pcrtur- 
bative, reduced order semiclassical equations. Specifi- 
cally, given the metric perturbations h^^ij (x) and h^^^{x), 
we would like to show that the ANEC integral along 
any complete, achronal, null geodesic in the spacetime 

(M, Tjab + ^h^ab '^'^'^^''ab) '^^ non-ncgativc to order e^. How- 
ever, this desired positivity property actually fails, as ex- 
plained in the introduction (see also below). Neverthe- 
less, we do obtain a positivity result involving a trans- 
versely smeared version of the ANEC integral, in which 
the null-null component of the stress tensor is averaged 
transversely to the geodesic as well as along the geodesic. 
This transversely smeared ANEC integral plays a key role 
in our main results. In this section we define a third in- 
tegral which we call a generalized ANEC integral, which 
is an integral over all of spacetime, and which reduces 
to the transversely smeared ANEC integral in a certain 
limit [c.f., Eq. (5.4) below]. This generalized ANEC in- 
tegral is useful as a technical tool in our proofs below. In 
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this section we define the transversely smeared and gen- 
eralized ANEC integrals in general spacetimes. We also 
derive the perturbative expansion in e of the usual and 
generalized ANEC integrals. 

A. Definition of the generalized ANEC integral 

Let 7 be any inextendible null geodesic in an arbitrary 
spacetime, {M,gab)- To begin, fix a smooth, positive 
function 5(x) on R^, with / S'(x) d^x = 1, which de- 
pends only on the magnitude |x| of x. This smearing 
function will control the transverse smearing. Let A be an 
afRne parameter for 7, and denote by A" the null tangent 
vector d/d\. Let V he a fixed point on the geodesic, and 
introduce an orthonormal basis A", C", a.tV, A = 1,2, 
where CCa = 0, X^Ca = -1, and e'^e^gab = Sab- Ex- 
tend this basis by parallel transport to all of 7. Intro- 
duce Fermi- Walker type coordinates x = (A, (, 
a neighborhood M of 7, such that the exponential map 
takes the vector -\- x^e\ at the point a;(A) on 7 to 
the point with coordinates (A, CiS^t)- Then the vector 
field A° = {d/d\Y is a vector field on M which is an 
extension of the tangent to the geodesic. Let x be some 
smooth function which is unity in a neighborhood of 7 
and which vanishes outside M . Fix lengths A, A^ and 
At and define the function 6-^ e C°°(M) by 

e,(A, C, XT) = exp |-i [aVA^ + CVAi] } 

x^5(xt/At), (5.1) 

where xt = (.t^,.t|,). By virtue of the tnmcating fmic- 
tion Xi the function B-y is well defined even where the 
Fermi- Walker coordinates do not exist. We define the 
generalized ANEC integral to be 

7«(A,Ai,AT)= / d^x^f^^Q^{x)Gab{x)X''X\ 
Jm 

(5.2) 

Clearly this quantity depends on our arbitrary choices 
of V, X, C°) etc. However, there are two separate cases 
in which this arbitrariness becomes unimportant. First, 
the limit 

lim lim 7,(A,Ai,AT), (5.3) 

A— >oo Aj:,,At— »0 

when it exists, should be independent of these arbitrary 
choices, and for sufiiciently well behaved stress tensors 
should reduce to the usual ANEC integral. Some of our 
results below will apply to the quantity (5.3). 

The second case in which we can obtain something 
which does not depend on our arbitrary choices of V and 
X is when we specialize the definition to perturbation 
theory about flat spacetime. Here we consider the trans- 
versely smeared ANEC integral 



Is{s;At)= lim lim Is{A,Al,At), (5.4) 

A— >oo Ai,— >0 

where At is fixed at a value of the order of the Planck 
length. In perturbation theory about flat space we can 
choose Af to be the entire spacetime and x to be unity, 
and obtain at each order in e quantities which depend 
only on (i) the length At, (ii) the smearing function S, 
and (iii) the choice of parallel propagated null vector 
along 7 with C"^a = ^1- In subsection VC below we 
derive explicit formulae for these quantities. 

B. Expansion of the ANEC integral 

We now derive the expansion in e of the usual ANEC 
integral; the expansion of the generalized ANEC integral 
will be considered in the following subsection. Suppose 

that "f{s) is a one parameter family of curves on M such 
that 7(e) is a null geodesic with respect to the metric 
Qabi^)- This one parameter family can be represented by 
a map F : R x (— eo,eo) M : (A,£) ^ x°-{X^e), where 
for each e, the parameter A is an affine parameter for the 
corresponding geodesic. To order e, this one parameter 
family of geodesies is characterized by the zeroth order 
geodesic 7 = 7(0) (a null geodesic in Minkowski space- 
time), and by the vector field v'^ = {d/deY on 7. This 
vector field cannot be completely freely specified on 7 
but must satisfy the equation 

A^A^VbVcw" + C*'^ ;,A''A" = 0, (5.5) 

in order that 7(e) be a geodesic to order e. Here the 
tensor C^"^^ is given by 

c(^)",, = -V"/i«-h2V(,/iW;. (5.6) 

Under the gauge transformation (3.24), we have F — > F = 
ip~^ o F, and correspondingly 

=v°-- (5.7) 

The gauge transformation properties (3.31) and (5.7) are 
consistent with differential equation (5.5). 
Now let 

I{e) = f dX Ga6(£)A"(£)A''(£), (5.8) 

where A"(e) is the tangent to 7(5). Let A" = A"(0) be 
the tangent to 7. Then it is easy to show that 

I{s) = sl^^^ + s'^I^'^^ + 0{e^), (5.9) 

where 

/(I) = / G^'(,^A"A^ (5.10) 

and 
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7(2) 



(1) 



(5.11) I?^ = [ e% e(0)(.)Gi^«A^ 



(2) 

Here denotes the Lie derivative, and GJ^^ is the com- 
plete second order perturbation in the Einstein tensor, 
given by 



G 



:2) 



(5.12) 



The second term in Eq. (5.11) can be thought of as re- 
flecting the fact that the metric perturbations cause a 

change in the geodesic. Note that the quantities J^'^-' and 
J^^-' are gauge invariant (in the "active" sense in which 
we are considering the gauge freedom here) , but that J*^^^ 
does depend on the gauge covariant vector field v"" on 7. 



C. Expansion of the generalized ANEC integral 

Consider now a corresponding analysis of the pcrtur- 
bative expansion of the generalized ANEC integral (5.2). 
For each finite e, the construction of the coordinate sys- 
tem described in Sec. VA yields scalar fields \{e), ({e) 
and ctt(£) = d^^x^x^ on the spacetime {M,gab{£)), 
and also the vector field A"(£) = {d/d\Y{e). Note that 
although these objects arc defined in terms of an algo- 
rithm to obtain a coordinate system, they arc themselves 
coordinate-independent scalar and tensor fields. Their 
domain of definition is however restricted to some neigh- 
borhood of the geodesic 7. They can be expanded as 

A(£) = A(")+£A(i)+0(£2) 
C(£)=c(°)+£C(i) + 0(e2) 

A«(e) = A(°)"+eA(i)" + 0(e2), (5.13) 



where the expansion coefficients on the right hand side 
are defined on all of M. Similarly the volume 4-form 
can be expanded as eabcd(.e) = e^°),^ + ^eiVcd + O(e^), 
where e^^^ = h^^^'^e^^^ /2. Inserting these expansions into 
Eq. (5.2) yields 



/,(£) = £/«+£^/f) + 0(£3). 



(5.14) 



where 



7W(A,A^,At)= / ei°2/ eW(a;)G«(a.)A«A\ (5.15) 

Here Q^\x) is given by the expression (5.1) with x re- 
placed by 1, and where the arguments of the expression 
(5.1) are A = X^°\C = C^^K etc., the inertial coordinates 
on Minkowski spacetime. 
The expression for /i^'' is 



+ 



/ 4de/GiV [2A(°) "^A'l) "6(0) + A(°) «A(°) "6 W 
Jm 



+ l/j«e;^(O)a;^(O)60W]. 



(5.16) 



Since the smearing function 5 in Eq. (5.2) depends only 
on the magnitude ax = of x^-, the function G^^^ that 
appears in Eq. (5.16) is given by 



eW(x) 



dQ 



dC 



(5.17) 

As in the previous subsection, the quantities /i^' and 
(2) 

Is are gauge invariant. We now specialize to a partic- 
ular choice of gauge, which is just that associated with 
the coordinate system described in Sec. VA. Specifi- 
cally, suppose that we are give a 1-parameter family of 
metrics gab^s) and a fixed choice of gauge for each s. 
Now apply the diffeomorphism (p^ given by identifying 
the coordinates A, (, x^ in the spacetimes (M, 5afe(0)) and 
{M.gabis)) [89]. This yields a choice of gauge in which 
we have 



AW- = A«=C(^)=a«=0. 
The expression for l'^^ in this gauge reduces to 

Jm 



(5.18) 



(5.19) 



In this gauge the vector field v'^ described in Sec. VB 
vanishes. Moreover it is straightforward to show that 
h^^^ vanishes along the geodesic 7, and that /i^^'a^A'' = 
throughout M . These consequences of our gauge choice 
will be used in Sec. VII A below. Finally, we note that in 
the limit A ^ 00, Al, Ay ^ of no transverse smearing, 
the expression (5.19) reduces to the previously obtained 



expression (5.11), since h 
gauge. 



(l)e 



vanishes on 7 in our chosen 



VI. THE ANEC INTEGRAL IN FIRST ORDER 
PERTURBATION THEORY 

In this section we establish the results concerning the 
first order contribution to the ANEC integral which were 
discussed in the introduction. We start in subsection 
VI A by characterizing the precise class of incoming states 
we are considering. In subsection VI B we obtain the so- 
lutions of the reduced order semiclassical equation (4.20). 
We derive a general formula for the first order perturba- 
tion (5.15) to the transversely smeared ANEC integral 
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for these solutions in subsection VIC. Finally, we show 
that this vanishes for pure incoming states in subsection 
VI D, and that for mixed states it is positive in the long- 
wavelength limit in subsection VIE. 

In Appendix B below we consider a specific subclass of 
solutions of the original, unmodified semiclassical equa- 
tion (3.4), given by the use of half advanced plus half 
retarded Greens functions. We show that for these so- 
lutions, in the region A > Acrit of parameter space, a 
transversely smeared ANEC integral vanishes for pure 
states and is always strictly positive for mixed states, 
even outside of the long wavelength limit. Although, as 
discussed in Sec. IV C above, we do not attribute any pre- 
ferred physical status to these solutions, expanding this 
result in \jo? and Xnajo? provides an alternative demon- 
stration of the results of this section for solutions of the 
reduced order equation in the long wavelength limit, at 
least for A > Acrit- 



A. Characterization of incoming states 

In the remainder of the paper we do not deal with 
fully general states. Rather, we restrict attention to a 
subclass of states whose two point functions have suit- 
able differentiability properties and fall-off properties at 
infinity when restricted to spatial slices. More specifi- 
cally, let S denote the hypersurface i = in Minkowski 
spacetime, and let V be the class of smooth functions 
on S X E all of whose derivatives are L-'^ on E x E. We 
consider states whose two point function perturbations 
i^/jj^-* (x, i; x', i') and Fl^^\-x.,t;x' ,t') and associated time 

derivatives dF-^^^/dt, dF^l^/dt', and d'^F^l^ /{dtdt'), for 
J = 1, 2 all lie in V when restricted to E x S. Our choice 
of class of states is dictated mostly by convenience and is 
not the most general class of states for which our results 
are valid; however it is a sufficiently large class of states 
to be interesting. 



We can express Jin = eF- 



(1) 



as 



Jin(x,ar') = j rf^k j d3k'/(k,k'). 



+ c.c. 



(6.1) 



Here as in Appendix C, bold face vectors are spatial, 
three dimensional vectors, while k = k'^ denotes a 4- 
vector. Also it is understood that k°- — (k, o^k) where 
Wk = |k|. Equation (6.1) essentially defines the func- 
tions / and g as suitable complex linear combinations of 
the spatial Fourier transforms of the four functions i^in, 
dFiJdt, dFin/dt' and d^Fi^/{dtdt') restricted to E x E, 
where E is the Cauchy surface t = 0; see Eq. (C3) below. 
Note that the functions / and g are formally related to 
the conventional creation and annihilation operators aj{^ 
and fik by 



/(k,k') oc (akCk') 
g{k,k') oc {al,ak)- 



(6.2) 



Note also that the part of the two point function that is 
purely positive frequency or purely negative frequency is 
given by the function /, while the "mixed-frequency" part 
(the part that is positive frequency with respect to 
one variable and negative frequency with respect to the 
other) is given by the function g. 
We expand / and g as 



/ = s/(i)+e2/(2)+0(^3) 



e5«+eV''+0(£^). 



(6.3) 



In the terminology introduced in Sec. IIIB above, we 
refer to the first order perturbed two point function F^^^ 
as "pure" if and only if g^^'^ = 0; otherwise F-j-^^ is said 
to be "mixed" . 

The functions / and g cannot be chosen arbitrarily 
but must satisfy the positivity condition (2.16), which 
for Fock space states is just the statement that 



($(u)$(u)t) > 



(6.4) 



for all complex smearing functions u. For general, alge- 
braic states cjin the positivity condition is equivalent to 
the positivity of the operator 



Ma 



9 f 

f* g* + J 



(6.5) 



where J(k, k') = ?i(5"^(k— k')/(167r'^a;k), in the sense that 
for all Uj e C^(M), 1 < j < 2, 

j d^k j dX Wi(k,a;k)*My(k,k')wj(k',a;k') > 0. (6.6) 

Thus, M is a positive operator on L^(R^) i^(R^), 
which implies, in particular, that 5 is a positive operator 
on L2(r3). 



j d^k j dX 



g{k,k')u{k)u{k')* > 0. 



(6.7) 



For Fock space states, Eq. (6.7) is just the statement 
that ($_|_(u)'I'§_|_(m)) > for any test function u. Note 
also that the corresponding classical positivity condi- 
tion, which requires that the two point function be the 
expected value of <I>(a;)$(j/) with respect to some posi- 
tive measure on the space of field configurations, is the 
stronger condition that 



g f 

r 9* 



> 0. 



(6.6 



We now insert the expansions (6.3) into the positiv- 
ity condition (6.5) and expand order by order in e to 
determine the restrictions on the incoming state pertur- 
bations. We obtain at first order that 



30 



7(1) > 



(6.9) 



in the sense of Eq. (6.7), while /(i) can be chosen arbi- 
trarily. At second order we obtain the restriction on g^"^^ 
that 



"P 0' 






/«] 




P 0" 


1 




/(I)* 


J 




1 



>o, 



(6.10) 



where P is the projection operator onto the kernel of g^^). 
In particular, when g^^^ = this reduces to 



(2) 



9 
/(I) 



/(I) 
J 



> 0. 



(6.11) 



Next, wc introduce an alternative, convenient set of 
coordinates on the light cone in momentum space. Recall 
that we have defined inertial coordinates A, Xrp ^ Xrp OH 
Minkowski space M, and also a null orthonormal basis 
A°, ej, 63, where A° is the tangent to the zcroth 
order geodesic 7. Thus, = AA" +xt, where = 
x^eA^ A = 1,2. We introduce corresponding coordinates 
7, /3, on momentum space such that 

r =7A"+/3C" + kT. (6.12) 

In these coordinates the positive light cone volume ele- 
ment can be written as 

2wk 



dp 



7 - 



dj. (6.13) 



Therefore solutions of the wave equation can be rep- 
resented by functions of /3 and of the two dimensional 
vector k.T. These null coordinates on momentum space, 
which are specially adapted to the given nidi geodesic 7, 
will be useful throughout our computations below. 

The first order two point function F-^^^^^ has an expan- 
sion analogous to Eq. (6.1) but with / and g replaced by 
/(I) and g^^\ We can rewrite this expansion in terms of 
the coordinates introduced above as 



/(/3,kT;/?',k^)e^'=-e*'='-' 



+5(/3,kT;/3',k^)e*-e 



+ C.C. 



Here / = Wk^k'/'^), 9 = ^^k^^k' g^'^\ and 



Wk 



2/3 



(6.14) 



(6.15) 



It it understood in these equations that /c = fc° is given 
by Eq. (6.12) with 7 = k|,/(2/3). Finally, as briefly dis- 
cussed in Appendix C, our assumed regularity properties 



on the incoming state imply that / and g arc continuous 
as functions of (/3, k^, /?', k^) and satisfy, for any integer 
N, 



max 



< 



C 



N 



(6.16) 



for some constant Cn, where Wk is given by Eq. (6.15). 

B. Solutions of the reduced order equation 

In Appendix C we show that for states in the above 
class the stress tensor (3.43) which acts as a source in the 
semiclassical equations is an tensor field on Minkowski 
spacetime. Consequently its spacetime Fourier transform 
exists as an function. In Appendix C wc show that 
the Fourier transform Sab is actually continuous every- 
where away from the light cone, and is bounded every- 
where except for an (integrable) divergence at the origin 
in momentum space. Thus, we may use Fourier trans- 
form methods to solve the (reduced order) semiclassical 
equations. 

The reduced order semiclassical equation (4.20) ex- 
presses the linearized Einstein tensor G^^j^[h^^^] in terms 
of the source tensor (3.43). In our analysis below we 
shall not need to solve (4.20) for the metric perturbation 

^ab'^ it will be suSicient to work directly with the lin- 
earized Einstein tensor. Wc now take the Fourier trans- 
form of Eq. (4.20). We use the following formula given 
by Horowitz [48] for the Fourier transform of the distri- 
bution Hx 



^A(fc) = -27r [In A^lfc^l 27 - 1 



i7rsgn(fc°)e(-fc2)] , 
(6.17) 



where 7 is Euler's constant and O is the step function. 
Here and below tildes denote Fourier transforms. Using 
Eqs. (3.33) and (3.36) we obtain 

KG«(fc) = Si{k)Sab + S2{k){kakb - Vabk^)^' , (6-18) 



where 



Si{k) = 1 + IGnaL^pk^Hxik), 



and 



S^ik) = (Stt) 



{-a + 2b)LlHx{k)+2pLl 



(6.19) 



(6.20) 



We have explicitly included factors of Gh = in these 

formulae. The tensor Sab is the spacetime Foiirier trans- 
form of the source tensor (3.43) discussed above. Note 
that Eqs. (6.18) - (6.20) could also be obtained by ex- 
panding the exact solutions (A5) (A8) given in Ap- 
pendix A in powers of h and hlii{h) (assuming A oc 
as discussed in Sec. HIE above). 
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Next wc rewrite the source tensor Sab in terms of the 
regularized two-point function of the incoming state. Us- 
ing Eqs. (1.13), (2.4) and (6.1), we find that 



SabQ) 

(2 



f^^\k,li')aab{k,k')S^{k + k' - I) 



+g^^\k,-k')aab{k, -k')S\k -k' -l) + c.c. 



The existence of the integral (6.25) follows from the reg- 
ularity properties of the functions f'^^^ and g'-^^ discussed 
in Appendix C: they are continuous away from the light 
cone, and have an integrable divergence cx l/(a;kWk') at 
the origin. 

Now by combining the alternative representation 
(6.14) of the two point function with Eqs. (6.25) and 
(6.21) (6.26), we obtain 



where 



lim /(i)(A, Al, At) = 



(^ab{k, k') 



{2^-l)k^ak'b) + {--2^)r,abkck' 
Hikah + k'X)- 



dp 

T 



dfT 



d'^k'j.f{k,k')J{k,k') 



(6-22) x\/2^Aexp [-(/3 + /3')'AV2] S [(kr + k'T)AT] 



Inserting Eq. (6.21) into Eq. (6.18) yields the null-null 
component of the linearized Einstein tensor 

KG(l)(0a6A"A'' = / d^k / d3k'/(l)(k,k')J(fc,fc')'5^(fc + fc' 



d^ 

T 



d'k. 



dfr_ 



d,^k'j-g{k,k')J{k,-k') 



x\/2^Aexp [-{j3 - l3'fh?/2] S [(ky - k^)AT] 
-^c.c. (6.27) 



+g^^\k,k!)J{k, -k')6\k -k'-l) 

+C.C. 

Here the function J is given by 
:^i^=5i(fc + fc')|(2^-l)(A.fc)(A.fc') 



Note that Eqs. (6.15) and (6.16) imply that for kr ^ 0, 
(6.2S|ie integrands vanish more rapidly than any power of /3 
as /3 — > 0, thereby assuring convergence of the integrals 
in Eq. (6.27) despite the factors oi 1/(3 coming from the 
light cone volume element. 



H [(A • kf + (A • k'f 

+(1 - 60(A -k + X- k'f{k ■ k') S2{k + k'). (6.24) 

By using the technicjucs of Appendix C, it is possible to 
show that the linearized Einstein tensor (6.23) has the 
same regularity properties as were proved in Appendix C 
for the source tensor (3.43). 

C. The first order ANEC integral: general formula 

We now calculate the first order perturbation to the 

generalized ANEC integral (5.15) and express it in terms 
of the functions Z'-"'^-' and g^^' characterizing the incoming 
state Wjn^ Combining Eqs. (5.15) and (6.23) yields that 



D. Pure incoming states 

We now show that in the limit A ^ cxd, the pure fre- 
quency contribution to the ANEC integral [the first term 
in Eq. (6.27)] vanishes. This can be seen from the fact 
that in the limit A ^ oo, the exponential factor in this 
term becomes i5(/3 -|- (3') and that therefore the entire ex- 
pression vanishes. Therefore, when £f*^^^ = we have 



lim lim 7^i)(A,Ai„At) = 0, 

A^oo At.Ar— >0 



(6.28) 



/(I) = d^kj rf^k' 



/W(k,k')J(fc,fc')e(°)(/e + fc') 



+g^^\kX)J{k,-k')m{k-k') 



+ C.C. (6.25) 



Here Q'^^ is the Fourier transform of the smearing func- 
tion (5.1): 



e(°)(7, /3, kr) = V2^K exp ^ [H^K^ + 7'Ai] 



X 5(kTAT). 



(6.26) 



and thus whenever the usual ANEC integral exists, it 
must vanish. 

There is a simple, intuitive way to understand this re- 
sult. The ANEC integral is the integral along a line in 
position space, and thus becomes the integral over a hy- 
perplane (the hyperplane X^ka = 0) of the Fourier trans- 
formed Einstein tensor G^^^ in momentum space. Now, 
for states whose first order perturbed two-point function 
is pure (i.e., for which 5^^^ = 0), the linearized Einstein 
tensor has support inside and on the light cone in mo- 
mentum space. This can be seen from Eq. (6.23) and 
the fact that if k and k' are future pointing null vectors, 
then A; -|- fc' is a future pointing null or timelike vector. 
Therefore, the only possible contribution to the ANEC 
integral must be concentrated on the null line fc" cx A'' in 
momentum space. Although this line is of measure zero, 
a priori there still could be a non-vanishing result since 
the (Fourier transformed) Einstein tensor could have a 
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distributional component on the line. However, the ar- 
gument in Appendix C shows that there is no sueh distri- 
butional component of the linearized Einstein tensor, so 
the ANEC integral must vanish, as our calculation above 
shows explicitly. 

As explained in the introduction, the result (6.28) that 
the ANEC integral vanishes to first order for incoming 
states that are pure to first order is one of the key results 
of this paper. It eliminates the counterexample to ANEC 
given in Ref. [18]. Moreover, the vanishing of the ANEC 
integral to first order is a necessary condition for the 
ANEC integral to be positive generally; any nonzero first 
order contribution for pure states could be arranged to 
be negative by choosing the sign of the first order state 
perturbation appropriately. Moreover, in this situation 
any transverse smearing could not help. 

The above rc^sult applies for solutions to the reduced 
order equation which are accurate to 0{h^). However, 
it is straightforward to extend this result to all orders in 
h, that is, to solutions of reduced order equations which 
are accurate to higher order in h. The only difference 
is that the functions and ^2 appearing in Eq. (6.24) 
are slightly altered, which does not affect the argument. 
More precisely, these functions are replaced by expan- 
sions to the appropriate order in h and ?i In ?i of the func- 
tions 1/Fi and F2/{FiFs), as can be seen from Appendix 
B below. 

Finally we remark that a limiting case of the above 
result in which the backreaction is dialed to zero is just 
the fact that the ANEC integral in Minkowski spacetime 
of a matrix element of the form {0\tI^^"^ \'>P)^^^ must vanish; 
see Eq. (3.18) above. This fact can also be deduced from 
the result previously established by Klinkhammar [16] 
and by Wald and Yurtsever [18] that the ANEC integral 

of the expected value of T^l^ in Minkowski spacetime is 
nonnegative for a large class of states. 



E. Mixed incoming states 

We now turn to the situation where we allow an arbi- 
trary, mixed incoming state. Using Eqs. (6.17), (6.19), 
(6.24) and (6.27) we find that 

7W(At)= lim lim /^^^(A, A^, At) 

A— >cx) Ai,— >0 



-I 



Here lu^ = I / {?,2tt^ aLf) , X = Acxp(7 - 1/2), and 7 is 
Euler's constant. Note that the continuity of g and the 
fall-off property (6.16) is sufficient to guarantee the exis- 
tence of the integrals (6.29) and (6.30). Also, note that 
there is no longer any dependence on the curvature cou- 
pling in Eq. (6.29), due to cancelations in Eq. (6.24) when 
k' = -k. 

The function K in Eq. (6.29) is essentially the fac- 
tor 5*1 that appears in Eq. (6.18) (and is also related by 
an expansion in h to the Greens function 1/Fi that ap- 
pears in Appendix A). The second term involving ^2 
in Eq. (6.18) does not contribute to the ANEC integral 
(6.29), because it has a tensorial structure in momentum 
space proportional to 



kakb Tjabk ■ 



(6.32) 



In the original Minkowski space coordinates this corre- 
sponds to the differential operator (P/dX^, which gives 
rise to a total derivative in the ANEC integral and gives a 
vanishing contribution. More precisely, the contribution 
to the first order perturbation I^^^ vanishes identically, 
and the the contribution to the first order perturbation 
of the generalized ANEC integral /i^^ vanishes once the 
limit A ^ 00 is taken. 

The formula (6.29) has a simple physical interpretation 
in terms of an ANEC integral computed in flat spacetime 
without backreaction, as we now describe. Let 7' denote 
the null geodesic in Minkowski spacetime obtained by 
displacing the geodesic 7 transversely by an amount — xy. 
Define 



If{xt) 



7' 



(6.33) 



This is just the ANEC integral obtained from the incom- 
ing state without including the effects of backreaction, 
i.e., in the test field approximation, and with no trans- 
verse smearing. It is independent of our choice of (, and 
hence it is a function on the two dimensional space of 
vectors perpendicular to A", where two vectors are iden- 
tified if they differ by a multiple of A''. Using Eq. (6.29) 
in the limit where K ^ 1 and 5 — > 1, and applied to a 
transversely displaced state, we obtain 



If{^t) 



(27r) 



2 T{AkT)K{AkT)S{KTAkT), (6.29) 



'£_ glAkT-XT 



^(AkT), 



(6.34) 



where 



i.e., If is just the Fourier transform of the function fi. 
Therefore we can rewrite the formula (6.29) as 



1 fi-TT"^ r r 
jr(AkT) = ^J^ dp J (fkr g{0, ky; /?, ky + Aky), 

(6.30) 



and 



^(Akr) = 1 - In Ak^) . (6.31) 



J(i)(At) 



lF{AkT)KiAkT) 

xSiAxAkr). 



(6.35) 



In other words, the ANEC integral with backreac- 
tion is just the test field ANEC integral /f(xt) in 
Minkowski spacetime convolved against the smearing 
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function 5(xt/At), and against the distribution K{'x.t) 
obtained from Kilix) by an inverse Fourier transform. 

Note that it follows from the analysis in Appendix D 
that the distribution K{:x.t) is given by a smooth func- 
tion away from the origin = 0, but not at the origin. 
However, the convolution (K o S){xt) is a smooth func- 
tion for all XT for our choice of smearing function given 
by Eq. (6.40) below. We also remark that Eqs. (6.15), 
(6.16) and (6.30) imply that the function |kTr:^(kT) 
is for any A'", which by Eq. (6.24) shows that the 
"test field ANEC integral" function /^■(xt) is smooth. 
This fact will be important in our analysis in Sec. VI E 2 
below. 



1. Nonnegativity of the ANEC integral in Minkowski 
spacetime 

From previous analyses by Klinkhammar [16] and by 
Wald and Yurtsever [18], it is known that the test field 
ANEC integral If{xt) is always nonnegative for a large 
class of states. This result forms a key element in our 
proof below that the smeared ANEC integral (6.35) is 
always nonnegative for suitable choices of the smearing 
function S. We give a short proof of the result here, in 
order to lay the foundations for later analyses. 

First we give a motivational non-rigorous argument, 
which applies only to states in the usual Fock space. Use 
the decomposition (3.21) of the field operator into pos- 
itive frequency and negative frequency parts. Then we 
obtain from Eq. (2.4) that, up to total derivatives with 
respect to A, 



: fi°) : A^A*- 

ao 



($'_)2 + ($^)2 + 2^'A+, (6.36) 



where primes denote derivatives with respect to A. The 
colons on the left hand side denote normal ordering. The 
first two terms on the right hand side integrate to zero 
when we integrate along the geodesic because this picks 
out the zero frequency part, and the last term is a man- 
ifestly nonnegative operator. Hence the ANEC "opera- 
tor" is nonnegative. 

We now give a rigorous proof of the positivity of 
/f(xt); it follows directly from our general formula 
(6.34) for If and the positivity condition on the two point 
function expressed in our coordinates (/3, kr). From 
Eq. (6.9), choosing u to be a suitable Gaussian and taking 
an appropriate limit allows us to deduce that 

J (fkT J d^k^g{Po,-kT;Po,K) >0, (6-37) 

for all /3q. Therefore, from Eqs. (6.30) and (6.34) the 
quantity If{0) is nonnegative. It is clear that the same 
is true for If{xt). 



2. The long wavelength limit 

We now specialize to the long wavelength limit dis- 
cussed in Sec. Ill E above. We also assume that g^-'^^ ^ 0, 
i.e., that the first order perturbed two-point function is 
mixed. First we show that the unsmeared ANEC in- 
tegral may be negative, and then show that for suit- 
able choices of the smearing function S, the transversely 
smeared ANEC integral is always positive (not merely 
nonnegative) in the limit £/Lp oo. Note that it would 
be inconsistent to analyze the solutions to the reduced 
order semiclassical equation (4.20) outside of this limit, 
as correction terms of 0{Lp/£^) were thrown away in the 
derivation of these equations. 

If we assume an incoming state of the form (3.41), then 
the test field ANEC integral varies as 



/i^(xT;a) a ^iFixr/a), 



(6.38) 



where If is the test field ANEC integral of the state 
W;^'. Consequently the Fourier transform scales as 

/F(AkT;a) = a~^/F(Q;AkT). Substituting this into 
Eq. (6.35) and making a change of variable in the in- 
tegral yields 



n^\Ar)=a-'J'^^F{A^.r)K 



X S 



(27r)2 
ATAkT 



Ak7 



a 



We now choose the smearing function to be 
5(xt) ^ 



(6.39) 



(6.40) 



1 + |xt|4" 

Its Fourier transform has the expansion for small kT 

5(kT) = 1 + uo]^ ln(|kT|) + i/ikl + 0(|kT|^), (6.41) 

where > [90]. Substituting Eqs. (6.31) and (6.41) 
into Eq. (6.39) and expanding in l/a^ yields that 



K7«(AT) = ^+i?— + 



A 

~3 



a 



,lna 



C_ 



O 



(lna)2 



where 



A = 



and 



B = i.o(A|-A2 t)(V?^^F)(0), 



AT.crit — StT 




(6.42) 



(6.43) 



(6.44) 



(6.45) 



Note that the sign of the coefficient a in the formula 
(2.10) for the anomalous scaling of the stress energy de- 
termines the sign of the coefficient B without smearing. 
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As previously mentioned, a is positive for the scalar field 
we arc considering, for all values of the curvature cou- 
pling [c.f.,Eq. (3.34) above], and is also positive for neu- 
trino and Maxwell fields [48]. If a had been negative, 
the coefficient B would have been positive without any 
transverse smearing. 

As we have just discussed, the leading coefficient (6.43) 
is always nonnegative. However, as we now show, there 
do exist states for which /f(0) vanishes, and this opens 
up possibilities for violations of ANEC. Let /i(/3, ky) be 
a function on the positive light cone satisfying 




d^kTlMAMI' =1, (6-46) 



and let \a) be the one particle state 

\a} = l^jj d^^T V^KP, kr) 4|0). (6.47) 

Choose the the incoming state perturbation Q^^^ to be 
that given by the density matrix perturbation 

=-|0)(0| + |a)(a|. (6.48) 

Then the function g is given by 

g{(3, kT;P', k^) (X h{P, kr) h{(3' , k'^)*. (6.49) 

If we now choose any smooth function v{xt) of compact 
support, and choose /i(/3, ky) = hi{(3)v{kT) for some 
suitable hi, then from Eqs. (6.30), (6.34) and (6.49) we 
find that 

Ipiy^T) oc lt^(xT)l^. (6.50) 

Hence in particular we can choose a state which achieves 
If{0) = 0. _ 

Now when A = /f(0) = 0, then the quantity V|Jf(0) 
is always positive or zero. This is because Jf(xt) > 
always, /f(0) = 0, and Ip is smooth as discussed above. 
It is clear that we can find states for which /f(0) = but 
VtJf(O) ^ 0. Therefore, we find that (i) in the limit of 
no smearing [i.e.. At 0, S{xt) (5^(xt)], B can be 
negative and therefore for sufficiently large a, the ANEC 
integral can be negative, and (ii) From Eq. (6.44), when 
the transverse smearing length is larger than the critical 
length (6.45), B is always nonnegative. 

So far we have shown that the smeared, first order 
ANEC integral is always non-negative to order 1/a^ be- 
yond leading order. However, from Eq. (6.50) it is clear 
that there exist incoming states for which B = 0; for 
instance one can choose the function w(xt) to vanish in 
a neighborhood of the origin. When A = B = 0, the 
expression for the coefficient C is 

C = ii.o(A| - A|,„h) / (0 ln[A|k2 ] kl J^(kT). 

(6.51) 



We now shall show that when A = i? = 0, for A^ > 
ATjcrit we have C > 0, with equality holding if and only 
if 5(1) = 0. 

To prove this result, we note first that ii A = B = 0, 
then we have /_f(0) = V^If{0) = 0. However, Ip is a 
smooth, nonnegative function, and hence it follows that 
all the derivatives of Ip at the origin up to and includ- 
ing third order vanish. In Appendix D, we show that 
this implies that the coefficient C is always strictly posi- 
tive unless the function Ip{xt) vanishes identically. We 
now show that /^(xt) cannot vanish identically unless 
g^^"^ ~ 0, i.e., unless the perturbed two-point function is 
pure to first order. Thus, in the mixed case, in order 
to establish positivity of the smeared ANEC integral in 
the long wavelength limit for nearly flat spacetimes, it is 
not necessary to continue the expansion (6.42) to higher 
powers in 1/a, nor is it necessary to go to second order 
in £. 

To prove that Ip{xt) vanishes identically if and only 
if g^^^ = 0, we start by writing 

/.oo 

Ip{kt) = / dpip{KT,l3) (6.52) 
Jo 

where Ip{xt, f3) is defined by Eqs. (6.30) and (6.34) with 
the integral over /3 omitted. Moreover Eq. (6.37) implies 
that Ip{xT,P) > always, and since If{xt) = Ip{xt) = 
we deduce that Ip{xt,P) = for all xy and (3. Hence 

j cfkr g{(3,kT;P,kT + Akr) = (6.53) 
for all Aky. 

Next we obtain a canonical representation for the func- 
tion g{k,k'). Let /C be the measure space (R^,d/u), with 
measure d/i = dfid^kx- Then g is a continuous complex 
function on /C x /C which satisfies g(k,k') = g{k',k)*. 
Moreover Eq. (6.16) implies that g is both and 
on K. X JC. The property of g together with its pos- 
itivity property (6.37) implies that it defines a positive, 
compact, self adjoint operator G on L^{1C) [91]. Hence 
by the Hilbert-Schmidt theorem, there is a complete or- 
thonormal basis ipn of L'^{1C) such that 

00 

g{k, k') = KVn{k)^n{k'Y, (6.54) 

for some A„ > 0, where the convergence is in the operator 
norm on the space of bounded operators on L^(/C). Now 
inserting the decomposition (6.54) into Eq. (6.53) and 
specializing to Ak^ = we find that 

00 „ 

Y^Xn d^kr \<fin{0, kr)!^ = 0. (6.55) 

n=0 

Since A„ > for all n, wc obtain that g = 0. 

By combining all the results in this subsection we 
find that for states whose perturbed two-point function 
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is mixed, the transversely smeared ANEC integral for 
At > Axjcrit will always be strictly positive for suffi- 
ciently large a, for solutions of the reduced order, first 
order semiclassical equation (4.20). 



VII. PURE INCOMING STATES AND THE 
SECOND ORDER PERTURBATION EQUATIONS 

The above results establish that for mixed incoming 
states, the leading order contribution to the smeared 
ANEC integral is always positive. However, for pure in- 
coming states the ANEC integral vanishes to 0{e), and 
therefore we need to investigate the second order per- 
turbation equation (4.23). In this section we calculate 
the O(e^) contribution (5.19) to the transversely smeared 
ANEC integral for solutions to Eq. (4.23), and show that 
it is positive in the long wavelength limit. 

To calculate the second order contribution to the 
ANEC integral, we shall need the second order contri- 
bution (5.12) to the Einstein tensor. Defining G^^ {x) = 
a^G^^^ix/a) we find from Eq. (4.20) that 



+0[{lnaf/a% (7.1) 
Here we have used the definitions 

h^^\x;a) = U^\x/a;a) = a'"^ x^^\x / a; a) (7.2) 

for j = 1,2, and have used the fact that x^^^ can be 
replaced with its leading order approximation x^^'*^-* to 
adequate accuracy, c.f., Eq. (3.52) above. The leading 
order metric perturbations x^^'*^^ and x^^'°^ are given by 
Eqs. (3.58) and (4.24). Finally, we insert the expansion 
(7.1) into Eq. (5.19) to obtain 

7(2) = J_ / 4°) 0f(-,— .-)A"A'' 
* na^ Jm ^ ^ a a a ' 

x{(Tr,c.L^Hi^Z«[x(-)] 



a" 



a." 



This equation when inserted into Eq. (5.10) will produce 
the expansion (6.42). Therefore it can be seen that the 
first three terms in Eq. (7.3) arc exactly analogous to 
those obtained from the first order analysis, and that 
consequently Eq. (7.3) can be rewritten as 



K/f)(AT)= lim lim Kjf)(A, Aj,, At) 

A— >oo Ai— >0 

A f\uoL (7 + AC ^ 

= ^ + + E — + O 

ct" or or 



+ ^x(^'°)^(Tr,4'H0[(ln«)V«l}. (7.3) 

We now exploit the close similarity between the first 
and second order perturbation equations. For the 
rescaled incoming state (3.41), the first order, reduced 
order equation (4.20) can be alternatively written in a 
form more closely parallel with Eq. (4.23): 

«G»'[/.'"i=^(rr.4;')+!^2ii'[x"'«i 

+^rt;'[x"'»>l,-,..«> + 0((!^)^ (7.4) 



(lna)2 



a' 



, (7.5) 



in analogy with Eq. (6.42). Here the coefficients A, B 
and C are functions only of the state perturbation Gj, 



(2) 



(2) 



and have the exact same functional dependence on 
as the coefficients A, B, and C in Eq. (6.42) have on 
Furthermore, the relevant positivity conditions on 



,-,(1) 

'^in' ^1^*^ same as corresponding conditions on 

Wjj^^; see subsection VII B below. The coefficient AC in 
Eq. (7.5) is defined to the contribution from the last two 

terms in Eq. (7.3), and depends both on cjj^^' and uj^^^ as 
well as the freely specifiable incoming piece of the metric 

perturbation xih'°''- 

Our strategy for proving the positivity of the second 
order ANEC integral is the following. First, in subsection 
VITA we show that AC = whenever A = B = 0. 
In subsection VII B, we show that the space of allowed 
second order state perturbations is effectively the same as 
that of the first order perturbations, and then to appeal 
to the first order analysis. 



A. Vanishing of the additional terms 

We now show that whenever A = B = 0, the last two 
terms in Eq. (7.3) vanish, and that consequently AC = 0. 
The arguments in this subsection will be mostly formal; 
we believe that these formal arguments could be trans- 
lated into rigorous arguments along the lines of the anal- 
ysis given in Sec. VI E 1 above, and using the positivity 
condition. However, we have not attempted to do so. 

It is not difficult to see that the last term in Eq. (7.3) 
vanishes. Since to the appropriate order in 1/a^ the 
transverse smearing is unimportant for this term, after 
the limit A ^ oo is taken it can be written as 



1 

2^ 



(1) 



(7.6) 



However, our choice of gauge guarantees that xiV = 
along 7. In particular this will be true to each order 
in the expansion of Eq. (3.52), so that x^aif^ also 

vanishes on 7, and therefore this term vanishes. 

The more interesting term in Eq. (7.3) is the second 
to last term, which can be shown to be proportional (at 
the relevant order in l/o?) to 
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' J '7 



7' 



where we have used Eq. (3.15). The meaning of the no- 
tation in the first term is that the function F^^^ is acted 
on by the first order change in the operator Vab induced 
by the metric perturbation x'^''^'^^ S'lid similarly for the 
second term. Now since A-;,'",],'"'' satisfies Eq. (3.58), it has a 
contribution both from the perturbation to the incoming 
state and from incoming classical gravitational radiation. 
Thus, Eq. (7.7) contains both a contribution quadratic 
in the incoming state perturbation uj\l^ and a cross term 

between uj^^ and the incoming classical gravitational ra- 
diation. 

Consider first the second term in Eq. (7.7). As men- 
tioned above we are assuming that the coefficients A and 
B vanish, since from Eq. (7.5) this is the only case in 
which the terms (7.7) are relevant. Therefore 



A = 



a;(2))A«A* =0. 



(7.8) 



The key idea that we now use is that this condition im- 

- (2) 

poses constraints on , which in turn imposes con- 
straints on u)l^^ sufficient to ensure that the second term 
in Eq. (7.7) vanishes. 

We now restrict attention to the case of Fock space 
states. Since the state is pure to first order, the first 
order perturbation to the density matrix is given by 



;3(^) = io)(vr^^+w|v>)(oi. 



(7.9) 



for some \tp)'"^^ S H with {0\tp)''^^ = 0. The most general 
second order density matrix perturbation is then 



= (1) (^| + |o)(v,|(2) 
+g-(l+trQ)|0)(0|, 



V)(0| 



(7.10) 



^2) (2.) ^ 

for some ' € H with (01-0) = 0, where Q is a 
positive, Hermitian trace class operator such that Q\0) = 
0. Since Q is trace class there will exist an orthonormal 
basis IV'j) of the space of states orthogonal to |0) such 
that 



into Eq. (7.8), then the contribution from the two terms 
in Eq. (7.10) involving \ ip) will vanish for the reason ex- 
i,ined in Sec. VI D. Therefore, using the formula (6.36) 
'and following the argument given in Sec. VIE 1, we find 
that 




0. (7.12) 



From the arguments given in Sec. VIE 1, it is clear that 
all the terms on the left hand side are individually non- 
negative, and hence they all vanish. It follows that 



(1) _ 



0, 



(7.13) 



where primes denote derivatives with respect to A and 
the notation means that the operator is evaluated on the 

geodesic 7. Next, using Eq. (7.9), we can rewrite the 
second term in Eq. (7.7) in an alternative notation as 

dA A"A^^O|V(„$(°V(,)$(^)|^^ , (7.14) 

Here <I>(^' is the first order change to the field operator 
induced by the metric perturbation Xa/^^ ■ Next, split 
up the zeroth order operator ^^^^ into its positive and 
negative frequency parts. The negative frequency parts 
will annihilate the vacuum on the left, and the positive 
frequency part can be commuted through the term 
(since the states |0) and IV')^^'' arc orthogonal), giving a 
result which vanishes by the condition (7.13). 

Consider now the first term in Eq. (7.7). To calculate 
this we need the explicit form of the operator Dab that 
enters into the point splitting prescription (2.19) for cal- 
culating the stress tensor. The general expression for Vab 
is given by [c.f., Eq. (2.4) above] 



V„bF^^\x,x') = Sab 

-aVaVbS - gabaS), 



-j^9ab9'"^Scd + ^GabS 



where 



S{x) = lim F.j^^\x,x') 



(7.15) 



(7.16) 



(7.11) 



for some qj > 0, j = 0, 1, 2 The overall density ma- 
trix will be positive to O(e^) whenever the perpendicular 

projections of IV")^^^ and \ip)^'^^ into the kernel of Q are 
orthogonal to each other. 

We now give a non rigorous argument for the vanishing 
of the second term in Eq. (7.7). If we insert Eq. (7.10) 



Sabix) = lim Va'VbF^^\x,x') 

= - lim Va^bF^^\x,x')+VaybS/2, 



(7.17) 



and Gab is the Einstein tensor. [Note that we are using 
a nonstandard notation in which the coincidence limit 
is implicitly understood in the symbol T^ab] Now the 
operator 'E'ab'fx^^'^^] in the expansion (3.1) will contain 



37 



pieces linear in Xab ^^'^ pieces that are linear in the 

derivative of x'i^b^^ ■ The pieces that are linear in xib'°^ 
will give a vanishing contribution to the first term in 
Eq. (7.7) since xii,^^ vanishes along 7 in our choice of 
gauge. The remaining piece of v'^^ yields 

+e /^^AG(V[X^''°)]5A«A^ (7.18) 

where C^J,'''^ is given by Eq. (5.6) with h^^jj replaced by 

xif,'"''- The first term in Eq. (7.18) can be seen to vanish 

using Eq. (5.6) and the fact that h^^fj = on 7 and that 

h^^^X'^X'' vanishes identically in our choice of gauge. Us- 
ing Eqs. (3.58), (6.36) and (7.9), we find that the second 
term in Eq. (7.18) contains the factor 

G'iV[x^''°^]A"A'' = (2/Ac)Re (o|$V(A,0,0)2|v)^'^ , 

(7.19) 

which vanishes by Eq. (7.13). Therefore the expression 
(7.7) should vanish. 

B. Smeeired positivity result 

We now explain how to adapt the perturbative smeared 

positivity result of Sec. VI to the present situation. The 
coefficients A, B and C in Sec. VI were expressed in terms 
of the function If{^t) (the flat spacetime ANEC integral 
along transversely displaced geodesies in the state ^^\^)^ 
and the only properties of this function necessary to prove 
the result was that it was a pointwisc positive smooth 
fimction. Correspondingly the coefficients A, B and C 
can be expressed in terms of the analogous function 

lf{^T)^ f (Ti°),a;L^))A<'A^ (7.20) 

This quantity can be expressed in terms of the mixed fre- 
quency part g^'^^ of the second order two point function 
by equations analogous to Eqs. (6.30) and (6.34). How- 
ever, from Eq. (6.11) we find that g^"^^ obeys a stronger 
positivity condition than g^^^ in the case we are consider- 
ing when g^^'i = 0. In particular, g^"^^ satisfies an analog 
of Eq. (6.9) and it follows that the function (7.20) must 
be nonnegative. 

Hence the coefficients A, B and C obey the same pos- 
itivity conditions as the coefficients A, B and C: A>Q, 
B > whenever A = 0, and C > whenever A = B = 0. 
Moreover, the case A — 13 = C — can be excluded 
in the following way. If all these coefficients vanish, an 
argument similar to that given in Sec. VIE 2 shows that 



the operator Q in Eq. (7.10) must vanish, and that con- 
sequently IV') = p^^^ = 0. Therefore, by defining e' = 
we see that we are really dealing with a first order per- 
turbation instead of a second order perturbation. 
^ Since we showed above that AC = whenever A = 
B = 0, we conclude from Eq. (7.5) that for sufficiently 
large a, for pure states the transversely smeared ANEC 
integral is positive. 

C. Second order vacuum polarization 

As we have explained, the second order vacuum po- 
larization term (T^^^ [/i'-^^ /i'^'], Win^o) is not explicitly 
known, but arises at sufficiently high order in our long- 
wavelength expansion that it can be neglected, when we 
assume that incoming gravitational radiation does not 
dominate the first order metric perturbation. However, 
if we drop the assumption on incoming gravitational radi- 
ation, and assume in addition an incoming vacuum state, 
then we can derive a condition that this unknown vac- 
uum term must satisfy in order for ANEC not to be vi- 
olated. This condition provides an additional test of the 
ANEC hypothesis which is independent of our analysis 
above. Moreover, if the condition is satisfied, then it can 
be shown that without any assumptions restricting the 
incoming classical gravitational radiation, that the trans- 
versely smeared ANEC integral is always nonnegative for 
solutions of the reduced order equations, for general, non- 
vacuum incoming states, in the long wavelength limit. 

The condition we find, by carrying out a reduction of 
order to the appropriate order of the perturbative semi- 
classical equations, is the following. Let h^^^j be any solu- 
tion of G''^fj[hM'>] = 0, and let h'f^ be any solution of 
Gablh^^^] + Gl^^\h<-^\ h^^^ = 0. Thus, the spacetime 

{M,rjab + sh^ab + ^"^^lib) Satisfies the vacuum Einstein 
equation to second order and consists of classical gravi- 
tational waves. Then the quantity 

(ri^)[ft(i),ft(i)],u;i„,o) + {Til\h^%u;,^,o) (7.21) 

describes the expected in-vacuum stress-energy tensor of 
the quantum field to second order on this spacetime. 
Moreover, the stress tensor (7.21) does not depend on 
which solution of cl^ijih^^^ + G''^^\h^^\ h^'^'>] = is 
chosen, and thus is a function only of /i^^^ . The condition 
is that the ANEC integral of the quantity (7.21) should 
always be nonnegative. We conjecture that this is the 
case. 



VIII. CONCLUSIONS 

In this section we recap briefly our main assumptions 
and assess the significance of our results. We have ex- 
amined the positivity of (transversely smeared) ANEC 



38 



integrals for solutions of the reduced order semiclassi- 
cal Einstein equation. Three small parameters have ap- 
peared in our analysis: e, measuring the deviation of the 
metric from the flat metric and of the quantum state 
from the incoming vacuum; h or equivalently L'^j C? ^ our 
"long wavelength" expansion parameter, and finally l/A'^, 
where A'' is the number of scalar fields coupled to gravity. 
We have calculated the ANEC integral to leading order 
in to quadratic order in £, and to the first three 

nonvanishing terms in C? . We restricted our analysis 
to the case where incoming gravitational radiation does 
not dominate the first order metric perturbation. Apart 
from this restriction, we have shown that the transversely 
smeared ANEC integral for nearly flat spacetimes is al- 
ways strictly positive along every null geodesic in the long 
wavelength limit, except in the trivial case of the vacuum 
state in Minkowski spacetime, where the ANEC integral 
vanishes identically. 

There were several independent places in our analysis 
in which, a priori, a violation of ANEC could easily have 
arisen. In particular, the ANEC integral for pure states 
need not have vanished at first order in e; the coefficient 
(6.51) need not have been of a defenite sign; or any one 
of the extra terms that appear at second order in e in the 
coefficient AC need not have been identically vanishing. 
Indeed, several times during the course of this work we 
believed to have discovered a serious violation of ANEC, 
only to find on more careful analysis that this was not the 
case. Therefore, we consider our results to be evidence 
in favor of the conjecture that ANEC comes sufficiently 
close to holding in general solutions of the semiclassical 
equations to rule out macroscopic traversable wormholes. 
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APPENDIX A: EXACT SOLUTIONS OF FIRST 
ORDER PERTURBATION EQUATIONS 

In this appendix we derive all of the exact solutions to 

the first order semiclassical equation (3.38) whose spatial 
Fourier transforms exist. Some of these exact solutions 
have been discussed by Horowitz [48] , in the special case 
of the homogeneous version of the equation, without the 
source term (3.43). Here we generalize the treatment 
of Horowitz to allow for first order perturbations to the 



quantum state. To solve Eq. (3.38) we can fix the in- 
coming state perturbation ij'^ (this is freely specifiable 
up to some regularity conditions discussed in Appendix 
C), and solve for the metric perturbation l'^'^^ . Because 
the equation depends on h^^^ only through its linearized 
Einstein tensor G^^^ , we can first solve for G^^^ and then 
use this to obtain h^^^ , as suggested by Horowitz [48]. 

The exact solutions to Eq. (3.38) are closely analogous 
to the solutions of the Klein-Gordon equation with neg- 
ative mass squared 



(□ - rr?)^{x) = p{x), 



(Al) 



where p{x) = p(x, t) is a source. We start by recalling the 

nature of the solutions of Eq. (Al). The general solution 
can be written in terms of the spatial Fourier transform 
p(k, t) of the source as $ = $< -|- where 



and 



.(x,i)= / 



|k|>m (27r)3 

/OO 
dt' G,. 
-co 



ikx 



A(k)e^'^''* + S(k)e-"^''* 



,{t - i';wk)p(k,t') 



(A2) 



(x,i)= / 



d^k 



+ 



(27r)3 

/OO 
dt' G.. 
-OO 



|k|<m 
cso 



C(k)e'"''* -I- Z)(k)e-'"''* 



sym,exp 



{t-t';Ki,)p{k,t') 



■ (A3) 



Here Wk = Vk^ — m^, Kk = Vrr? 



k2, 



sym,osc(^5 ^) — 

sin(w|t|)/(2w), Gsym,exp(t; k) = -e-''l*l/(2K)„ and the 
functions A(k), . . . , D(k) are arbitrary except for the re- 
ality conditions A{—'k) = A(k)* etc. Following Horowitz, 
we will refer to the portion (A2) of the solution as the 
tachyonlike or oscillatory part, and the portion (A3) as 
the exponential part. 

We explicitly display these solutions to the negative 
mass squared Klein Gordon equation because the solu- 
tions to the semiclassical equation (3.38) have a very 
similar character. In particular these solutions can be 
divided into "oscillatory" and "exponential" pieces. We 
will obtain the general solutions by spacetime Fourier 
methods. As background, we start by recalling how to 
obtain the solution (A2) - (A3) of Eq. (Al) by spacetime 
Fourier transforms (as opposed to merely spatial Fourier 
transforms). 

It is clear that the "tachyonlike" portion $> of the 
general solution can be straightforwardly obtained using 
Fourier transforms, when the source p(x, t) is sufficiently 
well behaved. The Fourier transform with respect to time 
of the Greens function obtained from Eq. (Al) has two 
poles on the real iv axis, and the choice of "ie" regular- 
ization prescription is equivalent to the choice of Greens 
function; for instance, as is well known, demanding that 
G{lu) be analytic in the upper half plane yields the re- 
tarded Greens function. The freely specifiable first term 
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in Eq. (A2) can clearly be written down by inspection 
using the location of the poles of the Greens function. 

The situation is slightly different for the exponentially 
growing portion (A3) of the solution. In this case the 
Greens function in the Fourier domain 



(A4) 



has no poles on the real axis, and hence Fourier trans- 
form methods produce a unique Greens function in the 
time domain, which is just the particular Greens function 



G. 



sym,exp 



{t) that we choose to use in Eq. (A3). Of course. 



other Greens functions appropriate to different boundary 



conditions do exist in the time domain, but Gg 



is the only one whose Fourier transform exists. Hence, 
solving Eq. (Al) using Fourier transform methods will re- 
produce the last term in Eq. (A3). The freely specifiable 
first two terms in Eq. (A3) are not directly obtained, but 
clearly again can be written down by inspection, using 
the location of the poles in G{oj). 

We now turn to a similar analysis of the semiclassical 
equation Eq. (3.38). We can obtain very general solutions 
to the equation by determining the analytic structure 
in oj of the appropriate Greens function in the Fourier 
domain. See Sec. VI B and Appendix C for a discussion 
of the existence of the Fourier transform of the source 
term (3.43). Using Eqs. (3.38), (3.33) and (6.17) we find 
the following formal expression for the general Fourier 
transformable solution 



l^^ab W- 



+ 



F2{k) 



{kakb Tjabk , 



Fi{k)Fs{k) 



where 



Fi{k) = 1 - IGnaL^pk^Hxik) 



F^{k) = (87r) 



{-a + 2b)LlHx{k)+2pLl 



and 



F3{k) = l + 6Llk''[/3 + bHx{k)]. 



(A5) 



(A6) 



(A7) 



(A8) 



While Eq. (A5) is not the general solution we are 

looking for, it is straightforward to write down the gen- 
eral solution by inspection, essentially by adding extra 
terms to Eq. (A5) that correspond to poles of the func- 
tions 1/Fi and l/(i^iF3). Consider first the function 
1/Fi. From Eqs. (6.17) and (A6) we can write Fi as the 
limit of a function analytic on the upper half u) plane: 
Fi{k) = lim£_>o+ Gi(k, (J -|- ie), where 



C?i(k,w) = 1 + 



In 



A2(a;2_^2 



(A9) 



= |k|, Lol = l/{32TT^aLl), A = Aexp(7 - 1/2), and 7 
is Euler's constant. The function Gi has branch cuts on 
the real axis at w > and lu < — Wk- The poles in the 
Greens function are just the zeros of Gi . These location 
of these zeros depend on A and Wk in the following way 
(see Fig. 2). Define 



Acrit = 4:TrLr,e '^V2a, 



(AlO) 



which is a fixed length of the order of the Planck length. 
Let Zi, 1 < i < 4, denote the four roots of the equation 

1 -|-zln[zA^cj,":] = 0. These complex roots depend only on 
the ratio A/ Acrit- Then the zeros of the Greens function 
(A9) are 



l<i<4. 



(All) 



When A > Adit, there are four complex zeros ±ll!,±lu*, 
where w = cto + if^o hcs in the first quadrant. When 
A < Acrit there are three separate cases. When 
exp(— cj^/a;^)/cj^ > A^, then there are two real 
roots ±ao and two imaginary roots ±i/3o- When 
exp(— a;^/a;^)/a;^ < A^, then if Wk > i^c there are four 
real roots ±ao, ict'o, and if Wk < Wc there are four imag- 
inary roots ±2/3o, aPo- 

The analytic structure of the function 1/F3(fc) is some- 
what simpler. There exists a unique A;o > depend- 
ing on A such that the locus of the zeros of F3{k) is 
k"^ = Lul - w'^ — k^. 

The general solution to Eq. (3.38) can be written as 

G^:^ = GiV' + gT"''' + GlV' (A12) 

We discuss these three terms in turn. The inhomoge- 
neous part of the solution is given by 



+ 



^ k ik-x 

(2^)4 

F2{k) 
Fi{k)Fs{k) 



Sab 



Fi{k) 

{kakb - Vabk'^)Sc 



. (A13) 



In this expression it is understood that the poles on the 
real lo axis in the functions 1 /Fi and I/F3 are regulated 
with the appropriate ie prescription to pick out the "half 
retarded plus half advanced" type contribution from each 
pole. [We choose this particular prescription for conve- 
nience, it would clearly be possible to use instead for ex- 
ample the "retarded" type contribution from each pole.] 
As discussed above, poles on the real oj axis will occur 
only for A < Acrit for the function 1/Fi, but will occur 
for all values of A for the function I/F3. Note also that 
in the conformally coupled case ^ = 1/6, the trace of the 
source tensor Sab vanishes, and hence the second term in 
Eq. (A13) vanishes. However, as we discuss below, the 
freely specifiable piece of the solution associated with the 
second term in Eq. (A13) [Eq. (A16) below] does not van- 
ish even for conformal coupling. 
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The second term G 



(1), free 



in Eq. (A12) is the freely 



specifiable, homogeneous piece of the general solution as- 
sociated with the Greens function 1/Fi{k). In the case 
A > Acrit5 it can be written as [92] 



Ca6(k)e- 



+£>„6(k)e+'«°*e 



-0ot 



+ C.C- 



(AM) 



Here C„5(k) and Dabik) are arbitrary except that they 
must be traceless and satisfy k°'Cab = l°'Dab = 0, where 
k°' = (k, ao+i/So) and = (k, — ao — «/3o)- The quantities 
ao and /3o depend on the mode frequency Uk = |k|, and 
on the Icngthscalc A as indicated by Eq. (All). This is 
purely an exponentially growing/ decaying type solution 
and has no oscillatory parts. 
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FIG. 2. An illustration of the dependence of the locations 

of the poles of the Greens function XjFx (k) on the parameters 
LOk, the frequency of the plane wave mode in question, and A, 
the additional free parameter with dimensions of length that 
appears in the quantum theory but not in the classical theory. 
In general there are four poles of the form ±«i, ±Z2, where 
the locations of the two poles zi and 22 arc as follows. In the 
hatched region below the curve, there is one real pole and one 
imaginary pole. Above the curve and below the dashed line, 
in the left region both poles are imaginary, and in the right 
region both poles are real. Finally, in the region above the 
dashed line, both poles arc complex. Acrit is a critical length 
of order the Planck length, and lOc is a critical frequency of 
order the Planck frequency. 



( 1) free 

The homogeneous term G^^ is more complicated 
in the case A < Acrit- In this case there are again freely 



specifiable terms for each of the poles of 1/i^i (sec Fig. 2). 
The poles off the real axis will give exponentially grow- 
ing/decaying terms, and the poles on the real axis will 
give "tachyonlike" contributions analogous to Eq. (A2). 
In particular, there will be a portion of the general so- 
lution which has the form of an integral over the values 
of k for which both poles are real. This expression will 
have the same form as of Eq. (A2), except that A 
and B arc replaced by transverse traceless tensors, and 
the quantity is replaced by either of two independent 
constants depending only on A (i.e., there is a sum of two 
terms). This portion of the general solution was obtained 
by Horowitz [48]. It can be written as 



tachyon 



{x)ab = i 

Jn 



Lab{k)e 



ik'X 



(A15) 



kl > 



Here TC is the union of two spacelike hyperbola k'^ 
and fc^ = fc| > in momentum space, where ki and k2 
are constants depending only on A. The tensor La,b{k) is 
a freely specifiable transverse traceless tensor on H which 

falls off sufficiently fast at infinity. 

The third term G^V'''*''''^ ''^^ ^q. (A12) is a purely 
transverse, freely specifiable piece of the solution which 
is associated with the Greens function I/F3. It can be 
written as [92] 



|k|<fe, 
+ 
-l-c.c 



/ d^k e*^ 



^(k)e'*''* + F(k)e-«''* 
/(k)e*''''* + J(k)e-*''''* 



{kakb - Vabk'^) 



{kakb - Vabk'^) 
(A16) 



Here = y/k'^ — k^, Kk = V^'o ~ the func- 

tions E, F, I and J are freely specifiable functions of 
k. Thus, there are both tachyonlike and exponentially 
growing modes of this type for all values of A. Also, as 
remarked above, this freely specifiable transverse piece 
of the solution does not vanish in the conformally cou- 
pled case ^ = 1/6, despite; the fact that the; the analo- 
gous transverse contribution to the inhomogeneous piece 
of the general solution [the second term in Eq. (A13) 
above] docs vanish for conformal coupling. 

Finally, we remark that exact solutions to the alterna- 
tive, rescaled version (3.45) of the first order perturbation 
equation can be obtained from the above analysis using 
the substitutions h^^^ij '^^^'^b' ~^ Lp/a, A A/c 

\-^ab 



and 



•Sab 



Sab 



APPENDIX B: ANEC INTEGRAL FOR EXACT 
SOLUTIONS 

In this appendix we consider the specific subclass S 

of the solutions of the exact, first order, scmiclassical 
equation discussed in Appendix A, given by using half 
advanced plus half retarded Greens function to obtain 
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the linearized Einstein tensor. In the case of exponen- 
tial type modes discussed in Sec. IV C above, the use 
of this Greens function to pick out a class of solutions 
is equivalent to throwing away the runaway solutions 
by hand. For the oscillatory type modes, this choice of 
Greens function yields a particular subclass of solutions. 
We shall show that the solutions in S have the property 
that their transversely smeared ANEC integral is always 
non-negative, even outside of the long wavelength limit, 
whenever A > Acrit- However, we also shall show that 
some exact solutions outside of this subclass do violate 
ANEC. 

As explained in Sec. IV C above, we see no reason to 
view solutions in S as being any more physically mean- 
ingful than any other subclass of exact solutions. How- 
ever, given a solution in <S, any other exact solution ob- 
tained from the same incoming state will have the same 
perturbative expansion in and In a/a^ [or equiva- 

lently, in h and h\n{h)]. This perturbative series up to 
any finite order also should coincide with what would bo 
obtained by carrying the reduction of order procedure 
of the semiclassical equation to the appropriate order in 
l/a^ and solving exactly the new, reduced order equa- 
tion. Thus, an expansion in l/a^ and Ina/a^ of the 
positivity result of this appendix provides an alternative 
proof of the results we established in Sec. VI above for the 
solutions of the reduced order equation (4.20), at least for 
A > Acrit- In other words, we can use the analysis of the 
solutions in <S as a mathematical tool to to establish a 
positivity result for solutions of the reduced order equa- 
tions. The alternative proof that this appendix provides 
also gives insight into the otherwise mysterious positiv- 
ity properties of the coefficients A, B and C (with and 
without smearing) discussed in Sec. VI E 2 above. 

We now turn to a proof of the positivity of the ANEC 
integral for this class of solutions. The analysis of the 
ANEC integral for the exact solutions to the unmodified 
semiclassical equation parallels that given in Sec. VI B 
for the solutions of the reduced order equation, with the 
only difference being that the functions 5*1 and ^2 of 
Eqs. (6.19) and (6.20) are replaced by the expressions 
and F2 / F3 respectively, where the functions Fi , F2 
and F^ were defined in Appendix A. Correspondingly we 
again obtain the formula (6.35), except that the function 
K{xt) is replaced by Ki{xt), where 



i^l(AkT) = 



Ak; 



In (a^ Ak^^ 



, -1 



(Bl) 



Note that if A > Acrit then KiCkT) is a continuous 

function; its Fourier transform Ki{'Xt) is and is con- 
tinuous everywhere away from the origin. The function 
Ki{xt) for A — 1.2 Acrit is shown in Fig. 3. For A < Acrit, 
the expression (Bl) blows up at some finite value of jkyl 
so that Ki is not even L^. In this case an appropri- 
ate ie rcgularization prescription should be understood 
to apply to the formula (6.35), the precise prescription 



being determined by the fact that we are choosing half re- 
tarded plus half advanced solutions. The rcgularization 
prescription yields a well defined distribution Ki, with 
well defined distributional Fourier transform Ki (xy) . We 
assume from now on that A > Acrit; however, it is possi- 
ble that the results in this appendix continue to hold for 

A < Acrit- 



1000 




x^l (in Planck lengths) 



FIG. 3. The ANEC integral for exact, "half advanced plus 
half retarded" solutions to the semiclassical equation can be 
expressed in terms of the function 7f(xt) [which is obtained 
from the ANEC integral in Minkowski spacctimo evaluated 
on transversely displaced geodesies] integrated against a par- 
ticular function ii'i(xT) with width of order of the Planck 
length, see Eq. (6.39) above. Here wc plot the function Ki 
as a function of the transverse distance |xt| from the original 
geodesic, in the case where A = 1.2Acrit. The fact that this 
function is negative for some values of its argument implies 
that there are incoming states for which the ANEC integral 
including backreaction is negative. 



From Eq. (6.35), the smeared ANEC integral is pro- 
portional to 



/ 



/F(xT)[i^lo5dil](xT), 



(B2) 



where S'dii(xT) = S{xt/At) is a "dilated" smearing func- 
tion and Ki o S^a denotes the convolution of Ki and 5dii. 
Since we know If{x.t) is positive, to make the integral 
(B2) positive it suffices to find smearing functions whose 
convolutions with Ki are pointwise positive everywhere. 
We now show that for the choice (6.40) of smearing fimc- 
tion, there exists some At > such that the function 
5(xt/At) will satisfy the positive convolution condition. 
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In fact, our proof below can be easily extended to apply 
to any positive smearing function which does not fall off 
at large xt more rapidly than Ki does (oc [xtI^^); in 
other words, any function which falls off slowly enough 
will be a suitable smearing function when sufficiently "di- 
lated". 

We now outline a proof of that the convolution is point- 
wise positive. We would like to show that there exists 
some A such that the function 



/i(A,xo) J rf^xJri(Ax)5(x + xo) 



(B3) 



is positive for all xq. [For the remainder of this appendix, 
we drop the subscript T on At-] Clearly, for any fixed xq, 
we can find a A such that the integral (B3) is positive, 
since A^i^i(Ax) <^^(x) as A ^ oo. Let Amin(xo) de- 
note the smallest number Aq such that the integral (B3) 
is positive for all A > Aq. Then we need just to show 
that A„iin(xo) is bounded above as a function of xq, so 
that there exists some positive A which works for all xq. 

To show this we suppose conversely that there exists 
some sequence ( Aj , xq, j ) such that Aj oo and 



/i(Aj,xo,j) = 0. 



(B4) 



Now suppose that the sequence x^j is bounded. Then 
by passing to a subsequence we can without loss of gen- 
erality assume that Xo,j — > xq for some xQ) so that 
{Aj^ ,xoj) (0,xo). However, the fimction /i(A,xo) 
is jointly continuous as a function of A~^ and Xq even at 
A~^ = 0, as can be deduced from the formula 



h= / dV^i(y)5(xo+A-V) 



(B5) 



and by using the properties of the functions S and 
Ki. Hence it follows that Ii{Aj,xoj) /i(oo,xo) = 
S'(xo) > 0. This contradicts Eq. (B4) above, and so we 
conclude that the sequence xq^j is unboimded. 

To exclude the possibility of the existence of such a se- 
quence with Xq, j unbounded, we now derive an estimate 
for the function 7i(A,Xo) for large Xq. First, we note 
that Ki{x) ~ — |x|~'' at large |x|, and in particular that 
there exists a fco > and eo > such that for all e < eo, 



l^^i(x)| <e 



(BG) 



whenever |x| > koe~^^^ [90]. Moreover, from the formula 
(6.40) for iS'(x) we can show that for any r] with < r] <1, 
there exists fci > such that 

(1 - r/)5(xo) < 5(x + xo) < (1 + v)Si^) (B7) 

whenever |x| < fci|xo|. Let be the integral of |i^i| 
over the domain where it is positive, and similarly define 

> 0, so that 



1 = j (fxKi{x) =K+-K- 



(B8) 



Consider now the integral (B3). For any e with < 
e < eo , we can split up this integral into three differ- 
ent parts: There is a contribution from the region where 
|x| > ko€~^^'^/A, which is bounded below by — CqcA^ 
from Eq. (B6), for some constant Cq. There is a contri- 
bution from the portion of the region |x| < koe~^^'^/A 
in which Ki is positive. Using Eq. (B7), this will be 
bounded below by - r])S{xo)/A'^ , ff 



fcoe"'/VA< fci|xo|. 



(B9) 



We can ensure that this condition holds by choosing e = 
(fci|xo|A/fco)~^; this requires that 



|xo| A > (A;o/fci)eo' 



(BIO) 



Finally, there will be a contribution from the correspond- 
ing region in which Ki is negative, which with the above 
choice of e is bounded below by —K~{1 + rj) S {xq) / A"^ . 
Combining these bounds and choosing r] = {2 + 4:K~)~^ 
yields that 



7i(A,xo) > 



,r4A2 



1 



2(1 + 4) 



(Bll) 



where Ci is a constant, whenever the condition (BIO) is 
satisfied. It clearly follows that the function Ai„i„(x()) 
is bounded above for large xq. This contradicts the un- 
boundedness of the above sequence xq, j, and we conclude 
that the function Amin(xo) is bounded above for all xq. 

We conclude that there is some fixed length A^ such 
that for the smearing function (6.40), the ANEC inte- 
gral (B2) is always non-negative. It is clear on dimen- 
sional grounds that the critical value of A^ is of order 
of the Planck length. We remark that this positivity re- 
sult would not hold for any smearing function which falls 
off more rapidly than \Ki\. Thus, for example, there 
is no positivity result for Gaussian transverse smearing. 
However, as mentioned above, our proof could easily be 
modified to apply to smearing functions that fall off more 
slowly than \Ki\. 

Finally, we consider the exact solutions of Eq. (3.38) 
outside of the subclass S of half retarded plus half ad- 
vanced solutions. We now show that when A < Acrit, 
there exist exact solutions outside of S which violate 
ANEC, even when transversely smeared. [It can be 
shown that violations of ANEC do not occur when A > 
Acrit , if we we discard the exponentially growing and de- 
caying pieces of the solutions.] Using the momentum 
space coordinates /3, 7, ky, the first order ANEC inte- 
gral (5.10) can be written as 



7(1) 



rf'kTG'«(/3^ 



0,7,kT)A»A''. (B12) 



For the tachyon type solution (A15), the linearized Ein- 
stein tensor can be written as 



GabiP, 7, kr) = Si-2jp + - k^) J^abiP, 7, ^t), 



(B13) 
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where Tab is transverse and traceless but otherwise freely 
speeifiable and = kr/jkrl- From Eqs. (B12) and 
(B13), and writing kr = (cos x, sin x), it follows that 



/oo /'27r -1 

dl dx;^-?-a6(0,7,x)A''A\ (B14) 

Since J-"ob(/?, 7, k^) is freely specifiable away from 2/?7 + 
Uq = 0, it is clear that in general I^^^ is non vanishing 
and can be of cither sign. Moreover, it is easy to see 
that a negative I'-^-' cannot always be made positive by 
transverse smearing in this case. 

Note that if h[^^{x; a) is a one parameter family of 
exact solutions such that ANEC is violated, then it will 
still be true that ANEC will be satisfied order by order 
in and Ina/a^. In other words, the violations of 

ANEC are non-perturbative in or equivalently in 

n. 



APPENDIX C: EXISTENCE OF SPACETIME 
FOURIER TRANSFORM OF EXPECTED 
STRESS-ENERGY TENSOR 

In this paper we have imposed a condition on the 
Minkowski spacetime in-states we consider, which re- 
quires that the two-point functions be smooth and have 
suitable falloff properties at spatial infinity. In this ap- 
pendix we derive some implications of these falloff as- 
sumptions, which are used in the body of the paper. 

We start by describing the class of states we are con- 
sidering. Let the two point bidistribution of the state 
be 



G{x, x') = Gq{x, x') + F{x, x'), 



(CI) 



where Go{x,x') is the vacuum two-point bidistribution 
in Minkowski spacetime (which was denoted Gin.o in the 
body of the paper). Throughout this appendix, we drop 
the subscripts "in" that appeared in the body of the pa- 
per. Also the conditions on F that we discuss below 
apply to both the functions F-j-^^ and F.j^^ that appear 
in the body of the paper, and our conclusions about the 
expected stress tensor in the state (CI) clearly then will 
apply to the source term (3.43) and to the first term in 
Eq. (3.54). 

The function F is determined by the restrictions of the 
four fimctions F, dF/dt, and dF/dt' and d'^F/{dtdf/) to 
S X E, where S is the hypersurface t = 0. We assume that 
all these functions lie in the class V defined in Sec. VI — 
that is, that they are smooth and that all of their spatial 
derivatives are on S x S. Note that imposing this 
condition on any other surface of constant time would 
yield the same class of states. [It would be sufficient 
for all of our results to assume only that the all spatial 
derivatives up to seventeenth order are L^; we have not 
investigated what the sharp falloff requirements are.] 



As in the body of the paper, we can express the two 
point function in terms of functions / and g via 



+C.C. 



(C2) 



where k = {k,uik), u)k = |k|, and "c.c." means the 
complex conjugate. This equation defines / and g to 
be suitable complex linear combinations of the spatial 
Fourier transforms of the above-mentioned four functions 
restricted to S x S: we have 

11 1 1 ] / ^(k,k') \ 

11-1-1 iFi{k.k')/uji, 

1-1-1 1 iF^t'{k,k')/u>w' 

-iJ \-F,tt'{k,k')/{uuwi,>) J 

(C3) 




1 -1 1 



Our assumptions on F imply that / and g are continuous 
away from k = and k' = and satisfy, for any integer 
N, 



ma^{|/(k,k')|,|5(k,k')|} < 7^^^^ 



(l-hk2+k'2)^' 

1 \ / 1 



(C4) 



for some constant Cn- [Note that the requirement that 
the total energy of the state be finite should imply an 
inequality of the form (C4) but without the 1 -|- 1/uiu 
factors and with N = 5/2 + e.] Moreover our assumptions 
imply the functions / = Wk^k'/ and g = Wk'^k'S defined 
before Eq. (6.15) are at least C^, a fact which is used in 
the body of the paper. 

Now the expected stress tensor in the state (CI) will 
be automatically smooth, from the relation 



Tab{x) ^VabF{x,x') 



(C5) 



and the fact that the initial data for F on S x E is smooth, 
so that F itself is smooth on M x M. Therefore its 
Fourier transform Tab{k) exists as a distribution. We now 
investigate the regularity properties of this distribution. 

In our notation below we make use of delta functions; the 
steps can be made rigorous by integrating all equations 
on both sides against a smooth test tensor field /"^(a;) of 
compact support on M . The expression for 2?ah is given 
in Eq. (7.15) above. We restrict attention below to the 
piece 

Sab{la) = [ d^xe-^'-^" lim VaVb'F{x,x') (C6) 

J x'^x 

of this expression; a similar analysis applies to the other 
type of term involving VaVfeS' in Eq. (7.15). From 
Eqs. (CI) and (C2) we obtain that 
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Sab{la) = Fabila) + FabHaT , (C7) 

where la = (1, w) and 

Fab{l,Lo) = yd^kwkWk-i (l,n)a 

X [6{u! - Wk + Wk-i)5(k, k - 1)(1, -n')6 
-S{uj - Wk - ^k-i)/(k, 1 - k)(l, n'),,] . (C8) 

Here n and n' are unit vectors in the directions of k and 
k — 1 respectively. We write the integral over k as 

I dkk^ dii d(pk, (C9) 
Jo J-i Jo 



where /i is the cosine of the angle between k and 1. The 
(5-function in the first term in Eq. (C8) can be rewritten 
as 



e[k-{w + l)/2] e{l-\co\)S{ii- iii{k,l,uj)) 



\k — ui\ 



kl 



(CIO) 



and similarly the ^-function in the second term becomes 

\k-Lo\ 



Q[l/2-\k-uj/2\] Q{uj - l)5{n- iJL2{k,l,uj))- 



kl 



(Cll) 

Here 8 is the step function, and jii and 1x2 are the ap- 
propriate values of /i that are enforced by the ^-function. 
Equations (C8)-(C11) yield 

Fab{\, W) OC ^^^"l'^'^ r dk ( difk k\k - Ujf 

X [(l,n)„(l,n')6 g(k,k-l)](/xi) 
^^"^-^l t^'^'\k I d<fikk^{k-wf 



I 



X [(1, n)„(l, -n')6 /(k, 1 - k)] (/.2). (C12) 

From the properties of the functions / and g and the form 
of Eq. (C12), it is clear that Fat and hence also 5*06(1, oj) is 
continuous everywhere away from the lightcone, and that 
|l|i^af,(l, uj) is bomidcd in a neighborhood of I = 0. Similar 
conclusions apply to the entire stress tensor Tab{k). 

Finally, we note that the stress tensor Tab{k) is L^, 
from which it follows that Tab{:f) is on Minkowski 
spacetime. We can prove this using the formula (C12) 
as follows. Let e"*" be any constant tensor, and let 
denote the first term in Eq. (C12). Then from Eq. (C4) 
and using the fact that the 5 functions in Eq. (C8) enforce 
|k — 1| = Iw — it follows that 



j \e"''F^,{l,u;)fdHdiJ <Co £ dl j' 



duj 



/ dk 

l(l+u^)/2 [1 + fc2 + (fc - wY] 



N 



(C13) 



where Co is a constant, which is finite for N > 5. A simi- 
lar analysis applies to the second term in Eq. (C12), and 
to the other types of term involving VaVftS in Eq. (7.15). 



APPENDIX D: FOURIER TRANSFORM OF THE 
FUNCTION KT^log[K|] 

In this appendix we prove the following result which is 
used in the body of the paper. 

Lemma Let h(x) be a non-negative, smooth, 
function on such that 

h{0) = hi{0) = hij{0) = hijkiO) = 0, (Dl) 

where the commas denote partial derivatives. Then 



I HI 



k2log(k2) >0, 



(D2) 



with equality iff /i = 0. 

Proof: The essential idea is that the Fourier trans- 
form of k^ log(k2) consists of a smooth positive function 
away from the origin, plus distributional contributions at 
the origin whose effects are unimportant because of the 
condition (Dl). 

Using polar coordinates k,ip,x,x defined by x = 
(a;cosx, a;sinx), k = {k cos{x + ^) , k sm{x + f)) , we can 
write the integral (D2) as 

/•2'ir POO 

/ dxlim / dxh{x,x)I{K,x), (D3) 

where 

j-K i-2v 

I{K,x)= dk dipk^\og{k^)e''""', (D4) 
Jo Jo 

and n = cos (fi. By evaluating the integral over k in 
Eq. (D4) and using the identities 



and 



dip 



(2m2 + 1)(1-M^) 



it can be shown that 

I{K,x) = 12tt/x^ 



d(p AJ{K, cc, n). 



(D5) 



(D6) 



(D7) 



Here A1{K, x, fj,) consists of linear combinations of terms 
of the form 



{log K)iKPiJ,'' 



cos{Kij.x) 



(D8) 



or 
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{logKYKPii" 



sva.{Kfjbx), 



(D9) 



where q,p,w and n are integers with 0<9<1, 0<p<5, 
< w < 8 and < n < 3. 

After integrating over x, the contribution from such 
terms becomes proportional to 



~g^,n{K^,){\ogKYKP^i- 



(DIO) 



where g-^^n{x) = h{x,x)/x'^ is the restriction of the 
smooth function /i(x)/x" to the line at angle x the 
X plane, and g^^n is the (one dimensional) Cosine or Sine 
transform of this function. The function g^j^-^n is smooth 
because of the condition (Dl) since < n < 3. It fol- 
lows that g^ji{K) will fall off at large K faster than any 
power of K, and hence the absolute value of the expres- 
sion (DIO) is bounded above by 



\\og K\'iKPiJL^Cp 
1 + {txK)P+^ ' 



(Dll) 



where we have chosen the power to be p -|- 1 and Cp is 

a constant. Carrying out the integral over Lp yields an 
expression which is bounded above at large K by 



\\ogK\iKPC' 
C" + KP+'^ ' 



(D12) 



where C and C" arc additional constants that depend 
on p and w. Now it can be seen that these terms give a 
vanishing contribution in the limit if ^ oo in Eq. (D3). 
Thus, from Eq. (D7) the original integral (D2) can be 
written as 



dx / dxh{x,x) 
Jo 



which is manifestly positive. 



12tt 



(D13) 



APPENDIX E: CONSISTENCY OF OUR 
RESULTS WITH EXAMPLES OF NEGATIVE 
(UNSMEARED) ANEC INTEGRALS IN 
SELFCONSISTENT SOLUTIONS 

As mentioned in the Introduction, in a very recent pa- 
per Visscr [43] has shown that the expected stress ten- 
sor of a scalar field in the Boulware vacuum outside a 
Schwarzschild black hole violates ANEC. Visser argued 
that this result suggests that similar violations would oc- 
cur in self-consistent solutions with backreaction. More- 
over, these violations of ANEC would not bo confined to a 
Planck-scale tube surrounding a particular null geodesic, 
but would instead occur over a macroscopic region. In 
this appendix we review Visscrs argument and show that 
— at least in the context of perturbation theory off of 
Minkowski spacetime — his class of examples is consis- 
tent with our positivity results. We also present a simple 



explicit model which illustrates this point, and discuss 
implications for the existence of traversable wormholes. 

A key element in Visser's argument is a method for ob- 
taining approximate self consistent solutions to the semi- 
classical equations, starting from solutions of the classical 
Einstein equation. Essentially, the idea is an extension 
of our long- wavelength or small h expansion beyond the 
context of perturbation theory about flat spacetime. We 
now describe this approximation scheme in a language 
similar to that we have used in the body of the paper. Let 
^'ci be some classical field with stress tensor T^l [i>ci, 5cd]) 
and let ,g^°\ ^^^''(o) g, solution of the classical Einstein 
equation 



(El) 



Now consider an additional quantum field and let us 
seek solutions of the semiclassical equation 



KGab[gcd] = [*cl, <7cd] + {fab)[g 



'cd\- 



(E2) 



Now if the state of the quantum field is the incoming vac- 
uum state, then the second term on the right hand side 
above will be in order of magnitude ^ Lp/C^, whereas 
the first term should be ~ where £ is the length- 

scale determined by the classical background solution 
Sah- Therefore if C ^ Lp, the quantum stress tensor 
can be treated as a small perturbation, and a leading or- 
der approximation to the self consistent solution will be 
given by 



gab 



„(0) , „(i) 

9ab ^ 9ab 
^cl,(0) _^ ^cl,(l)^ 



(E3) 



where g^f^' and ^f^^'^^) are calculated from the linearized 
version of Eq. (E2) with source (7a(,)[p^°^] as well as the 
linearized field equation for At leading order, in 

regions of the spacetime where \[''^^'(°) and vj/'^''^^) vanish, 
the Einstein tensor for this self consistent solution will be 
just the test field quantum stress tensor on the classical 
background, as claimed by Visser [43] . 

Consider now the application of this scheme to the 
Boulware vacuum outside a Schwarschild black hole. 
This state is unphysical because the expected stress ten- 
sor diverges on the horizon. However, Visser argues that 
the expected stress tensor far from the horizon will likely 
be approximately the same as the stress tensor for the 
static vacuum state oiitside a spherical star with radius 
close to its Schwarschild radius. [The stress tensors will 
not be exactly the same because the expected stress ten- 
sor in the static vacuum state docs have a non-local de- 
pendence on the spacetime geometry.] Therefore, the 
approximate self consistent solution obtained from the 
above scheme starting from a compact star and a test 
field in the static vacuum state should violate ANEC over 
macroscopically large regions far from the star. 
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It might appear that this violation of ANEC is qual- 
itatively different from the cases treated in our analysis 
in the body of the paper, and, thus, that new possibil- 
ities are open for severe violations of ANEC. However, 
we shall now argue that this is not the case by show- 
ing that examples of this type also exist in the context 
of perturbation theory off of Minkowski spacetime, and 
that these examples satisfy our smeared version of ANEC 
provided only that stress-energy of the classical matter 
itself satisfies ANEC. To see this, consider perturbing 
Minkowski spacetime by adding a classical stress tensor 

-^ifc^ '^'- ^'-'^ example this could be the linearized stress 
tensor for a static spherical star. The first order metric 
perturbation in the semiclassical theory will then satisfy 
a modified version of Eq. (3.4) wherein the term T^l'^''^^ 
is added to the right hand side. When one goes to the 
long wavelength limit and performs a reduction of order, 
one finds that, to lowest order, the metric perturbation 
satisfies Eq. (4.20) with hsab replaced by T^l^''^\ [Note 
that the resulting equation simply says that k times the 
Einstein tensor of the semiclassical solution is just the 



sum of the classical stress tensor T 



(l),cl 



ab 



and the vac- 
uum polarization stress tensor in the spacetime r]ab + hfij, 
where g[jj^[/i'^'] = T^l^''^\ in agreement with our discus- 
sion above.] However, our proof of the positivity of the 
smeared ANEC integral given in Sec. VIE above used 
only the positivity of the ordinary Minkowski ANEC inte- 
gral for Sab- Thus, if T^b^' ''^ satisfies ANEC in Minkowski 
spacetime, our analysis shows that the smeared ANEC 
integral in the semiclassical spacetime will be positive to 
first order. 

A simple example will illustrate these points. Consider 

a point mass m moving along the geodesic x°'{t) = tu"" 
in Minkowski spacetime. Let /i^'^ be the linearized grav- 
itational field of this point mass, which is given by 



27r mUaUb ^, 



(E4) 



Next, using Eq. (3.32) we can calculate the expected 
stress tensor in the incoming vacuum state on the space- 
time with metric rjab + hf^,, which is given by 



2ai?UaUb 



-{2a/3 + 2b){kakb-r]abk''ke) 



(E5) 



From this stress tensor we can obtain a second linearized 
metric perturbation Finally, we calculate the 

ANEC integral along the geodesic a;°(A) = A"" + AA° of 

the Einstein tensor of the spacetime rjab + h^ab + /^ob'^"*'^™' 
which is given by 

/ = I d^k Gab{k)X'X'' e^'^^ '' 5{k ■ A). (E6) 

Using Eqs. (6.17), (D13) and (E5) we obtain for At ^ 



7 = - 



6ma 

(27r)4K 



(E7) 



where At is the length of the component of A" per- 
pendicular to A° and u", i.e., the impact parameter of 
the null geodesic with respect to the point mass. Equa- 
tion (E7) is essentially the large r, weak gravity limit of 
Visser's counterexample to ANEC involving a static star. 

Thus, the ANEC integral is negative along all geodesies 
away from the point mass in the self consistent solution. 
Nevertheless, this is consistent with our result of positiv- 
ity of smeared ANEC, since in calculating the smeared 
ANEC integral for a geodesic a distance At from the 
point mass, the negative contribution (E7) will be com- 
pensated by the positive contribution from the point 
masscis stress tensor itself, which also scales like A^* be- 
cause of our smearing function (1.12) [93]. 

This example illustrates that our results are consistent 
with having negative ANEC integrals over a macroscopic 
region. The price one must pay is that the amount of 
ANEC violation is restricted to be very small compared 
to distant mass scales. We now give a crude argument 
which suggests that the restriction is easily sufficient to 
prevent the existence of macroscopic traversable worm- 
holes. Let us characterize the region with negative ANEC 
integral by the quantity with dimensions of mass 



M_ = 



(E8) 



and suppose that there is a region of positive ANEC inte- 
gral a distance ~ A away characterized by the mass M+. 
Then, in order of magnitude, our result implies that 



M_ + 



m 



Consider now a static macroscopic wormhole, and sup- 
pose that the wormhole can be characterized by one 
lengthscale jC. On dimensional grounds, the energy den- 
sity required to hold it open should be be of order 
~ —1/C^. [This is confirmed by explicit calculations 
in specific examples by Ford and Roman [14].] Conse- 
quently we have M_ ~ — £, and therefore 

A 



Since the distance A to the positive mass region should 
be > C, we obtain M+ > 10^^° Mq{jC /Icm)^, which is 
a ridiculuously large mass. Moreover, the natural re- 
quirement that the A be larger than the gravitational 
radius Af_|_ of the positive mass yields the restriction 
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of T does not commute with making the change of vari- 
ables. This issue is not particular to the reduction of order 
proceedure; it arises whenever one has an equation that 
is known only up to first order in some small parameter. 
More precisely, the finite size corrections will be larger 
than this by the factor ~ C/{ct) > 1, c.f., Eq. (4.1) 
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This seemingly somewhat convoluted method of specify- 
ing the gauge is necessary to translate the specification 
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f 

Jo 



dx Jo{kx) 



kKi{k)/2, 



where Ki is a modified Bessel function {not the function 
Ki discussed in the body of the paper) , and by using the 
expansion 

kKi{k) = l+fc2log(fc)/2+[(27- l)/4-log(2)/2] k^+0{k^ 

M. Reed and B. Simon Functional Analysis (Academic 
Press, San Diego, 1980), p. 206. 

Strictly speaking, although the exponentially growing 
pieces of Eq. (A14) satisfy the equation of motion (3.38), 
the exponentially decaying pieces do not because for 
these pieces the integral over the past light cone in the 
non-local term in the equation of motion does not con- 
verge. However, from the time reversal invariance of the 
semiclassical theory it is apparant that tfie exponentially 
growing and exponentially decaying solutions are equally 
physical. If we had chosen to write the equations of mo- 
tion in terms of "out states" at J"^ instead of "in states" 
at J~ , then we would have obtained a first order semi- 
classical equation which posseses exponentially decaying 
solutions but no exponentially growing solutions. 
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[93] The formula (E7) is only valid for At / 0; there are 
additional distributional components of the ANEC inte- 
gral as a function of At at the origin, coming from both 
the classical and quantum metric perturbations, h'^f, and 
^quantum ^YiBu Calculating the smeared ANEC integral, 
it is necessary to include these components in order to 
obtain a finite result; a naive application of the formula 
(E7) yields an infinite answer. The classical contribution 
is m|A''wa|5''(AT). 
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TABLE I. A summary of our results for the ANEC integral in the different cases. "Pure to first order" indicates that the 
two point function of the scalar field is pure to first order in e. "Always > 0" means that the smeared ANEC integral is always 
strictly positive for all solutions of the equations except for the trivial, fiat spacetime/vacuum solution. 

Pure to first order Mixed at first order 

No smearing Smearing No smearing Smearing 

First order in e = =0 Can be < Always > 

in long wavelength limit 

Second order in e Can be < Always > N/A N/A 

in long wavelength limit 
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TABLE II. In this table, for the aid of the reader, we list in alphabetical order some of the symbols that appear in the paper. 
We do not list symbols whose meaning is very conventional, or symbols which are used only in the immediate vicinity of where 
they are introduced. For each item listed, wc give a brief description, and also a reference to the equation in the text where 
the symbol first appears, or after which the symbol is first introduced. Except in special cases, we do not list separately the 
following variants of symbols: symbols with tildes or symbols with the superscripts (0), (1) or (2). The former always denote a 
Fourier transform, and the latter superscripts always denote expansion coefficients in an expansion of a quantity in powers of 
e. 



Symbol 


Meaning 


Equation in 
first appears 


which Other 
Equati 


a 


Coefficient that controls the anomalous scaling 


(2.10 ) 


(3.34 ) 


A 


Expansion coefficient in long wavelength expansion 


(6.42 ) 




Aab 


Fourth rank local curvature tensor 


(2.7 ) 


(2.10 ) 




Linearized local curvature tensor for metric perturbation h^^fj 


(3.33 ) 


(3.32 ) 




Operator that acts on metric perturbations to yield a local curvature tensor 


(3.45 ) 


(3.46 ) 




Local curvature tensor at zeroth order in expansion 


(4.20 ) 


(3.46 ) 


b 


Coefficient that controls the anomalous scaling 


(2.10 ) 


(3.34 ) 


B 


Expansion coefficient in long wavelength expansion 


(6.42 ) 




Bab 


Fourth rank local curvature tensor 


(2.8 ) 


(2.10 ) 




Linearized local curvature tensor for metric perturbation h^^ij 


(3.33 ) 


(3.32 ) 




Operator that acts on metric perturbations to yield a local curvature tensor 


(3.45 ) 


(3.46 ) 


Bab 


Local curvature tensor at zeroth order in expansion 


(4.20 ) 


(3.46 ) 


c 


Expansion coefficient in long wavelength expansion 


(6.42 ) 






D'Alembertian type wave operator (wrt x or y) 


(2.17 ) 






Operator arising in point splitting prescription, consisting of 


(2.19 ) 


(7.15 ) 




differential operator followed by (implicit) coincidence limit 






First order change in Vab due to metric perturbation h^^^ 


(3.1) 






Operator that acts on metric perturbations giving first order change in Dab 


(3.1) 






Operator taking sources Sx and Sy to bisolution of wave equation 


(3.13 ) 




Fin 


Regularized two point function of incoming state Win. 


(2.24 ) 




Fin 


Regularized two point function of rescaled incoming state Win. 


(3.41 ) 


(3.42 ) 


Fi(fc), F2{k), Fz(k) Functions arising in exact solutions 


(A6) - (AS) 


(A5 ) 


/ 


"Pure frequency" part of two point function, in momentum space 


(6.1 ) 


(C2 ) 


/ 


Same as f^^^ but multiplied by lk| |k'| 


(6.15 ) 


(6.14 ) 




Linearized Einstein tensor for metric perturbation 


(3.6 ) 






Operator that acts on metric perturbation to give linearized Einstein tensor 


(3.6 ) 




'^ab 


Second order part of Einstein tensor, depending on h^^^ and h^ab ■ 


(5.12 ) 




^ab 


Rescaled version of second order part of Einstein tensor 


(5.12 ) 






Operator that acts on pairs of metric perturbations 


(3.6 ) 




9 


"Mixed frequency" part of two point function, in momentum space 


(6.1) 


(C2) 


9 


Same as g''^^ but multiplied by |k| |k'| 


(6.15 ) 


(6.14 ) 


G 


Two point bidistribution of the state w on {M,gab)- 


(2.13 ) 




Go 


Two point bidistribution of the in vacuum state lvq on (M, gab)- 


(2.20 ) 




Gin 


Two point bidistribution of the in state Win on {M,riab)- 


(2.20 ) 




Gin.O 


Two point bidistribution of the vacuum state Win.o on {M,r]ab)- 


(2.20 ) 




"■ab ' "-ab 


Metric perturbations 


(3.1) 




"'ab ' "ab 


Rescaled metric perturbations in treatment of long-wavelength limit 


(3.44 ) 


(7.2) 




Horowitz distribution that enters the expression for vacuum polarization 


(3.35 ) 


(3.32 ) 


n 


Usual Fock space of states on Minkowski spacetime 


(3.5 ) 






ANEC integral and its expansion coefficients 


(5.8 ) 




r 7-(l) 7-(2) 

J-S'i J-s -1 -'- s 


Generalized ANEC integral and its expansion coefficients 


(5.2 ) 




f r{l) f(2) 


Limiting form of generalized ANEC integral, with only transverse smearing 


(5.4 ) 




/f(xt) 


ANEC integral without backreaction as a function of transverse displacement 


(6.33 ) 


(6.34 ) 


7f(xt) 


Same as If but for rescaled incoming state Win 


(6.38 ) 




J{k,k') 


Function that arises in solutions of reduced order semiclassical equations 


(6.24 ) 


(6.27 ) 


k, r,kT 


Coordinates on momentum space 


(6.12 ) 






Distribution that describes the effect of backreaction on the ANEC integral 
for solutions of reduced order equations 


(6.31 ) 


(6.35 ) 


Ki{yiT), ^i(kT) 


Function that describes the effect of backreaction on the ANEC integral 


(Bl) 


(B2) 
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for solutions of original scmiclassical equations 





Fixed, Planck scale frequencies controlling tachyon type solutions 


(A15 ) 


A16 ) 


L 


Lcngthscale of incoming state or of semiclassical solution 


(2.6 ) 




N 


Number of scalar fields in the 1 /N expansion 


(2.6 ) 




Sab 


Source tensor in linearized semiclassical equation 


(3.39 ) 


(IIB) 


5(xt) 


Transverse smearing function 


(5.1) 


(5.2) 


'-'061 >3 


Comcidence limits of derivatives of two point function 


(7.16), (7.17) 




Tab bed] 


Linear map on states on (M,riab) returning conserved stress tensor on {M,gab) 


(3.2 ) 




rp(0) 
^ ab 


Usual Minkowski spacetimc stress tensor (linear map on states) 


(3.3 ) 




^ab 1 J' ^ab L ' J 


Expansion coefficients of Tat [fled (e)] that act on metric perturbations 
or pairs of metric perturbations and return linear maps on Minkowski 
spacetime states 


(3.3 ) 




XT 


Two dimensional transverse coordinate on Minkowski spacetime 


(5.1 ) 


(6.12) 


Zab 


1 ' 1 ' C C 1 1 11 1 A A A 1 A J 'A 

Uombmation of fourth rank local curvature tensors that enters into 
anomalous scaling of stress tensor 


(3.57 ) 


(2.10 ) 




Linearized form of Zab for metric perturbation ft^^^ 


(3.46 ) 


(3.32), 


Operator that acts on metric perturbations to yield anomalous scaling tensor 


(3.46 ) 




a 


Parameter used to characterize long wavelength limit 


(3.41 ) 




a 


Unknown numerical coefficient of local curvature term in stress tensor 


(2.9 ) 


(3.32) 


ao 
P 


Real part of complex frequency in oscillatory modes of exact solutions 
Coordinate on momentum space 


(All ) 
(6.12 ) 




/3 


Unknown numerical coefficient of local curvature term in stress tensor 


(2.9 ) 


(3.32) 


/3o 


Imaginary part of complex frequency in oscillatory modes of exact solutions 


(All ) 




7 


Null geodesic in spacetime (M,gab), also zeroth order geodesic in {M,riab) 


(5.2 ) 




7 


Coordinate on momentum space 


(6.12 ) 




e 


Basic expansion parameter of perturbation expansion 


(1.7) 


(1.6) 


e 


Generic small parameter in discussion of higher order time derivative 
equations of motion 


Sec. IV C 




c 


Null coordinate on Minkowski spacetime, also a Fermi- Walker coordinate 
on (M, gab) in a neighborhood of the null geodesic 


(5.1) 


(5.15) 


e 


Step function 


(6.17 ) 






Function entering definition of generalized ANEC integral 


(5.1) 


(5.15) 




Inverse of Newtons constant 


(2.1) 




A 


AfRne parameter along geodesic 7, also Fermi- Walker coordinate on (M, gab) 
in a neighborhood of 7, also null coordinate on Minkowski spacetime 


(5.1) 


(6.12) 


A 


Undetermined lengthscale appearing in expression for linearized stress tensor 


(3.32 ) 


(6.17) 


A 


Rescaled version of above lengthscale A 


(6.31 ) 


(A9) 


Acrit 


Critical value of lengthscale A 


(AlO ) 




A" 


Tangent vector to geodesic 7, also vector field on Minkowski spacetime 


(5.1) 


(6.12) 


A, At, Ai 


Lengthscales entering definition of generalized ANEC integral 


(5.1 ) 


(5.2) 


A 


Used m Appendix a instead of At. 


(B3 ) 




AT,crit 


Critical value 01 transverse smearing lengthscale 


(6.45 ) 






Numerical coefficient in expansion of Fourier space smearing function 


(6.41 ) 






Curvature coupling coefficient 


(2.1 ) 






Expansion coefficients of incoming state for Fock space states 


(3.7 ) 




^(1) ^(2) 
P , P 


Expansion coemcients 01 rescaled incoming state tor i^ock space states 


(7.9), (7.10) 




T 


Timescale characterizing radiation reaction effects 


(4.2 ) 




T 


("ieneric small parameter in discussion of reduction of order 






* 


Evolution timescale for solutions of radiation reaction equation 






^+5 ^— 


Po^sitivp flTiH Tipcrativp rrpmiPTicv riipcp<5 of fiplH onpratrtr 


(3.21 ) 






Rescaled versions of metric perturbations , h^^jj 


( p;n ^ 
l^o.ou ) 


(7 ')^ 


(1,0) (2,0) 


Leading order approximations to xiV ' xit' ^ong wavelength expansion 


(3.52 ) 


(3.58), 




Vacuum state on Minkowski spacetime 


(1.6) 




(1) (2) 
'^in = '-'in 


Expansion coefficients of incoming, Minkowski spacetime state 


(1.6) 




-(1) -(2) 
^in ' ^in 


Expansion coefficients of rescaled, incoming, Minkowski spacetime state 


(3.42 ) 






Critical frequency of the order of the Planck length 


(6.31 ) 


(A9) 


{;■} 


Product notation for stress tensors acting on Minkowski spacetime states 


(3.2 ) 
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